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ABSTRACT 


The  five-year  research  program  sponsored  by  this  contract  has 
produced  technical  findings  in  several  areas:  characterization  of  semi¬ 
conductor  devices  and  integrated  circuits;  electronic  ceramics;  noise 
within  solids  and  semiconductor  devices;  quantum  and  electronic  properties 
of  elemental  and  compound  semiconductors;  and  magnetic  films.  This  report 
begins  with  a  summary  of  these  findings.  It  then  describes  our  most  re¬ 
cent  results,  as  follows.  It  is  shown  how  the  falloff  of  transistor  gain 
and  cutoff  frequency  at  high  currents  depends  on  the  device  operating 
conditions;  this  description  enables  a  new  interpretation  to  be  given  to 
experimental  results  previously  reported.  The  decline  observed  in  the 
response  of  p-i-n  photodiodes  for  short  wavelength  excitation  is  explained 
by  accounting  for  carrier  transport  and  recombination  in  the  p-type  region. 
We  describe  a  method  for  the  numerical  calculation  of  the  energy  structure 
and  cuarge  distribution  resulting  from  various  defects  in  the  solid  state 
including  clustering  and  amorphous  disorder.  A  recombination  and  trapping 
model  is  developed  to  interpret  photoconductivity  data  observed  on  high 
resistivity  3As2Se3~2Sb2Se3  amorphous  films.  A  study  of  the  influence  of 
crystalline  microstructure  on  the  electrical  characteristics  of  Li20-2Si0 
crystals  shows  that  for  volume  fractions  of  crystals  exceeding  twenty 
percent  the  dc  conductivity  and  ac  dielectric  losses  are  markedly  reduced. 


SUMMARY 


The  five-year  research  program  sponsored  by  this  contract  has 
produced  technical  findings  in  several  areas:  characterization  of  semi¬ 
conductor  devices  and  integrated  circuits;  electronic  ceramics;  noise 
within  sol  .ds  and  semiconductor  devices;  quantum  and  electronic  properties 
of  elemental  and  compound  semiconductors;  and  magnetic  films.  This  report 
begins  with  a  summary  of  these  findings.  It  then  describes  our  most  re¬ 
cent  results,  as  follows. 

For  transistors  having  epitaxial  collectors,  the  falloff  of  tran¬ 
sistor  gain  and  cutoff  frequency  at  high  current  may  depend  predominantly 
on  phenomena  occurring  in  the  collector  region.  We  show  in  detail  how 
the  mechanism  responsible  for  falloff  depends  on  device  operating  condi¬ 
tions  and  the  structural  and  material  parameters  of  the  device  make-up. 

We  illustrate  the  usefulness  of  this  description  in  understanding  observed 
device  behavior  and  show  how  it  enables  a  new  interpretation  to  be  given 
to  experimental  results  previously  reported.  Important  implications  for 
transistor  modeling  are  also  discussed. 

An  expression  for  the  quantum  yield  of  p-i-n  photodiodes  is  formulated 
that  takes  into  account  the  effects  of  carrier  transport  and  recombination 
in  the  p-type  region  and  consequently  describes  the  decline  in  the  res¬ 
ponse  observed  for  short  wavelength  excitation.  Applied  to  a  silicon 
p-i-n  photodiode,  this  expression  yields  theoretical  predictions  in  ex¬ 
cellent  agreement  with  experimental  data  over  the  entire  spectral  range. 

Many  of  the  properties  that  make  semiconductor  devices  useful  in  both 
electronic  and  optical  applications  arise  from  the  defect  structure  in  the 
solid  state.  Previous  methods  available  for  calculating  the  details  of 
the  energy  levels  and  charge  distribution  resulting  from  defects  apply 
only  to  such  simple  defect  centers  as  shallow-level  donors  and  acceptor 
impurities.  We  describe  here  a  method  enabling  calculations  for  various 
other  kinds  of  defects  including  vacancies,  vacancy- impurity  clustering 
and  amorphous  disorder.  Using  this  method  we  have  done  calculations  for 
certain  specific  defects  in  silicon:  boron  and  phosphorous  impurities, 
the  vacancy  and  boron  and  phosphorous  complexes,  and  impurity  clusters. 

The  pertinence  to  radiation-damaged  silicon  is  discussed. 


We  report  the  results  of  photoconductivity  experiments  done 
on  high- resistivity  3As2Se2_ 2Sb2Se^  amorphous  films  in  a  temperature 
range  from  278°K  to  308°K  as  functions  of  the  light  intensity. 

At  intermediate  light  levels  the  photocurrent  varies  as  the  light 
intensity  to  the  0.7  power  and  at  high  light  intensities  it  varies 
in  direct  proportion  to  the  light  intensity.  To  explain  the 
observed  data  a  recombination  trapping  model  is  suggested. 

Previous  studies  in  this  program  have  shown  that  it  is  possible 
to  control  the  volume  fraction  of  crystals  in  pure  Li20-2Si02  glasses 
by  nucleation  and  growth  heat  treatments.  The  resulting  glass- 
ceramics  exhibit  improved  mechanical  and  thermo  properties.  The 
present  study  shows  the  influence  of  the  crystal  in  microstructure 
on  the  electrical  characteristics  of  the  glass-ceramics.  If  the 
volume  fraction  of  crystals  exceeds  twenty  percent  the  dc  conductivity 
and  ac  dielectric  losses  of  the  material  decline  markedly.  The 
temperature  dependence  of  the  electrical  processes  increases  with 
the  volume  fraction  of  crystals.  Interpretation  of  these  results 
is  given  based  on  an  application  of  heterogeneous  dielectric  mixing 
laws  to  glass-ceramic  materials. 
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I .  Introduction 

The  original  general  objective  of  this  research  program  was  to 
establish  at  the  University  of  Florida  a  "Center  of  Competence  in 
Solid  State  Materials  and  Devices"  capable  of  doing  research  of  benefit 
to  the  Department  of  Defense.  At  this  stage,  five  years  following 
initiation  of  contract  support,  ample  evidence  exists  that  the 
University  of  Florida  has  indeed  become  such  a  center  of  competence. 
Among  the  evidence  supporting  this  conclusion  are  the  technical 
findings  fostered  by  this  contract.  These  are  documented  in  the  nine 
previous  reports  in  this  series1-9  and  in  118  related  papers  published 
in  research  journals. 

During  the  past  year  the  objectives  of  the  research  program  have 
shifted  decisively.  In  accordance  with  an  agreement  made  with  the 
contract  monitor  in  January  1971,  we  have  focussed  the  objectives  to 
give  greatest  emphasis  to  the  following  problem  areas: 

1.  The  development  of  new  types  of  detectors  with  stress 
given  to  utilizing  the  Dember  and  photomagnetoelectric 
effects  in  semiconductors  and  to  utilizing  new  mask  geometries 
for  Schottky  barrier  detectors  compatible  with  conventional 
silicon  technology. 

2.  The  incorporation  of  the  effects  of  irradiation  into 
the  characterization  of  device  behavior;  junction  field- 
effect  devices  are  to  receive  special  attention.  Emphasis 
is  placed  on  the  development  of  characterizations  enabling 
computer  simulation  and  design  of  circuits. 

3.  Full  static  and  dynamic  characterization  of  bipolar  and 
field-effect  devices  aiming  both  toward  improved  device 
design  and  toward  improved  circuit  implementation.  The 
characterization  is  to  deal  with  operation  over  a  wide  range 
of  temperatures,  with  emphasis  on  the  range  between  ordinary 
room  temperature  and  cryogenic  temperatures. 

4.  Noise  studies  proceeding  collaterally  with  each  of  the 
efforts  listed  above. 

5.  Radiation  studies  of  amorphous  semiconductors. 


The  content  of  the  present  report  reflects  this  change  in  emphasis. 
The  report  begins  with  a  summary  of  major  technical  findings  produced 
during  the  five  years  of  this  contract,  and  then  describes  our  most  recent 
results. 
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11 •  Summary  of  Research  Program  (September  1967  to  September  1972) 

The  five-year  research  program  sponsored  by  this  contract  has  produced 
technical  findings  in  several  areas:  characterization  of  semiconductor  devices 
and  integrated  circuits:  electronic  ceramics;  noise  within  solids  and  semi¬ 
conductor  devices;  quantum  and  electronic  properties  of  elemental  and 
compound  semiconductors;  and  magnetic  films.  These  findings  appear  documented 
in  over  one  hundred  technical  papers  published  in  research  journals  and  in 
a  comparable  number  of  papers  presented  orally  at  national  and  international 
research  meetings.  Here,  following,  we  describe  the  major  results. 


A.  CHARACTERIZATION  OF  SEMICONDUCTOR  DEVICES  AND  INTEGRATED  CIRCUITS 
A-l  Introduction 

Semiconductor  devices  and  circuits  form  the  core  of  the  great  majority 
of  electronic  systems.  At  present,  a  substantial  part  of  the  design  of 
these  devices  and  circuits  is  based  on  laboratory  experimentation  and  iterative 
cut-and-try  procedures.  To  reduce  the  cost  and  time  of  such  procedures  and 
to  secure  superior  system  performance,  much  effort  in  this  nation  now  aims 
toward  the  development  of  computer  techniques  to  aid  design.  Essential  to 
the  success  of  such  techniques  is  the  proper  characterization  of  semiconductor 
devices  and  integrated  circuits.  As  its  main  objective,  this  portion  of  our 
research  program  has  focussed  attention  upon  fundamental  problems  related 
to  securing  these  characterizations. 

One  can  view  the  starting  point  in  semiconductor  circuit  and  system 
design  as  the  processing  steps  used  in  fabrication:  mask  generation, 
diffusion  or  implantation  of  impurities,  etc.  The  end  point  is  a  circuit 
that  meets  specifications  in  system  performance  and  cost.  In  between  are 
four  modeling  problems,  as  Fig.  IJ-l  Illustrates. 

The  first  of  these  problems  seeks  to  describe  for  each  device  in  the 
circuit  the  relationship  between  the  processing  and  the  resulting  physical 
make-up:  the  spatial  distribution  of  impurities  and  defects,  the  corresponding 
distribution  of  the  pertinent  material  parameters,  etc. 


3 


4 


Fig.  II-l  Four  modeling  problems  involved  in  modern  circuit  design. 


In  the  second  problem  one  deals  with  the  relevant  partial  differential 
equations,  often  by  numerical  methods,  to  relate  the  physical  make-up  to 
the  spatial  and  temporal  dependence  of  such  internal  variables  as  free 
carrier  densities,  charge  densities,  current  density,  and  electric  field. 
Response  is  sought  for  excitation  pertinent  to  the  proposed  system  application. 
The  understanding  thus  gained  directly  guides  the  design  of  the  devices 
comprising  the  circuit. 

Moreover,  it  supplies  insight  vital  in  the  third  modeling  problem, 
in  which  one  links  the  physical  make-up  and  the  current-voltage  behavior  at 
the  terminals  of  each  device.  The  vehicle  describing  this  behavior  is 
usually  called  a  circuit  model  or  an  equivalent  circuit.  The  availability 
of  adequate  circuit  models  is  essential  in  the  design  of  any  circuit. 

The  fourth  modeling  problem  involves  selection  from  the  reservoir  of 
available  models  of  a  circuit  model  to  represent  each  device  in  the  circuit 
being  designed. 

This  research  program  has  made  contributions  mainly  in  the  last  three 
of  these  modeling  problems,  as  is  described  below. 

A-2  Relating  Device  Make-up  to  Internal  Device  Variables 

In  this  second  modeling  problem  we  have  contributed  in  several  areas. 

A-2.1  THE  EFFECTS  OF  HIGH  DOPING  IN  SEMICONDUCTOR  DEVICES:  Although  the 
fabrication  techniques  used  for  two  decades  have  produced  regions  of  semi¬ 
conductor  devices  such  as  bipolar  and  field-effect  transistors  in  which 
high  concentrations  of  impurities  (exceeding  5  x  1017  per  cubic  centimeter) 
are  present,  traditional  semiconductor  device  theory  and  design  practice 
largely  ignore  the  consequent  effects  on  device  operation.  Aside  from  the 

first-order  phenomenological  treatment  by  Kauffman  and  Bergh  and  by  Buhanan, 
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and  the  early  similar  work  by  Kroemer,  transistor  theory  took  no  account 

of  the  presence  of  heavy  doping.  Urder  the  sponsorship  of  this  contract 
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Kleppinger  and  Lindholm  ’  pior  ered  an  effort  to  describe  the  effects  of 

high  doping  in  a  way  useful  to  the  semiconductor  device  designer.  Based  on 

this  contract-supported  work,  which  is  described  below,  Van  Overstraeten  et  al.,® 

DeMan  and  Mertens7(of  the  University  of  Louvain  in  Belgium)  and  Mock®  (of 

the  IBM  Corporation)  have  demonstrated  that  heavy  doping  takes  a  dominant 

role  in  determining  the  current  gain  and  frequency  response  of  bipolar 
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transistors.  Before  the  availability  of  these  results,  efforts  toward 
optimizing  transistor  design  were  frustrated  by  the  poor  capability  of 
existing  theory  to  predict  these  two  parameters  from  knowledge  of  the 
physical  make-up  yielded  by  the  fabrication  processing.  The  new  results, 
based  on  contract-supported  research,  enabled  greatly  improved  prediction, 
thus  providing  a  basis  for  improving  the  design  of  bipolar  transistors  and 
the  electronic  systems  in  which  they  are  used.  Recently  DeMan  and  Mertens7 
report  the  incorporation  of  the  effects  of  heavy  doping  into  their  computet 
program  for  the  computer-determined  design  of  semiconductor  circuits. 

The  main  objective  of  the  research  by  Kleppinger  and  Lindholmf’^ 
supported  by  this  contract,  was  to  combine  and  phrase  existing  theories 
available  in  the  physics  literature  in  a  form  amenable  to  the  semiconductor 
device  theorist  and  designer.  In  the  first^  of  the  two  papers  published 
from  this  research,  Kleppinger  and  Lindholm  developed  and  justified  an 
approach  to  the  problem  of  describing  the  dependence  on  impurity  concentra¬ 
tion  of  the  density  of  quantum  states  in  a  semiconductor.  The  approach 
applies  generally  to  any  semiconductor,  including  those  widely  used  in 
device  fabrication:  silicon,  germanium  and  gallium  arsenide.  The  approach 
combined  previously  separate  theories  concerning  the  various  contributors 
to  the  dependence  of  the  quantum-state  density,  utilizing  an  approximation 
that  enabled  application  of  the  superposed  theory  over  the  range  of  impurity 
concentrations  of  technological  Interest.  To  demonstrate  the  practical  utility 
of  the  approach,  the  prediction  yielded  by  its  application  to  gallium  arsenide 
was  compared  to  an  experimentally  determined  density  of  states.  Agreement 
was  satisfactory.  Throughout,  the  paper  stressed  the  implications  to  semi¬ 
conductor  device  design. 

In  their  second  paper, Kleppinger  and  Lindholm  applied  the  approach 
explicitly  to  silicon,  which  is  the  semiconductor  of  greatest  importance  in 
the  technology  of  integrated  circuits.  They  calculated  the  density  of  states 
and  the  Fermi  level,  and  showed  the  calculated  density  states  to  agree  with 
pertinent  experimental  results. 

As  a  consequence  of  this  research,  collaboration  began  between  Dr. 

Lindholm  and  Dr.  Roger  Van  Ovcrstraeten,  Professor  and  Head  of  the  Solid 
State  Electronic  Laboratory,  University  of  Louvain,  Belgium.  In  May  of  1972 
Dr.  Van  Overstraeten  began  a  month-long  visit  in  Gainesville,  Florida  to  do 
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research  with  Dr.  Lindhoim.  This  visit  produced  a  research  paper9  which 
places  on  firmer  ground  the  description  of  current  within  heavily  doped 
material.  The  paper  includes  the  derivation  of  an  expression  describing 
current  in  a  semiconductor  with  an  arbitrary  quantum  density  states, 
enabling,  as  a  special  case,  a  description  of  the  effects  of  the  warping 
•of  the  quantum  density  of  states  which  accompanies  heavy  doping.  Special 
attention  is  given  to  the  Einstein  ratio  of  diffusivity  to  mobility. 

In  semiconductor  device  theory,  this  ratio  has  gained  importance  for  two 
reasons.  First,  it  appears  whenever  carrier  drift  and  diffusion  coexist, 
which  in  device  materials  frequently  occurs.  Second,  mobility  is  much  more 
easily  measured  than  is  diffusivity,  and  knowledge  of  the  Einstein 
ratio  thus  yields  a  value  for  diffusivity  without  need  for  direct  measure¬ 
ment.  For  nondegenerate  semiconductors  in  thermal  equilibrium,  the  Einstein 
ratio  equals  the  thermal  voltage.  However,  rigorous  analysis  of  device 
behavior  requires  understanding  of  the  Einstein  ratio  in  semiconductors  that 
are  both  highly-doped  and  in  nonequilibrium.  Early  efforts  supported  by 
this  contract  *  contributed  to  this  understanding,  but  ignored  the 
distortion  of  the  quantum  density  of  states  that  occurs  with  high  doping. 

The  computations  done  in  our  more  recent  study9  show  a  dependence  of  the 
Einstein  ratio  on  impurity  concentration  that  contrasts  markedly  with 
earlier  predictions. 

REFERENCES  Vindicates  research  sponsored  by 
this  contract) 

1.  W.  L.  Kauffman  and  A.  A.  Bergh,  "The  Temperature  Dependence  of  Ideal 

Gain  in  Double  Diffused  Silicon  Transistors,"  IEEE  Transactions  on 
Electron  Devices.  ED-15,  732  (1968).  ' 

2.  D.  Buhanan,  "Investigation  of  Current-Gain  Temperature  Dependence  in 
Silicon  Transistors,"  IEEE  Transactions  on  Electron  Devices.  ED-16. 

117  (1969).  - - —  - 

3.  H.  Kroemer,  "Theory  of  a  Wide-Gap  Emitter  in  Transistors,"  Proceedings 

I.R.E.,  45,  1535  (1957).  - 

*  4.  D.  D.  Kleppinger  and  F.  A.  Lindhoim,  "Impurity  Concentration  Dependence 

of  the  Density  of  States  in  Semiconductors,"  Solid  State  Electronics. 

14,  199  (1971).  - 

*  5.  D.  D.  Kleppinger  and  F.  A.  Lindhoim,  "Impurity  Concentration  Dependent 

Density  of  States  and  Resulting  Fermi  Level  for  Silicon,"  Solid  State 
Electronics.  14,  407  (1971). 
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6.  R.  Van  Overstraeten,  H.  De  Man,  and  R.  Mertens, "Transport  Equations 
in  Heavy  Doped  Silicon,"  IEEE  Transactions  on  Electron  Devices 
(March  1973). 

R.  Van  Overstraeten,  H.  DeMan,  and  R.  Mertens,  "Calculation  of  the 
Emitter  Efficiency  of  a  Bipolar  Transistor,"  IEEE  Transactions  on 
Electron  Devices  (to  appear). 

7.  H.  DeMan  and  R.  Mertens,  "SITCAP — A  Simulator  of  Bipolar  Transistors 
for  Computer-Aided  Circuit  Analysis  Program,"  Digest  of  Technical 
Papers:  1973,  IEEE  International  Solid-State  Circuits  Conference,  16, 
104  (1973). 

8.  M.  S.  Mock,  "Transport  Equations  in  Heavily  Doped  Silicon,  and  the 
Current-Gain  of  a  Bipolar  Transistor,"  Solid  State  Electronics 

(to  appear) . 

*9.  R.  Van  Overstraeten  and  F.  A.  Lindholm,  "Generalized  Non-Equilibrium 
Transport  Equations  and  Einstein  Ratio  for  a  Semiconductor  with  an 
Arbitrary  Density  of  States,"  Solid  State  Electronics  (to  appear). 

*10  F.  A.  Lindholm  and  R.  W.  Ayers,  "Generalized  Einstein  Relation  for 
Degenerate  Semiconductors,"  Proceedings  of  the  IEEE,  56,  371  (1968). 

*11  S.  S.  Li  and  F.  A.  Lindholm,  "Alternative  Formulation  of  Generalized 
Einstein  Relation  for  Degenerate  Semiconductors,"  Proceedings  of  the 
IEEE,  56,  1256  (1968). 

A-2.2  TRANSPORT  IN  SEMICONDUCTORS  WITH  LOW  SCATTERING  RATE  AND  AT  HIGH 
FREQUENCIES :  The  conventional  hydrodynamic  transport  equations  used  in 

the  Analysis  of  the  electrical  characteristics  of  semiconductor  devices 
require  modification  when  the  number  of  recombination-generation-trapping¬ 
tunneling  events  become  comparable  to  the  number  of  scattering  events. 

These  modifications  may  become  important  in  describing  the  transport  of 
minority  carriers  across  the  thin  base  region  of  a  microwave  transistor 
and  in  devices  made  of  compound  semiconductors  with  high  mobilities  and 
low  lifetimes.  If  the  signal  frequency  becomes  comparable  to  the  collision 
frequency,  additional  modifications  may  be  needed. 

This  study  gives  a  rigorous  derivation  of  the  hydrodynamic  carrier 
transport  equations  by  extending  the  moment  method  to  include  the  moment 
equation  of  the  collision-free  path,  which  is  generalized  to  contain  both 
the  scattering  and  the  generation-recombination-trapping-tunneling  events. 
This  leads  to  a  description  of  the  flow  of  electrons  or  holes  that  contains 
an  inductive  term  neglected  in  previous  characterizations  in  semiconductor 
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device  physics.  If  the  signal  frequency  becomes  comparable  to  the 
reciprocal  of  the  average  collision  relaxation  time,  this  inductive  term 
can  take  an  important  role. 

REFERENCES 

1.  C.  T.  Sah  and  F.  A.  Lindholm,  "Transport  in  Semiconductors  with  Low 
Scattering  Rate  and  at  High  Frequencies,"  Solid  State  Electronics 
(to  appear) .  ™ 

A-2.3  THE  EFFECT  OF  A  BURIED  LAYER  ON  THE  COLLECTOR  BREAKDOWN  VOLTAGES 
OF  BIPOLAR  JUNCTION  TRANSISTORS;1  Bipolar  transistors  designed  for  high- 
frequency  or  high-speed  applications  are  very  shallow  structures.  For 
example,  devices  used  for  small-signal  microwave  amplification  may  have 

O 

emitter  and  collector  junctions  that  lie  within  5000  AngstrOm  units  from 
the  surface;  typically,  the  devices  are  fabricated  by  successive  diffusions 
into  N/N+  epitaxially  grown  silicon.  The  N-type  epitaxial  material  is  made 
thin,  less  than  10  microns,  to  reduce  collector  series  resistance.  Similar, 
devices  are  used  commonly  in  high-speed  digital  integrated  circuits.  The 
standard  buried  collector  technology  presently  employed  by  the  Bell  System 
uses  epitaxial  layers  which  are  7  to  9  microns  thick. 

In  this  study,  we  describe  a  premature  punch-through  effect  occurring 
in  this  type  of  transistor  that  can  lower  the  collector  breakdown  voltages 
below  the  values  predicted  by  existing  theory.  We  express  the  premature 
punch-through  voltage,  BVppT>  which  characterizes  the  effect,  in  terms  of 
parameters  of  the  fabrication  processing,  and  from  this  derive  an  expression 
for  the  thickness  of  the  epitaxial  collector  that  maximizes  BVppT  and,  hence, 
the  usable  region  of  device  operation.  Measurements  are  presented  that 
support  the  theoretical  predictions. 

REFERENCES 

1.  F.  W.  Hewlett,  Jr.,  F.  A.  Lindholm,  and  A.  J.  Brodersen,  "The  Effect 

of  a  Buried  Layer  on  the  Collector  Breakdown  Voltages  of  Bipolar  Junction 
Transistors,"  Solid  State  Electronics  (to  appear). 

A-2.4  HIGH  CURRENT  REGIMES  IN  TRANSISTOR  COLLECTOR  REGIONS;1  This  work  is 
described  in  Section  III  of  this  report. 

REFERENCES 

1.  D.  L.  Bowler  and  F.  A.  Lindholm,  "High  Current  Regimes  in  Transistor 
Collector  Regions,"  IEEE  Transactions  on  Electron  Devices  (March  1973). 
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A-3  Generation  and  Unification  of  Device  Circuit  Models 

This  phase  of  our  work  deals  with  the  third  modeling  problem  defined 
in  Figure  II-l.  The  emergence  in  recent  years  of  the  computer  as  a  tool 
for  the  analysis  and  design  of  semiconductor  devices  and  circuits  has 
focussed  renewed  attention  on  the  models  that  characterized  the  terminal 
behavior  of  devices.  In  particular,  because  the  computer  has  afforded 
for  the  first  time  a  general  means  of  accomplishing  the  detailed  analysis 
of  large-signal  transient  behavior,  it  has  awakened  special  interest  in 
large-signal  (nonlinear)  dynamic  models.  Our  work  has  focussed  mainly  on 
the  large-signal  dynamic  circuit  models  for  the  bipolar  transistor  and 
for  various  types  of  field  effect  transistors:  insulated-gate  (IGFET  or 
MOSFET) ,  junction-gate  (JFET) ,  and  metal-gate  (MESFET) . 

In  the  historical  perspective  of  the  research  done  under  the  sponsor¬ 
ship  of  this  contract,  there  were  two  main  issues  that  needed  initial  study: 

(a)  generation  of  new  and  more  completely  adequate  circuit  models 
using  advances  made  in  the  second  modeling  problem  defined  in 
Fig.  II-l,  concerning  semiconductor  device  physics;  and 

(b)  evaluation  of  existing  circuit  programs  for  accuracy  and 
suitability  in  circuit  analysis  codes,  making  modifications 
where  needed. 

Concerning  the  first  of  these  issues,  each  of  our  advances  outlined 
above  in  Section  A-2  potentially  enables  the  generation  of  new  and  superior 
circuit  models.  To  the  present,  our  studies  concerning  highly-doped  semi¬ 
conductors,  summarized  above,  have  had  the  most  impact.  The  effects  of  the 
warping  of  the  quantum  density  of  states  that  accompanies  high  doping  has 
recently  been  included  in  a  new  and  comprehensive  circuit— analysis  program 
developed  by  DeMan  and  Mertens  (reference  7,  Section  A-2). 

For  the  second  of  these  issues  our  work  toward  evaluating  the  existing 
circuit  models  for  semiconductor  devices  began  with  an  invited  research 
paper  at  the  1968  International  Symposium  on  Circuit  Theory.  This  paper 
tried  to  assess  where  the  technical  community  stood  in  1968  in  relation  to 
device  models  for  computer-aided  design.  Subsequently,  it  led  to  the 
following  more  comprehensive  statement  of  the  considerations  involved. 

A-3.1  DEVICE  MODELING  FOR  COMPUTER-AIDED  DESIGN  AND  ANALYSTS  OF  TNTFCRATFn 

2  '  — — - 

CIRCUITS:  A  list  is  proposed  of  properties  considered  to  be  desirable  for 
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models  for  integrated  devices.  These  properties  serve  as  a  guide  in 
modeling  per  se  and  as  vehicle  for  assessing  the  adequacy  of  any  given 
model.  In  illustration,  the  adequacy  is  examined  of  a  generalized  Ebers- 
Moll  Model  for  the  large-signal,  dynamic  behavior  of  the  integrated 
bipolar  transistor.  This  examination  discloses  two  related  unsolved 
problem  areas:  the  accounting  for  high-level  behavior  while  retaining 
an  accuracy-simplicity  compromise,  and  the  determination  of  model 
parameter  values  for  device  operation  over  a  wide  range  of  currents. 

Three  approaches  for  attacking  these  problems  are  suggested  and  compared. 

Attention  then  turned  to  the  need  for  modifying  conventional  transistor 
models  lo  adapt  them  for  use  in  existing  network-analysis  programs,  which 
led  to: 

A-3.2  INTEGRATED-CIRCUIT  TRANSISTOR  AND  DIODE  MODELS  FOR  NETWORK-ANALYSTS 
—  —  _____  _____ - - - - - 

PROGRAMS :  A  set  of  new  models  of  varying  complexity  and  accuracy  is 

derived  for  the  large-signal  transient  behavior  of  bipolar  transistor  and 
diode  structures  found  in  integrated  circuits.  The  derivation  starts  with 
the  defining  equations  of  the  so-called  nonlinear  model  for  the  four-layer 
integrated-circuit  transistor  and  proceeds  by  reformulating  several  terms 
in  these  equations  to  remove,  in  effect,  current  sources  in  the  symbolic 
model  dependent  on  the  complex  frequency  variable  4.  The  removal  of  this 
dependence  yields  models  containing  resistors,  capacitors,  and  nonfrequency- 
dependent  sources.  With  the  availability  of  these  models,  network-analysis 
programs  unable  to  accommodate  frequency-dependent  sources  become  directly 
applicable  to  the  computation  of  the  large-signal  transient  behavior  of 
bipolar  integrated  circuits.  Considerations  pertaining  to  practical 
application  guide  the  choice  of  approximations  used  in  the  derivation.  The 
limits  of  validity  of  these  approximations  and  of  others  that  set  the  range 
of  validity  of  the  resulting  models  receive  careful  attention. 

Related  to  this  general  problem  of  tailoring  transistor  circuit  models 
to  the  needs  of  computer-aided  circuit  analysis  programs,  we  reported  a 
simple,  pragmatic  method  of  describing  an  aspect  of  device  physics  neglected 
in  previous  circuit  models? 
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A"3-3  INCOFPORATION  of  the  early  effect  in  the  EBERS-MOLL  MOm?T..4»5 
The  addition  of  a  single  parameter,  termed  the  Early  voltage,  enables 
incorporation  of  the  Early  effect  in  the  Ebers-Moll  model  for  the  large- 
signal  static  behavior  of  transistors.  The  consequent  voltage  dependences 
of  collector  current,  common-emitter  current  gain,  and  small-signal  output 
conductance  are  set  forth. 

To  this  point,  attention  had  focussed  mainly  on  the  bipolar  transistor. 
Advances  in  fabrication  technology  had  made  it  plain,  however,  that  the 
insulated-gate  field-effect  transistor  (IGFET)  would  soon  challenge  the 
bipolar  transistor  as  the  active  device  most  widely  used  in  electronic 
systems.  For  the  IGFET  and  related  devices,  a  clear  need  existed  to 
demonstrate  the  relations  among  the  available  circuit  models.  This 
motivated  the  following  study. 

A-3.4  UNIFIED  MODELING  OF  FIELD-EFFECT  DEVICES : 8 

A  unifying  approach  is  described  to  the  derivation  of  models  for 
field-effect  devices.  The  most  general  model  yielded  by  this  approach 
contains,  as  special  cases,  models  for  the  small-  and  large-signal  behavior 
of  the  two  major  types  of  field-effect  transistor  and  of  other  related 
devices  such  as  the  semiconductor  current  limiter  and  the  pinch  resistor. 
From  this  most  general  model  emerge,  by  application  of  particular 
approximations  and  constraints,  several  models  proposed  earlier  for  use 
in  computer-aided  circuit  analysis.  The  unifying  approach  thus  demonstrates 
the  relationships borne  among  various  of  the  existing  models  for  field-effect 
devices.  In  addition,  it  provides  a  wide  diversity  of  new  models  of 
^ erent  complexities  and  degrees  of  accuracy. 

Full  utilization  of  the  computer  for  the  design  of  semiconductor 
circuits  requires  that  the  parameters  of  circuit  models  for  devices  be 
rapidly,  accurately,  and  economically  determined.  During  the  last  year 
of  this  contract  and  continuing  to  the  present,  a  portion  of  our  research 
effort  has  yielded  considerations  pertinent  to  the  design  of  a  facility  for 
automating  the  determination  of  the  parameters  of  bipolar  and  field-effect 
transistor  circuit  models.  The  results  achieved  to  date  ^’8  bear  however, 
not  only  on  this  problem  but  also  on  the  broader  problem  of  the  systematic 
design  ol  electronic  circuits. 
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The  automation  facility  desired  must  supply  parameters  enabling 
computer  simulation  of  circuit  behavior  by  any  of  the  major  TREE  analysis 
codes.  Of  central  concern  is  the  effect  on  device  and  circuit  performance 
of  neutron  dosage,  gamma  rate  and  pulse  power.  The  simulations  of 
interest  must  include  the  capability  of  predicting  response  to  large- 
signal  and  time-varying  excitation. 

In  planning  a  system  capable  of  achieving  these  objectives,  one  must 
first  define  the  model  whose  parameters  are  to  be  measured.  Upon  this 
requirement,  in  a  surprinsingly  subtle  way,  hinges  to  a  large  degree  the 
success  or  failure  of  the  desired  measurement  system.  The  approach  taken 
in  our  effort  can  be  understood  only  in  appropriate  perspective,  which 
we  now  briefly  try  to  supply. 

Two  decades  of  effort  have  made  available  many  different  transistor 
models,  which  provide  different  compromises  between  accuracy  and  mathematical 
complexity  and  which  contact  to  different  degrees  with  the  processing  used 
in  device  fabrication.  The  existence  of  this  vast  reservoir  of  available 
models  helps  yet  also  hinders  the  designer  of  a  system  for  parameter 
measurement. 

Two  approaches  are  apparent.  Either  one  can  select  from  the  reservoir 
the  single  existing  model  that  at  present  appears  to  be  best  suited  to 
the  needs  of  the  radiation-effects  community,  and  design  a  facility  capable 
of  automatically  yielding  its  parameters.  Or,  alternatively,  one  can  try 
to  design  a  facility  able  to  measure  the  parameters  of  any  of  the  existing 
models  and,  if  possible,  of  any  model  likely  to  be  developed  in  the  future. 

History  argues  in  favor  of  this  second  approach.  The  past  two  decades 
have  seen  continual  development  of  new  and  increasingly  more  adequate 
models,  and  no  reason  exists  to  suppose  that  this  development  will  soon 
cease.  Indeed,  the  smaller  and  more  sharply  defined  geometries  promised 
by  emerging  technologies  such  as  electron  beam  techniques  and  ion  implanta¬ 
tion  create  a  pressing  need  for  new  and  better  models.  Thus  a  parameter 
measuring  system  tailored  to  any  single  existing  transistor  model  is  likely 
doomed  to  early  obsolescence.  It  is  by  far  the  wiser  approach  to  design 
the  system  to  be  compatible  with  all  of  the  existing  models  and,  insofar 
as  possible,  with  all  models  of  the  future. 
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In  the  way  of  this  approach  stands  the  chaos  of  the  existing  large- 
signal  transistor  models.  Discussions  concerning  the  many  alternative 
models  available  are  widely  scattered  in  the  literature;  and,  for  the  most 
part,  the  treatments  leave  unclear  the  relationships  among  the  different 
models.  To  the  circuit  designer  and  analyst,  an  understanding  of  these 
relationships  and  of  the  range  of  model  complexity  and  accuracy  available 
to  his  use  is  highly  desirable.  To  the  designer  of  an  effective  and 
sufficiently  general  facility  for  automatic  parameter  measurement,  such 
an  understanding  is  essential. 

The  main  thrust  of  our  effort,  therefore,  has  been  to  simplify, 
reduce,  and  unify  the  results  of  previous  studies  to  a  form  that  clarifies 
the  relationships  among  the  models,  the  hierarchy  of  model  complexities 
and  accuracies,  and  the  limits  of  model  validity.  We  have  demonstrated 
that  all  existing  models  are  special  cases  of  a  general  model,  and  have 
given  reasons  for  the  belief  that  future  models  will  likewise  be  special 
cases.  The  general  model  is  described  by  a  matrix.  The  problem  of 
designing  a  sufficiently  gener  il  system  for  automated  parameter  measurement 
thus  hinges  on  the  problem  of  determining  the  parameters  that  define  the 
elements  of  this  matrix. 

The  origin  of  the  chaos  that  persists  in  transistor  modeling,  despite 
decades  of  intense  and  well-financed  effort,  lies  in  the  difficulty  of 
the  boundary-value  problem  describing  transistor  behavior.  That  is  where 
our  study  has  begun.  We  examine  the  approximations  that  have  conventionally 
been  made  in  dealing  with  this  boundary-value  problem  and  are  thereby  led 
to  the  general  model  and  to  the  physical  basis  and  approximations  under¬ 
girding  it.  We  then  show  that  the  general  model  contains  as  special  cases 
all  existing  device  circuit  models. 

The  availability  of  the  general  model  unifies  the  phenomenological 
inclusion  of  radiation  effects;  it  proves  to  be  as  easy  to  include  them 
in  any  one  model  as  in  any  other. 
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A-4  Selection  of  Circuit  Models 

For  many  years,  research  in  semiconductor  device  characterization 
focussed  nearly  exclusively  on  the  first  three  of  the  modeling  problems 
indicated  in  Fig.  II-l.  In  1971,  under  the  sponsorship  of  this  contract, 
we  pioneered  research  concerning  the  fourth  problem  which  deals  with 
selection  of  models. 

Success  in  the  simulation  and  automated  design  of  integrated  circuits 
depends  critically  on  the  models  chosen  to  represent  the  behavior  of 
active  devices.  The  basic  consideration  involved  in  the  choice  is  to 
obtain  required  accuracy  from  a  model  of  minimum  complexity.  Basing 
calculations  on  a  model  of  insufficient  accuracy  will  normally  lead  to 
meaningless  results.  On  the  other  hand,  use  of  a  model  more  complex  than 
needed  threatens  the  capability  of  computer  memory,  increases  CPU  time, 
increases  the  chance  of  numerical  errors,  and  limits  the  size  of  circuits 
that  the  simulation  program  can  accommodate.  Moreover,  the  more  complex 
the  model  the  more  time  and  effort  must  be  spent  determining  the  model 
parameters.  A  design  using  optimization  techniques  compounds  the  need 
for  model  simplicity  because  such  a  design  requires  many  circuit  analyses. 

To  counter  these  problems  we  described  a  method  *  that  selects,  for 
each  transistor  in  a  circuit,  the  model  of  least  complexity  that  will 
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give  acceptable  accuracy.  At  the  heart  of  the  method  is  the  capability 
to  assess  model  adequacy.  This  capability  derives  from  a  self-consistency 
test  in  which  the  values  of  currents  and  voltages  computed  in  a  simulation 
of  the  circuit  behavior  are  compared  with  onset  parameters.  By  this 
comparison  one  determines  whether  these  computer  values  are  consistent 
with  the  approximations  underlying  the  device  models  used  in  the  simulation. 
The  onset  parameters  for  a  model  are  the  terminal  currents  and  voltages 
above  or  below  which  the  model  fails  to  give  a  satisfactory  representation 
of  device  behavior.  We  set  forth  the  onset  parameters  for  the  Ebers-Moll 
model  and  discuss  their  determination  by  terminal  measurement  and  by 
calculation  based  on  the  transistor  makeup.  To  illustrate  the  method  of 
model  assessment,  we  simulated  the  performance  of  an  operational  amplifier; 
the  effectiveness  of  the  method  in  selecting  model  complexity  was 
demonstrated. 

Two  invited  research  papers  followed  concerning  this  same  subject. 

The  first  revised  the  specification  of  one  of  the  onset  parameters 
described  earlier  and  assessed  the  applicability  of  the  method  to  the 
simulation  of  circuit  behavior  in  a  radiation  environment.  The  second 
paper  made  additional  extensions  of  the  method:  to  accommodate  the 
field-effect  transistor,  rather  than  solely  the  bipolar  transistor;  and 
to  accommodate  the  selection  of  models  representing  time-varying  response, 
rather  than  solely  static  behavior. 

The  core  idea  of  our  method  for  model  selection  has  been  incorporated 
into  one  of  the  most  comprehensive  computer-aided  programs  for  circuit 
analysis  presently  available."5 

REFERENCES  Vindicates  support  by  this  contract) 

*1.  F.  A.  Lindholm  and  S.  W.  Director,  "Assessing  Model  Adequacy  in 

Integration  Circuits  Simulation,"  1971  IEEE  International  Solid  State 
Circuits  Conference:  Technical  Digest.  17  (1971).  - 

*2.  F.  A.  Lindholm,  S.  W.  Director  and  D.  L.  Bowler,  "Assessing  Model 

Adequacy  in  Selecting  Model  Complexity  in  Integrated-Circuit  Simulation," 
IEEE  Journ.  Solid  State  Circuits,  SC-6,  213-222  (1971)  [invited  research 
paper] . 

*3.  F.  A.  Lindholm,  "Device  Characterization  for  Computer  Analysis  of  Large 
Circuits,"  1971  IEEE  Annual  Conference  on  Space  and  Radiation  Effects: 
Digest,  106  (1971)  [invited  research  paper].  ~ 

F.  A.  Lindholm,  "Devices  Characterization  for  Computer  Analysis  of  Large 
Circuits,"  IEEE  Transactions  on  Nuclear  Science.  18,  206-210  (1971) 

(invited  research  paper) .  ~ 


16 


*4.  F.  A.  Lindholm,  "Generation  and  Selection  of  Dynamic  Large-Signal 
Models  for  Bipolar  and  Field-Effect  Transistors,"  1972  International 
IEEE  Symposium  on  Circuit  Theory.  £,  385  (1972)  [invited  paper]. 

5.  H.  J.  DeMan  and  R.  Mertens,  "SITCAP — A  Simulator  of  Bipolar 
Transistors  for  Computer-Aided  Circuit  Analysis  Programs,"  1973 
IEEE  Int.  Solid  State  Circuits  Conference:  Digest.  16,  104  (1973) . 


17 


8.  ELECTRONIC  CERAMICS 


8-1  Introduction 

The  general  objective  of  the  electronic  ceramics  program  sponsored 
in  this  research  contract  was  to  establish  the  relative  importance  of 
microstructure  on  the  electronic  behavior  of  a  potential  new  generation 
of  electronic  materials  and  to  determine  the  influence  of  radiation 
environments  on  the  properties  and  microstructure  of  the  materials.  Two 
major  classes  of  new  electronic  materials  were  selected  for  investigation; 
semiconducting  glasses  and  glass-ceramic  insulators.  Both  types  of 
materials  had  received  little  systematic  investigation  of  microstructural- 
properties-irradiation  relationships  at  the  outset  of  our  research  program. 
In  addition  each  type  of  material  appeared  to  have  considerable  promise 
for  highly  radiation  resistant  components.  Glass-ceramics  were  potentially 
radiation  hard  insulators  or  encapsulants.  Semiconducting  glasses  were 
promising  radiation  insensitive  resistor  elements  and  showed  some  potential 
for  active  device  replacement. 

The  aim  of  our  research  was  to  establish  critical  processing  conditions 
and  to  understand  the  development  of  a  controlled  microstructure  in  these 
materials.  Following  identification  of  nucleation  processes,  we  established 
the  kinetics  of  crystallization  and  other  aspects  of  microstructural 
evolution  with  thermal  processing.  The  relationships  between  thermal 
processing  and  microstructure  to  the  electrical  conductivity,  dielectric 
properties,  mechanical  strength  and  thermal  properties  was  determined. 

And  finally,  the  effect  of  fast  neutron  irradiation,  and  in  some  cases 
gamma  irradiation,  on  the  metrical  properties  of  the  materials  was  measured 
and  related  to  changes  in  the  microstructure  of  the  material. 

A  total  of  32  research  publications  and  four  edited  books  resulted 
from  this  research.  In  addition  approximately  sixty  technical  presentations 
were  made  at  various  international,  national,  and  research  laboratory 
meetings  throughout  this  country  and  abroad.  Seven  advanced  degrees  were 
awarded  for  research  of  students  involved  in  this  program.  Two  summer 
consultantships  for  the  U.  S.  Army  amorphous  semiconductor  program  at 
Picatinny  Arsenal,  Physical  Sciences  Branch  were  held  by  Professor  L.  L. 
Bench.  This  research  project  made  it  possible  for  him  to  introduce  through 
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Che  summer  program  a  materials  processing  and  microstructural  approach 
to  the  large  amorphous  semiconductor  program  in  progress  at  Picatinny 
Arsenal. 

A  summary  of  the  various  aspects  of  the  electronic  ceramics  research 
sponsored  by  this  project  follows  with  appropriate  citations  of  papers 
published  from  the  work  included. 

B-2  SEMICONDUCTING  GLASSES 
B-2.1  Microstructural  Effects 

The  relative  importance  of  thermal  treatment  and  microstructure  on 
the  properties  of  amorphous  semiconductors  was  first  reported  from  research 
sponsored  in  this  project.  (1-5)  Glasses  in  the  V^-P^-I^O  systems, 
where  R20  is  any  alkali  metal  oxide,  were  studied.  These  papers  establish 
that  a  cast  semiconducting  glass  will  possess  electrical  properties  which 
are  a  function  of  the  casting  temperature,  quench  rate,  and  annealing  time. 
The  thermal  variables  alter  a  microphase  separation  in  the  glass  which 
influences  the  electronic  structure.  Variations  in  electronic  structure 
give  rise  to  alterations  in  both  the  DC  electrical  conductivity  and  the  AC 
dielectric  spectra  of  the  materials.  Subsequent  thermal  treatments  after 
casting  result  in  both  a  coarsening  of  the  phase  separation  and  the  onset 
of  crystallization  in  the  semiconducting  glasses.  These  changes  produce 
very  large  alterations  in  the  dielectric  properties  of  the  semiconducting 
glasses  as  a  result  of  the  initiation  of  Maxwell-Wagners-Sillars  hetero¬ 
geneous  dielectric  losses  due  to  the  precipitation  of  high  conductivity 
crystallites  in  the  low  conductivity  glassy  matrix. 

B-2. 2  Irradiation  Effects 

The  first  studies  of  irradiation  effects  involving  bulk  amorphous 
semiconductors  were  performed  in  this  project. These  studies  were 
also  of  glasses  in  the  V205_P205_R2°  system,  where  R20  is  any  alkali  metal 
oxide.  The  research  has  shown  that  the  exposure  of  the  semiconducting 
glasses  to  y-irradiation  increases  the  electrical  conductivity  by  10%  when 

Q 

fluences  in  the  range  of  2  x  10  rads  are  achieved.  This  effect  is 
apparently  due  to  ionization  of  the  multivalence  transition  metal  ions 
in  the  glass.  Fast  neutron  exposures  in  the  range  of  5  x  1016  to  5  x  1017  nvt 
show  that  bulk  semiconducting  glasses  in  the  as-cast  state  are  relatively 
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little  affected  by  fast  neutron  Irradiation.  However,  semiconducting 
glasses  that  have  been  thermally  treated  to  enhance,  or  modify,  the 
electrical  conductivity  are  strongly  influenced  by  the  fast  neutron 
irradiation.  Crystallites  precipitated  in  the  glass  are  disordered 
by  the  fast  neutron  irradiation  exposure  and  the  dielectric  properties 
associated  with  the  isolated  high  conductivity  phase  are  destroyed. 

The  general  effect  of  the  fast  neutron  irradiation  is  to  make  thermally 
treated  semiconducting  glasses  appear  more  like  bulk  semiconducting 
glasses.  Application  of  these  results  to  predicting  the  radiation 
resistance  of  amorphous  semiconductor  switching  devices  led  to  the 
conclusion  of  a  1  x  10^  nvt  damage  threshold. 

B-2,3  Theory 

Analysis  of  microstructural  effects  on  the  transport  properties  of 
amorphous  semiconductors  was  also  examined  from  a  theoretical  standpoint 
in  the  project.  *  ^  The  heterogeneous  semiconductor  theory 

evolved  from  this  study  has  shown  that  many  of  the  AC  characteristics  of 
semiconducting  glasses  are  due  to  the  presence  of  phase  discontinuities 
in  the  materials. 

A  cluster  calculation  for  amorphous  solids  was  also  developed . 

The  quantum  mechanical  calculation  of  the  electronic  bond  structure  of 
amorphous  solids  was  done  with  a  CNDO  semiempirical  model.  Clusters  as 
large  as  50  plus  atoms  were  evaluated.  This  investigation  has  shown  that 
there  is  only  a  very  small  tailing  of  band  states  in  the  forbidden  gap 
even  with  large  degrees  of  disorder  in  the  solid.  Consequently,  there 
is  little  theoretical  basis  for  the  Cohen-Fritsche-Ovskinsky  band  model 
proposed  for  amorphous  semiconductors,  which  involves  an  overlap  of  band 
states  in  the  middle  of  the  forbidden  energy  band.  The  model  proposed 
by  Mott  is  much  closer  to  being  consistent  with  the  calculation  results 
obtained. 

B-3  GLASS-CERAMICS 

The  glass-ceramics  research  in  this  project  was  largely  concerned 
with  the  L^O-SiC^  glass  system.  This  binary  glass  system  is  basic  to 
many  of  the  commercial  L^O-A^O^-SiC^  glass-ceramic  systems  marketed  by 
the  glass  industry.  Consequently,  understanding  the  mechanisms  of 
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nucleation  and  crystallization,  microstructural  factors,  and  irradiation 
damage  in  the  binary  Li20-Si02  system  has  led  to  a  greater  understanding 
of  the  behavior  of  the  commercial  glass-ceramics.  Of  even  more  importance 
is  that  the  binary  glass  system  is  simple  enough  that  it  is  possible  to 
understand  the  fundamental  mechanisms  involved  in  the  development  and 
control  of  microstructure  and  the  interpreation  of  the  irradiation  damage 
effects. 

B~3.1  Nucleation  and  Crystallization  Kinetics 

Identification  of  a  metastable  nucleation  process  in  the  Li20-Si0 

glasses  was  reported  in  a  series  of  papers . (15_20)  A  metastable  lithium 

metasilicate  phase  appears  in  the  very  early  stages  of  heat  treatment  of 

°  ° 

the  glasses.  The  phase  appears  as  200A  by  50A  needle-like  precipitates 
in  the  glass.  The  metastable  phase  occurs  in  very  small  quantitites,  in 
a  range  of  .01  volume  percent  or  less.  Identification  of  this  phase 
required  the  use  of  a  large  number  of  investigative  tools  including  hot 
stage  transmission  electron  microscopy,  x-ray  small  angle  scattering, 
quantitative  dielectric  relaxation,  and  x-ray  diffraction  analysis  using 
a  Guinier-DeWolff  camera.  Applying  the  combination  of  these  specific 
tools  to  this  problem  together  in  a  concerted  effort  to  identify  the 
fundamental  nucleation  species  was  the  first  such  effort  in  the  glass 
literature.  The  results  obtained  have  helped  to  clarify  many  of  the  more 
qualitative  and  less  detailed  investigations  reported  by  a  number  of  other 
authors.  The  results  have  also  been  applicable  to  understanding  the 
mechanicsms  of  nucleation  in  a  number  of  multicomponent  glasses  reported. 

By  identifying  the  nucleation  mechanism  it  was  possible  to  apply 
this  understanding  to  controlling  the  microstructure  of  crystallized 
glass-ceramics  in  the  binary  Li20-Si02  system  without  resorting  to  the 
addition  of  heterogeneous  nucleating  agents.  The  kinetics  of  crystallization 
were  measured  using  both  quantitative  x-ray  diffraction  and  quantitative 
optical  microscopy  methods.  Again  the  confirmation  of  the  kinetics  using 
a  variety  of  tools  was  the  first  such  effort  in  the  glass-ceramics  literature. 
The  results  obtained  explain  a  number  of  the  inconsistencies  of  the 
activation  energies  reported  for  crystallization  in  this  system  by  other 
invest  Lgators . 
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B-3 . 2  Microstructural  Effects  on  Propertiss 

Having  determined  the  nucleation  and  crystallization  kinetics  of 

the  Li20-Si02  glass  made  it  possible  to  produce  a  series  of  Li20-Si0 

glass-ceramics  with  controlled  volume  franctions  and  particle  sizes  of 

the  Li2Si20,-  crystalline  phase.  The  mechanical  and  dielectric  properties 

of  the  controlled  microstructure  glass-ceramics  were  measured  and  the 

data  analyzed  using  quantitative  theories  of  mixing  effects  of  second 
(21-22) 

phases.  It  was  shown  that  a  major  factor  in  controlling  the 

strength  of  class-ceramics  was  the  development  of  a  dense  glassy  matrix 
prior  to  the  crystallization  of  the  bulk  of  the  glass.  If  the  matrix 
densif ication  does  not  occur,  large  interfacial  stresses  were  developed 
during  the  crystallization  process  and  cracks  resulting  from  the  inter¬ 
facial  stresses  markedly  reduced  the  strength  of  the  glass-ceramics. 
Quantitative  microstructural  effects  on  the  dielectric  properties  showed 
that  the  properties  were  markedly  changed  when  the  crystallized  phase 
became  connected. 

In  both  of  these  studies  the  systematic  development  of  a  controlled 
microstructure,  the  quantitative  analysis  of  that  microstructure,  and 
the  measurement  of  its  effect  on  properties  were  the  first  such  studies 
reported  in  the  glass-ceramic  literature. 

B-3. 3  Irradiation  Effects 

The  influence  of  fast  neutron  irradiation  on  nucleation  processes 
in  the  glass-ceramics  were  investigated. (7,8,23)  These  studies  showed 
that  the  metastable  precipitates  became  disordered  when  exposed  to  a 
neutron  fluence  of  1  x  1017  nvt  or  greater.  As  a  result,  the  relative 
effectiveness  of  the  nucleation  process  was  decreased.  However,  we 
also  discovered  during  the  course  of  this  investigation  that  exposure 
of  the  glass  to  the  neutron  irradiation  resulted  in  the  formation  of 
structural  heterogeneties  which  increased  the  efficiency  of  subsequent 
thermal  nucleation  treatments. 

A  major  finding  in  the  irradiation  investigation  is  that  all  crystal¬ 
lites  of  dielectric  solids  become  progressively  disordered  as  a  result 
of  the  neutron  irradiation.  The  progressive  disordering  decreases  the 
electrical  conductivity  and  eliminates  progressively  the  AC  dielectric 


22 


losses  associated  with  the  presence  of  an  isolated  second  phase.  To  our 
knowledge  these  studies  of  neutron  irradiation  effects  in  glass-ceramics 
are  still  the  only  investigations  in  this  field. 

B-4  CHARACTERIZATION  METHODS 

In  order  to  pursue  the  research  topics  discussed  above  it  was  necessary 

to  develop  or  extend  several  analytical  tools  to  serve  as  characterization 

methods  for  glasses  and  ceramic  materials.  A  general  review  of  a  philosophy 

(24) 

of  characterization  was  published  and  a  book  was  edited  on  the  subject 

(2 

of  quantitative  characterization  of  ceramics  during  a  portion  of  this  project. 
B-4.1  X-ray  Methods 

Limitations  of  sample  preparation  methods  on  the  accuracy  of  x-ray 
spectrochemical  analysis  were  examined  in  order  to  eliminate  this  trouble¬ 
some  source  of  efforts  in  the  analysis  of  our  experimental  samples  in  the 
(  26} 

research  program.  A  major  effort  was  expended  in  the  application  of 

a  new  type  of  x-ray  powder  camera,  called  the  Guinier-DeWolf f  x-ray 

(27) 

diffraction  camera.  This  camera  was  used  tc  determine  the  volume 

fraction  of  very  small  quantities  of  second  phases  in  the  semiconducting 
glasses  and  the  glass-ceramics  being  investigated.  This  study  resulted 
in  the  quantitative  use  of  the  Guinier  cameras  which  made  possible  the 
identification  of  the  nucleation  mechanisms  and  structural  changes  in  the 
electronic  materials  being  studied. 

B-4. 2  Microscopy 

Sample  preparation  techniques  were  developed  for  obtaining  measure¬ 
ments  of  nucleation  and  crystallization  processes  of  glass-ceramics  in 
the  electron  microscope.  These  measurements  were  done  using  the  hot  stage 
electron  microscope  which  made  it  possible  to  see  the  sequence  of  the 
reaction  mechanisms  operating.  Possible  errors  in  obtaining  quantitative 

analyses  of  microstructures  in  the  semiconducting  glasses  and  glass-ceramics 

(28) 

were  delineated  by  a  study  of  thickness  error  effects  in  optical  microscopy. 

This  study  made  it  possible  to  apply  mathematical  corrections  to  the  quanti¬ 
tative  microscopy  measurements  when  needed. 
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B-4.3  Surface  Properties 

Many  electronic  materials  are  processed  using  powder  ceramic  techniques. 
Understanding  the  surface  behavior  of  these  powders  is  at  a  primitive  stage 
of  evolution  at  this  time.  Some  of  the  characterization  methods  success¬ 
fully  used  in  other  parts  of  this  research  project  appear  to  be  suitable 
for  understanding  surface  powder  characterization.  The  potential 
application  of  the  dielectric  and  various  x-ray  techniques  and  electron 
microscopy  techniques  were  explored  and  reported  in  this  project. (29-3°) 

B-4.4  Dielectrics 

Dielectric  relaxation  is  especially  sensitive  to  the  presence  of 

isolated  heterogeneous  phases.  Consequently,  the  AC  spectra  of  many  of 

the  materials  being  investigated  is  strongly  influenced  by  the  presence 

of  nucleated  and  crystallized  second  phases.  This  development  of 

quantitative  dielectric  relaxation  as  a  microstructural  analysis  tool  is 

one  of  the  major  outcomes  of  the  characterization  part  of  this  research 
(17  29-32) 

program.  ’  This  tool  is  now  used  by  a  number  of  investigators 

throughout  the  country  to  understand  the  early  stages  of  second  phase 
development  in  exploratory  materials. 

B-5  SILICON  DEFECT  AND  IRRADIATION  DAMAGE  CALCULATIONS 

The  computer  based  cluster  calculations  developed  for  understanding 
the  electronic  structure  and  properties  of  an  amorphous  semiconductor 
are  also  suitable  for  analyzing  the  many  complicated  defect  structures 
possible  in  crystalline  silicon  semiconductors.  A  series  of  progressively 
more  complex  defects  were  examined  in  this  study,  including  impurity  pairs, 
impurity-vacancy  complexes,  and  multi-impurity  complexes. ^33-35^  The 
energy  levels  associated  with  the  impurities  obtained  from  the  cluster 
calculations  correlated  extremely  well  with  the  literature  when  literature 
values  were  available  for  comparison.  Many  defects  calculated  have  been 
theorized  in  the  literature  without  any  calculational  basis  for  the  energy 
levels  reported.  Application  of  this  method  to  understanding  various 
processing  problems  and  irradiation  effects  in  low  doped  and  degenerate 
silicon  semiconductor  devices  are  possible  from  this  research. 
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B-6  CONFERENCES  AND  BOOKS 


The  first  conference  on  the  Physics  of  Electronic  Ceramics  was 
organized  and  hosted  at  the  University  of  Florida  during  the  course  of 
this  project.  The  Conference  brought  together  over  150  industrial, 
university,  and  government  scientists  to  review  the  status  of  theory  of 
electronic  ceramics.  Two  books  resulted  from  the  Conference;  "Physics 
of  Ceramics,"  vol.  A  and  "Physics  of  Ceramics,"  vol.  B.  Both  were  edited 
by  L.  L.  Hench  and  D.  B.  Dove. 

A  timely  book,  "Characterization  of  Ceramics,"  was  edited  by  L.  L. 
Hench  and  R.  W.  Gould  during  the  period  of  the  contract.  Many  chapters 
in  this  book  provide  guidelines  important  to  the  characterization  of 
all  types  of  electronic  materials. 

Also,  during  the  project  an  international  symposium  on  nucleation 
and  crystallization  in  glasses  was  organized  by  L.  L.  Hench  for  the 
American  Ceramic  Society.  This  led  to  Special  Publication  #5  by  the  ACS 
titled  "Advances  in  Nucleation  and  Crystallization  in  Glasi  ,‘S,"  L.  L. 

Hench  and  S.  W.  Freiman,  editors.  This  volume  provides  a  thorough  summary 
of  the  major  developments  in  the  field  of  nucleation  and  crystallization 
of  glasses  in  the  1960 's  and  early  70’ s. 

These  four  books  constitute  a  significant  library  of  up-to-date 
information  on  physical  ceramics  and  glasses.  This  research  project 
provided  the  impetus,  motivation,  and  much  of  the  intensive  research 
climate  within  the  University  of  Florida  electronic  materials  group 
required  to  make  the  conferences  a  success  and  the  publications  a  reality. 
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C.  QUANTUM  AND  ELECTRONIC  PROPERTIES  OF  SEMICONDUCTORS  AND  RELATED  DEVICES 


C-l  Introduction 

During  the  period  of  contract  support,  a  substantial  effort  has  focussed 
on  the  experimental  and  tneoretical  studies  in  semiconductor  materials 
and  devices  that  bear  particular  interest  to  DOD  applications.  These 
studies  include: 

(a)  calculation  of  the  diffusivity-mobility  ratio  and  the  density  of 
states  function  for  degenerate  and  nondegenerate  semi¬ 
conductors  ; 

(b)  experimental  and  theoretical  investigation  of  the  recombina¬ 
tion  and  trapping  mechanisms  of  the  injected  carriers  in 
the  defect  centers  of  such  semiconductors  as  silicon,  gold- 
doped  silicon,  indium-arsenide,  chromium-doped  gallium 
arsenide,  nickel-doped  germanium  and  indium-doped  cadmium- 
sulfide.  Photomagnetoelectric  and  photoconductive  methods 
have  been  used  extensively  to  yield  values  of  important 
parameters  such  as  carrier  lifetimes,  capture  rate  and 
mobilities  over  a  wide  temperature  range;  and 

(c)  experimental  and  theoretical  studies  of  photodetectors. 

This  study  resulted  in  the  development  of  a  new  grating 
type  silicon  Schottky  barrier  photodiode  which  shows 
superior  performance  characteristics  over  0.4~1.06ym 
spectral  range.  In  addition,  it  gave  a  new  understanding 
of  the  quantum  yield  of  metal-semiconductor  photodiodes 
and  p-i-n  photodiodes. 

As  a  result  of  this  research  effort,  twenty-four  papers  have  been 
published  in  national  and  international  technical  journals  and  six 
papers  were  presented  at  the  national  or  international  conferences. 

A  summary  follows  of  the  most  important  results. 
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C-2  STUDIES  OF  DIFFUSIVITY-MOBILITY  RATIO  AND  DENSITY  OF  STATES  FUNCTION 
FOR  DEGENERATE  SEMICONDUCTORS 

C-2.1  Generalized  Einstein  Relation  ^ 

An  expression  has  been  derived,  formulated  in  terms  of  tabulated 
functions,  for  a  generalized  Einstein  relation  valid  for  degenerate  as 
well  as  nondegenerate  semiconductors.  To  facilitate  the  use  of  this 
expression  in  device  design,  a  graph  has  been  constructed  displaying 
the  dif fusivity-mobility  ratio  against  carrier  concentration;  the  graph 
holds  for  any  temperature  and  for  any  materials  whose  density  of  quantum 
states  varies,  to  a  reasonable  approximation,  as  the  square  root  of  energy. 
Separate  curves  have  been  plotted  for  the  special  cases  of  silicon  and 
germanium.  For  any  level  of  concentration,  these  curves  afford  a  comparison 
with  the  constant  value  of  the  dif fusivity-mobility  ratio  predicted  by  the 
conventional  Einstein  relation,  and  display  the  error  deriving  from  its 
misuse. 

C-2. 2  Density  of  States  Function  for  Degenerate  Semiconductors  ^ 

Theoretical  calculation  of  density  of  quantum  states  as  a  function  of 
impurity  concentration  has  been  made  for  silicon  and  germanium.  The  result 
indicates  that  deviation  of  the  density  of  states  function  from  conventional 
theory  is  negligible  for  electron  energy  greater  than2.kT  above  the  conduction 
band  edge.  However,  the  tailing  effect  of  density  of  states  function  in  the 
forbidden  gap  becomes  eminent  for  ND  >  5xl019cm“3  for  Si  and  N  >  1019Cm“3  for 
Ge. 

C-3  STUDIES  OF  RECOMBINATION  AND  TRAPPING  MECHANISMS  IN  SEMICONDUCTORS  3Y 
USING  PHOTOMAGNETOELECTRIC  (PME)  AND  PHOTOCONDUCTIVE  (PC)  METHODS 

C-3.1  Photomagnetoelectric  and  Photoconductive  Effects  in  Phosphorous-Doned 
Silicon  at  Low  Temperatures  d  ~  — c — 

The  PME  and  the  PC  effects  in  phosphorus-doped  silicon  have  been  investigated 

as  functions  of  illumination  intensity  and  magnetic  field  at  21.8°K  and  4.2°K. 

Concurrent  measurements  of  the  PME,  PC  and  the  Hall  effects  in  this  sample 

yield  the  ambipolar  diffusion  length,  ambipolar  diffusion  constant  and  the 

carrier  lifetimes.  Theories  for  the  PME  and  PC  responses  are  developed  for 

the  case  of  high  injection  (strong  illumination)  and  arbitrary  magneti c  field. 
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The  experimental  results  are  found  in  good  agreement  with  the  theoretical 
predictions.  Study  of  the  PME  effective  lifetime  versus  the  excess  carrier 
densities  has  been  made  for  this  sample.  The  results  indicate  that  the 
Shockley-Read  recombination  model  serves  adequately  for  interpreting  the 
observed  data.  The  Quantity,  t^,  as  determined  from  S-R  models,  is  found 
in  excellent  agreement  with  that  obtained  from  PME  theory.  In  addition, 
it  is  found  that  the  effect  of  trapping  can  be  neglected  over  the  whole 
injection  range  for  this  sample. 

C-3.2  Recombination  and  Trapping  Processes  of  the  Injected  Carriers  in 
Gold-Doped  Silicon^ 

Steady-state  recombination  and  trapping  processes  of  the  injected 
carriers  in  gold-  and  phosphorus-doped  silicon  have  been  studied  for  the 
case  that  the  doping  and  temperature  range  are  such  that  the  equilibrium 
carrier  densities  are  small  compared  with  the  impurity  density.  The 
fluctuation  of  the  charge  in  the  impurity  centers  as  a  result  of  injection 
is  taken  into  account  in  the  charge-neutrality  condition.  The  relationship 
between  the  injected  carriers  (electrons  and  holes)  for  charge  neutrality 
is  established  in  the  form  of  power  laws  (i.e.,  Ap  =  TAn01)  that  hold  in 
different  ranges  of  injection.  Numerical  analyses  of  Ap  versus  An  and  t 
versus  An  are  presented  for  the  case  of  n-type  silicon  overcompensated  and 
undercompensated  by  the  gold  impurities. 

C-3.3  Low-Temperature  Photomagnetoelectric  Properties  of  Gold-Doped  n-Type 
Silicon^  ~~  '  ~  ~ 

The  photomagnetoelectric  (PME)  effect  was  observed  between  21  and  84°K 
in  silicon  doped  with  gold  and  phosphorus.  A  generalized  diffusion  equation 
is  formulated  and  solved  with  the  impurity  centers  included  in  the  charge- 
balance  equation  and  arbitrary  intensity  of  steady  illumination.  When  the 
carrier  densities  are  much  smaller  than  the  gold  density,  the  charged 

impurities  maintain  charge  neutrality,  and  this  results  in  a  range  of 

2 

injection  where  Ap  =  T.An  +  T  An  .  If  the  quadratic  term  dominates, 
the  PME  current  is  proportional  to  the  y  power  of  the  photoconductance. 

This  power  law  is  observed  in  our  measurements,  from  which  we  deduce  values 
for  effective  carrier  lifetimes  between  36  nsec  and  63  ysec. 

C-3.4  Photomagnetoelectric  and  Photoconductive  Effects  in  n-Type  InAs^ 

A  study  of  the  photomagnetoelectric  (PME)  and  photoconductive  (PC) 
effects  in  n-type  InAs  single  crystals  has  been  made  between  4.2  and  300°K 
under  small-signal  conditions.  It  was  found  that  for  T  >  77°K  the  PME  open- 


31 


circuit  voltage  is  proportional  to  the  magnetic  induction,  and  for 
T  >  21°K  the  PME  open-circuit  voltage  saturates  at  large  magnetic 
induction,  which  is  in  accord  with  the  large-Hall-  angle  PME  theory. 
Carrier  lifetimes  were  computed  from  the  PME  and  PC  data  between  4.2 
and  300°K  ,  and  the  hole  mobility  was  estimated  from  the  PME  data  at 
21  and  4.2°K.  The  observed  linear  relationship  between  the  PME  open- 
circuit  voltage  and  photoconductance  at  low  light  levels  manifests 
that  the  effect  of  trapping  is  essentially  nil  over  the  entire  tem¬ 
perature  range  from  4.2  to  300°K. 


C-3.5  Investigation  of  the  Recombination  and  Trapping  Processes  of 
Photoinjected  Carriers  in  Semi-Insulating  Cr-Doped  GaAs  Usinp 
PME  and  PC  Methods^ 


Investigation  of  the  recombination  and  trapping  processes  of 

g 

photoinjected  carriers  in  semi-insulating  Cr-doped  GaAs  (p  ~  10  ft 
cm)  has  been  made  at  80  and  300°K,  using  the  photomagnetoelectric 
(PME)  and  photoconductive  (PC)  effects  under  large-injection  condi¬ 
tions.  A  generalized  theory  for  the  PME  and  PC  effects  is  developed, 
taking  into  account  the  variation  of  carrier  lifetimes  with  injected- 
carrier  density  (i.e.,  xr  =  <Anp)  and  of  the  trapping  of  holes  in 
the  Cr  levels  (i.e.,  Ap  =  TAn,  T<1),  by  using  a  simple  recombination 
and  trapping  model.  The  result  yields  a  power-law  relationship 
between  the  PME  short-circuit  current  and  the  photoconductance,  in 
a  form  that  Ipj^g  varies  as  AG^'^+B^  valid  for  different  injection 
ranges.  Two  well-defined  ranges  of  injection  were  observed  from  the 
present  results:  In  region  I  (high  injection),  I  is  directly 

FMcj 

proportional  to  AG  (i.e.,  B=0) ,  and  in  region  II  (intermediate  injection), 

1.2  1 

IDMP.  varies  with  AG  ’  (i.e.,  B  =  -  -?) .  Numerical  values  of  X  ,  X  , 

trie-  3  n  p 

and  Ipj^g  as  functions  of  AG  are  given.  The  results  are  in  good 
agreement  with  the  theoretical  predictions. 

C-3.6 

Although  room-temperature  electron  and  hole  capture  rates  for 
nickel  atoms  in  germanium  have  been  previously  reported,  no  conclu¬ 
sive  result  has  as  yet  been  reached.  This  is  due  to  the  fact  that 
carrier  capture  rates  for  nickel-doped  germanium  reported  previously 
were  determined  mainly  from  the  temperature  dependence  of  lifetimes, 
which  is  further  complicated  and  prone  to  obtain  incorrect  results 
when  the  capture  rates  become  temperature  dependent.  In  this  report 
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'*i  lave  demonstrated  that  by  including  the  excess— carrier  trapping 
effect  in  the  theory  of  photomagnetoelectric  (PME)  and  photoconduc¬ 
tance  (PC)  effects  and  by  elaborately  controlling  the  compensation 
ratio  of  nickel  and  shallow- impurity  densities  in  germanium,  we  are 
able  to  determine  accurately  three  of  the  four  capture  rates  for 
nickel  atoms  in  germanium  from  room-temperature  PME  and  PC  measure¬ 
ments.  In  addition,  with  the  present  results,  we  are  able  to  inter¬ 
pret  the  previous  results  consistently  at  room  temperature.  The 
present  results  are  Cnl  =  1.4  X  10_9,  =  8  X  10"9,  and  C  2=4X10~8 

in  units  of  cm8  sec  3  at  T  =  300°K.  ^ 

C-3.7  Electron  Mobility  and  Shallow  Impurity  Levels  in  In-Doped  and 
Cu-Doped  Cadmium- Sulfide 

Electrical  conductivity  and  Hall  effect  measurements  were  made 

between  300°K  and  4.2°K  on  In-doped  (~6xl017cnf  3)  CdS  single  crystals. 

Normal  band  conduction  was  observed  between  300°K  and  20°K,  while  the 

impurity  hopping  conduction  in  the  second  shallow  donor  state  (E  ) 

D2'  ’ 

dominated  the  conduction  mechanism  for  T  <  20°R.  In  the  normal  band 

conduction  region,  the  best  theoretical  fit  was  obtained  by  using  a 

two-shallow- impurity  level  model.  The  result  yields  E--  =  0.033  eV 

17  o  D1 

for  In-level  with  NQ1  =  6  x  10  cm  J  and  ED2  -  0.007  eV  (which  is 
unknown)  with  NQ2  =  1.3  x  1017cm  3,  assuming  =  1  x  1017cm~3  and 
me  =  °*2  mQ*  The  Hall  mobility  versus  temperature  was  studied  for 
both  Cu-doped  and  In-doped  CdS  between  4.2°K  and  300°K.  It  was  found 
that  the  longitudinal  optical  phonons  dominate  the  scattering  of 
electrons  for  T  >  200°K  and  with  yL  =  110  [exp(300  ±  10) /T  -  l]cm2/V-sec. 
A  brief  discussion  of  the  scattering  mechanism  for  T  <  100°K  is  given 
for  both  In-doped  and  Cu-doped  CdS. 

3.8  Investigation  of  Recombination  and  Trapping  Processes  in 

^•morPh°us  Films  by  Photoconductivity  Method3 ® 

Photoconductivity  experiment  has  been  performed  on  high  resis- 

9 

tivity  (p  =  10  f2-cm) 3As2Se^-2Sb2Se^  amorphous  films  in  the  temperature 
range  from  278°K  to  308°K,  as  a  function  of  light  intensity,  I  .  The 
results^are  that  at  intermediate  light  levels  the  photocurrent  varies 
as  *0  and  at  light  levels  the  photocurrent  is  directly  propor¬ 

tional  to  Iq.  A  simple  recombination  and  trapping  model  is  formulated 
to  interpret  the  observed  photoconductivity  data.  From  the  temperature 
dependence  of  the  photocurrent,  an  effective  trap  level  located  0.3eV 
below  the  band  edge  is  deduced. 
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C-3 . 9  Radiative  Recombination  in  Q^-Doped  N-Type  GaAs  at.  Low  Temperatures  1 1 
Band  to  band  radiative  recombination  in  C^-doped  n-type  GaAs  has 
been  investigated  at  20.8°K  and  4.2°K  using  photomagnetoelectric  (PME) 
and  photoconductive  (PC)  methods.  Capture  rates  determined  from 
the  present  experimental  data  are  in  excellent  agreement  with  that 
computed  from  the  Hall’s  direct  radiative  recombination  formula.^ 

The  results  are  Br  =  1.15  X  10~8cm3/sec  at  20.8°K  and  1.23  X  10-7cm3/sec 
at  4. 2°K. 

C-4  THEORETICAL  AND  EXPERIMENTAL  STUDIES  OF  METAL- SEMICONDUCTOR 
SCHOTTKY- BARRIER  AND  PME  PHOTODETECTORS 

12 

C-4.1  A  New  Grating  Type  Gold-N  Type  Silicon  Schottky  Barrier  Photodiode 
A  new  grating  type  Au-nSi  Schottky-barrier  photodiode  has  been 
fabricated  and  examined.  The  device  is  constructed  with  a  new  struc¬ 
ture  of  contact  mask.  The  mask  forms  a  grating  of  gold  film  over  the 
contacting  silicon  substrate.  The  resulting  device  shows  features 
such  as  broad  bandwidth,  high  responsivity ,  high  quantum  yield  and 
low  noise  characteristics.  Report  is  given  here  for  the  I-V  charac¬ 
teristics,  the  spectral  dependence  of  quantum  yield,  responsivity, 
noise  performance  and  the  technology  used  to  fabricate  such  device. 

Without  optimization,  the  present  device  shows  a  bandwidth  of  ~lGHz, 
a  responsivity  of  0.44  A/W  at  0.63ym  and  0.63  A/W  at  X  =  0.9ym. 

C-4. 2  Photodetection  Using  Photomagnetoelectric  and  Dember  Effects  in 
Gold-Doped  Silicon  13 

Photodetection  in  the  visible  and  near  IR  wavelength  region  using 

photomagnetoelectric  and  Dember  effects  in  gold-doped  silicon  has  been 

demonstrated.  Response  time  t  ■'f  0.2  Vis  and  noise  equivalent  power 
-10  r 

(NEP)  of  6  X  10  for  a  photomagnetoelectric  (PME)  cell,  and  t  <0.03 

-10  r 
ds  and  NEP  =  2.5  X  10  for  a  Dember  cell,  were  obtained  from  the 

present  results. 

14 

C-4. 3  Quantum  Yield  of  Metal-Semiconductor  Photodiodes 

A  theoretical  expresssion  for  the  quantum  yield  of  metal-semi- 
conductor  photodiodes  is  derived  that  takes  into  account  the  effects 
of  the  presence  of  an  inversion  region.  This  expression  predicts  a 
sharp  decrease  of  the  quantum  yield  in  the  short  wavelength  regime, 
in  contrast  with  previous  theoretical  predictions.  Applied  to  a 
Au-nSi  Schottky-barrier  photodiode,  our  theory  displays  the  functional 
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dependence  with  which  the  quantum  yield  increases  with  increasing 
doping  densities  for  short  wavelengths  and  decreases  with  increasing 
doping  densities  for  long  wavelengths.  Excellent  agreement  is  demon¬ 
strated  between  theoretical  predictions  and  experimental  data  in  the 
short-wavelength  regime  (i.e.,  for  0.4  <  A  <  0.6  ym)  . 

C-4.4  Quantum  Yield  of  P-I-N  Photodiodes  ^ 

An  expression  for  the  quantum  yield  of  p-i-n  photodiodes  is 
formulated  that  takes  into  account  the  effects  of  carrier  transport 
and  recombination  in  the  p-type  diffused  layer  and  consequently  describes 
the  decline  in  the  response  observed  for  short  wavelength  excitation. 
Applied  to  a  silicon  p-i-n  photodiode,  this  expression  yields  theo¬ 
retical  predictions  in  excellent  agreement  with  experimental  data 
over  the  entire  spectral  range. 

C-4.5  Analyses  of  Transient  Capacitance  Experiments  for  Au-GaAs  Schottky 
Barrier  Diodes  in  the  Presence  of  Deep  Impurities  16 

Transient  capacitance  experiments  have  been  performed  for  Au-n  type 

GaAs  Schottky  barrier  diodes  between  285  and  316°K.  The  GaAs  substrates 

are  doped  with  oxygen  and  chromium  respectively.  From  the  transient 

dark-capacitance  and  photocapacitance  measurements,  the  thermal  emission 

rate,  the  doping  densities  and  the  thermal  activation  energies  of  oxygen 

and  chromium  impurities  are  obtained  at  T  =  300°K. 

As  a  continuation  of  our  past  efforts  our  present  research  activities 
include;  Cl)  study  of  the  radiation  damage  in  bulk  GaAs  single  crystals 
by  using  the  methods  of  photomagnetoelectric  and  photoconductive  effects 
and  the  transient  capacitance  techniques,  C2)  study  of  electronic  trans¬ 
port  and  optical  properties  in  In  AsP  alloy  systems ;and  (3)  experimental 
and  theoretical  study  of  infrared  photodetectors  using  InAsP  alloys. 
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D.  NOISE  IN  SOLIDS  AND  IN  SOLID-STATE  DEVICES 


During  the  five  years  of  contract  support,  our  comprehensive  labor¬ 
atory  facilities  for  the  measurement  of  noise,  impedance  and  complete 
static  and  small-signal  frequencies  characteristics  have  enabled  research 
that  has  produced  many  useful  technical  findings.  As  an  integral  part 
of  our  overall  efforts  here  in  solid-state  electronics,  the  noise  research 
program  has  contributed  in  the  following  ways: 

(a)  As  a  tool  for  helping  evaluate  the  quality  of  newly- 
developed  devices.  This  capability  proved  valuable, 
for  example,  in  the  development  of  the  new  grating-type 
Schottky  barrier  photodiode  by  Li  and  Wang,  described 
in  Section  C. 

(b)  As  a  means  to  provide  insignt  not  obtainable  in  any 
other  way.  For  example,  study  of  noise  in  breakdown 
diodes  leads  to  a  detailed  understanding  of  the 
conduction  mechanisms  predominating  in  reverse-biased 
silicon  pn  junctions. 

(c)  As  a  means  for  corroborating  data,  reinforcing  the 
results  of  other  types  of  measurements  and  of  theoretical 
predictions. 

(d)  As  a  direct  measure  of  the  quality  of  a  solid-state 
device  acting  as  a  switch  and  of  its  sensitivity  as 
an  amplifier  or  a  detector. 

Our  research  program  has  yielded  results  in  each  of  these  various 
areas;  di  „criptions  of  the  most  important  of  these  follow. 

D-l  Low  Frequency  Noise  in  Bipolar  Junction  Transistors 

D-l.l  This  project  involved  the  experimental  study1  of  low  frequency 
noise  in  bipolar  junction  transistors.  In  addition  to  the  basic  under¬ 
standing  provided  by  the  research,  it  has,  from  a  practical  viewpoint, 

provided  a  characterization  of  the  transistor  that  is  useful  in  the 

2  3 

design  of  integrated  amplifiers.  ’  The  major  results  are:  (a)  the 
characterization  of  burst  noise,  demonstrating  that  the  source  of  burst 
noise  in  bipolar  junction  transistors  resides  predominantly  in  the 
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emitter-base  depletion  region  and  in  the  depletion  region  at  the  surface 
of  the  base  region  near  this  junction;  (b)  the  demonstration  that  the 
characterizaticnof  1/f  noise  requires  two,  rather  than  one,  generators 
in  the  equivalent  circuit  —  one  generator  associated  with  the  surface 
near  the  emitter-base  junction  and  the  second  with  the  active  base 
region;  (c)  the  development:  of  a  complete  equivalent  circuit  for  the 
low-frequency  noise;  (d)  the  explanation  of  structuring  of  the  low 
frequency  noise  spectra  observed  but  previously  unexplained. 

A  theoretical  model  that  describes  the  burst  noise  was  also  developed. 
This  model,  based  on  assuming  random  occupancy  of  flaw  centers  associated 
with  crystallographic  defects,  explains  the  characteristics  of  the  burst 
noise  seen  in  bipolar  transistors.  In  contrast  to  previous  explanations, 
this  new  model  takes  account  of  the  influence  of  defects  at  the  surface 
as  well  as  those  in  the  bulk.  For  many  transistors  the  defects  at  the 
surface  dominate  in  determining  the  observed  burst  noise.  Experimental 
results  on  specially  fabricated  transistors,  some  of  which  were  gold-doped, 
compare  well  with  the  predicted  behavior.  From  the  experimental  data  one 
can  infer  values  of  some  of  the  material  parameters  that  describe  the 
defects. 
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D~2  Noise  in  Phototransistors 

This  study  showed  the  utility  of  noise  measurements  in  determining 
the  parameters,  Beta  and  f^,  which  are  crucial  in  characterizing  a  photo¬ 
transistor.  For  transistors  fabricated  with  only  two  leads,  this  capability 
proves  of  special  value.  Detectivity  and  noise  equivalent  power  were 
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related,  providing  a  rare  example  in  which  the  noisf  characterization 
of  the  electrical  engineer  was  joined  to  the  sensitivity  measures  of  the 
detector  specialist. 
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D“3  Noise  in  Bipolar  Transistors  as  a  Function  of  Temperature1 

We  measured  the  noise  in  bipolar  transistors  at  room  temperature 
and  at  the  temperature  of  liquid  nitrogen  (77°K).  At  room  temperature 
excellent  agreement  between  theory  and  experiment  prevailed,  which  the 
literature  abundantly  suggests.  At  77°K,  however,  many  transistors  show 
excess  white  noise.  Our  studies  demonstrated  that  part  of  this  noise 
arises  from  generation-recombination  processes  within  the  device 
materials;  but  the  source  of  some  of  the  excess  noise  remains  unexplained. 
Research  presently  under  way  seeks  an  understanding  of  this  source. 
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D-^  Noise  in  Junction  Field-Effect  Transistors  and  in  MOS  Field-Effect 
Transistorsl  ~  ~  - - 

This  study  revealed  aspects  of  the  limitations  of  these  devices. 

From  noise  and  static  measurements  we  showed  that  leakage  occurring  in 

the  substrate  junction  of  an  MOS  field-effect  transistor  proved  to  be 

more  degrading  to  device  performance  than  the  gate  leakage  current  of 

junction  field-effect  transistors.  These  results  held  valid  at  high 

temperatures  for  all  applications  except  those  involving 

extremely  low  source  resistance.  The  study  gave  a  detailed  understanding 

of  the  effect  of  substrate  noise  current  upon  device  performance. 
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D-5  Noise  xn  Zener  and  Avalanche  Breakdown 

In  silicon  MOS  structures,  avalanche  breakdown  occurs  in  the  deple¬ 
tion  region  separating  the  electrical  channel  from  the  drain  contact. 

Noise  measurements  demonstrated  that  only  one  type  of  current  carrier 
predominates  in  causing  multiplication.1  Because  secondary  minority 
carriers  are  collected  at  the  substrate  terminal,  they  do  not  contribute 
to  the  multiplication  process.1 

Our  present  research  concerning  the  noise  properties  of  breakdown 
diodes  is  yielding  insight  concerning  the  conduction  mechanisms  in 
reverse-silicon  pn  junctions.  Here  both  electrons  and  holes  share  the 
responsibility  for  the  occurrence  of  multiplication.  Noise  measurements 
reveal  the  relative  contributions  of  Zener  aid  avalanche  currents. 
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D~6  Eff-ect  of  Radiation  on  Noise  in  Junction  Field-Effect  Transistors 
Damage  produced  by  even  moderate  radiation  tends  to  affect  severely 
the  low-frequency  noise  performance.  Measurement  of  the  noise  spectra 
as  a  function  of  temperature  yields  information  concerning  energy  levels 
associated  with  the  various  defect  centers  caused  by  the  radiation  damage. 
In  our  study  now  under  way,we  have  measured  the  noise  spectra  of  junction 
field-effect  transistors  before  and  after  irradiation  in  the  University 
of  Florida  nuclear  reactor.  Analysis  of  the  data  is  proceeding,  and 
the  insight  provided  will  be  described  in  the  next  Scientific  Report. 

D-7  Noise  in  Amorphous  Semiconductors1*^ 

The  band  structure  and  transport  properties  of  chalcogenide  semi¬ 
conductor  materials  have  attracted  interest  in  recent  years,  especially, 
following  the  announcement  of  new  electronic  devices  made  of  these 
materials.  By  comprehensive  electrical  and  optical  measurements,  combined 
with  measurements  of  noise  spectra,  we  contributed  to  the  understanding 
of  transport  in  these  materials.  Data  secured  both  from  optical  and 
electrical  measurements  were  consistent  with  the  prediction  of  a  theoretical 
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model  containing  a  reasonably  well-defined  energy  gap.  Noise  measure¬ 
ments  showed  a  strong  1/f  chracteristic,  which  is  expected:  lack  of 
structure  gives  rise  to  a  wide  distribution  of  time  constants. 
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E.  MAGNETIC  FILMS 

Part  of  the  research  under  this  contract  was  directed  at  digital 
memory  applications  of  magnetism,  a  subject  directly  relevant  to 
military  applications.  During  the  tenure  of  the  contract,  the  emphasis 
shifted  from  an  investigation  of  properties  relevant  to  random-access 
memory  with  thin  films  to  topics  related  to  sequentially-accessed  memory 
with  circular  magnetic  domains.  There  follows  a  synopsis  of  the  main 
results. 

E-i  Thin  Films 

We  were  interested  in  fast  switching  (reverse)  processes,  as 
particularly  relevant  to  memory  systems  under  development  elsewhere. 

At  the  initiation  of  the  study  the  physical  process  was  unclear  by  which 
fast  reversal  takes  place,  due  to  confusion  because  of  the  anomalously 
high  threshold  at  which  the  fast  performance  began.  We  identified^ 
threshold  conditions  in  a  consistent  way  and  related  them  to  measure¬ 
ments  of  dispersion. 

2 

An  improved  technique  was  developed  for  the  analysis  of  eddy- 
current  interactions:  particularly  those  interactions  between  film 
and  stripline  arising  from  a  transient  trapped-flux  torque.  By  a 
simple  R-C  circuit  approximation  to  the  magnetic  diffusion  equation, 
we  took  geometrical  factors  into  account. 


41 


REFERENCES 


1.  J.  K.  Watson,  "Effect  on  Magnetic  Film  Switching  of  Crystallite 
Interactions,"  IEEE  Trans.  Magnetics.  4,  722  (1968). 

2.  J.  K.  Watson  and  C.  T.  Leis,  "Analog  Approximations  of  Transient 

Flux  Diffusion,"  INTERMAG  Conference  paper  21,5.  Washington,  D.  C. 
(1970);  IEEE  Trans.  Magnetics.  6,  655  (1970);  IEEE  Trans.  Magnetics. 
2,  304  (1971).  “  - - 


III.  High  Current  Regimes  in  Transistor  Collector  Regions 
(D.  L.  Bowler  and  F.  A.  Lindholm)  ” 


1 .  INTRODUCTION 

A  number  of  authors  have  investigated  the  falloff  of  transistor  f 

T 

and  B  arising  from  high  current  in  the  epitaxial  collector  region.  The 
possibility  that  these  effects  are  due  to  base-widening  was  first  suggested 

1  2  o 

by  Messenger  and  first  published  by  Kirk  .  Beale  and  Slatter  have 
studied  devices  which  are  biased  in  such  a  way  that  they  appear,  at  the 
terminals,  to  be  in  the  active  region  but  in  which,  at  high  collector 
currents,  the  large  internal  voltage  drops  lead  to  saturation  of  the 
intrinsic  transistor.  Clark  ,  and  also  Hahn  ,  ,  have  discussed  the  dc 
characteristics  of  transistors  operating  in  this  same  mode. 

In  contrast,  van  der  Ziel  and  Agouridis^  have  suggested  the  possibility 
that  multi-dimensional  effects  are  important  and  have  proposed  a  model 
incorporating  lateral  injection.  Whittier  and  Tremere8  adopted  the  idea 
that  above  some  critical  current  space-charge-limited  flow  would  set  in 
and  then,  in  an  extensive  experimental  study,  attempted  to  determine  whether 
the  falloff  above  that  critical  current  behaved  more  like  the  van  der  Ziel 
and  Agouridis  model  or  more  like  the  Kirk  model. 

It  is  difficult  to  reconcile  these  various  suggestions  of  the  internal 
phenomena  responsible  for  observed  transistor  behavior  for  they  appear  to 
be  alternative  explanations  of  the  same  thing.  It  is  the  purpose  of  this 
paper  to  demonstrate  thay  they  are  not,  and  to  suggest  a  framework  within 
which  the  factors  governing  the  dominance  of  one  phenomenon  relative  to 
another  and  the  consequent  implications  on  the  functional  dependences  of 
transistor  gain  and  frequency  lesponse  are  revealed. 

In  particular,  we  will  show  that  the  high  current  behavior  of  devices 
having  long,  low  conductivity  epitaxial  collectors,  depends  on  the  applied 
collector-to-base  voltage.  For  voltages  lower  than  some  critical  value  the 
depletion  region  will  collapse  as  the  collector  current  is  increased  and  at 
some  critical  current  the  device  will  enter  either  the  saturation  or  quasi¬ 
saturation  mode  of  operation.  For  voltages  higher  than  the  critical  value, 
increasing  collector  current  will  cause  the  depletion  region  to  expand  until 
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the  collector  is  fully  depleted  (if  it  was  not  already  so  depleted  at  zero 
current  by  the  applied  bias)  after  which  further  increase  in  the  current  to 
some  critical  value  will  result  in  a  tendency  for  space-charge-limited  flow 
to  start.  To  explain  currents  higher  than  the  critical,  some  mechanism 
involving  multi-dimensional  flow  must  ve  invoked.  We  illustrate  the  useful¬ 
ness  of  these  conclusions  to  the  understanding  of  measured  device  performance 
by  several  examples  and  consider  their  implications  in  device  modeling  for 
computer-aided  circuit  analysis  and  design. 

2  .  THEORETICAL  DEVELOPMENT 

We  consider  a  one-dimensional  npn  transistor  in  which  the  uniform 
collector  doping  (N^)  is  much  less  than  the  base  doping.  This  assures  that 
the  depletion  region  is  essentially  all  on  the  collector  side  of  the  collector- 
base  junction.  We  assume  that  conduction  in  the  undepleted  portion  or  the 
collector  is  ohmic  and  that  hole  current  in  the  collector  is  negligible. 

We  begin  the  discussion  by  noting  that  there  are  two  situations  of  interest. 

In  one,  the  device  voltages  and  currents  are  such  that  the  collector  is  not 
fully  depleted.  In  the  other,  the  voltages  and  currents  are  such  that  full 
depletion  takes  place. 

A.  COLLECTOR  NOT  FULLY  DEPLETED 

In  Figure  1  we  show  the  situation  when  the  collector  is  not  fully  depleted. 
The  base-to-collector  metallurgical  junction  is  at  x=0,  the  collector  width 
s  W  and  xq  indicates  the  edge  of  the  depletion  region  when  the  collector 
current  is  zero.  The  plane  x ^  which  may  be  greater  or  smaller  than  xq, 
indicates  the  depletion  edge  in  the  presence  of  collector  current. 

If  the  collector-to-base  voltage,  V  is  held  constant  and  the  collector 

current  is  increased,  two  tendencies  will  develop.  One  will  tend  to  cause  the 

edge  of  the  depletion  region  to  expand  towand  W  ,  the  other  will  tend  to  cause 

c 

it  to  collapse  toward  x=o.  In  the  first  case,  as  the  collector  current 
increases,  the  electrons  which,  comprise  the  current  cause  a.  reduction  in  the 
net  charge  in  the  depletion  region  and  therefore,  in  order  to  maintain  the 
applied  potential,  the  depletion  region  will  tend  to  expand;  i.e,  x ^  will 
move  toward  Wc*  On  the  other  hand,  an  increasing  collector  current  gives 
rise  to  an  increasing  voltage  which  the  charge  in  the  depletion  region 
must  maintain  and  tends  to  cause  the  collapse  of  the  depletion  region;  i.e. 

*2  moves  toward  x=0.  We  wish  to  determine  under  what  conditions  each  of 
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these  tendencies  will  prevail  . 

B.  COLLECTOR  FULLY  DEPLETED 

If  the  tendency  of  the  depletion  region  to  expand  dominates,  the  situation 

ultimately  becomes  that  shown  in  Figure  2 (a.).  Here  we  show  the  case  in  which 

the  collector  region  is  fully  depleted,  either  by  the  applied  V  in  the 

cb 

absence  of  current  or  by  some  appropriate  comvination  of  collector  voltage 
and  collector  current.  Again  we  assume  that  Vcb  is  held  constant.  As  the 
collector  current  is  increased,  the  electrons  will  again  cause  a  reduction  in 
the  net  charge  density  in  the  collector  region  and  this  means  that  in  order  to 
support  the  applied  potential  the  electric  field  distribution  in  the  collector 
must  change. 

We  show  qualitatively  in  Figure  2(b)  how  the  electric  field  behaves  under 

these  conditions.  The  field  distribution  marked  (a)  corresponds  to  the 

collector  current  at  which  the  depletion  region  just  extends  to  W  .  As  the 

c 

collector  current  increases  above  this  value,  the  field  distribution  changes 

in  such  a  way  as  to  keep  the  area  under  the  curve  constant  since  V  is 

cb 

assumed  constant.  This  is  illustrated  by  curves  (b)-(d).  At  a  sufficiently 

high  value  of  collector  current  the  electric  field  at  the  collector-base 

junction  will  become  zero  as  indicated  by  curve  (e) .  We  call  the  current 

at,  which  this  occurs  Ic_jlcc>  the  onset  current  for  phenomena  associated 

with  high  collector  currents.  Normally,  for  the  transistor  to  carry 

currents  larger  than  this,  it  would  be  necessary  for  space-charge-limited 

flow  to  begin.  However,  in  the  usual  circuit  environment  where  there  are 

resistive  loads  in  the  collector,  and  also  in  the  situation  presently  under 

consideration  where  Vcb  is  constrained  to  remain  constant,  space-charge- 

Und-ted  flow  cannot  occur  so  other  mechanisms  must  account  for  increased 

currents.  We  will  refer  to  operation  at  currents  greater  than  I  as 

c-hcc 

operation  in  the  mode  associated  with  space-charge-limited  flow. 

C.  OPERATION  AT  LOW  COLLECTOR  VOLTAGE 

In  Figure  3  we  illustrate  qualitatively  common  emitter  characteristics 
typical  of  high  voltage  epitaxical  transistors  which  have  a  collector  doping 
substantially  lower  than  the  base  doping.  These  characteristics  show  what 
has  beem  described  as  two-region  saturation^.  At  the  lowest  values  of  V 

ce 

the  conductivity  of  the  collector  is  fully  modulated  by  the  Injection  of 
minority  carriers  from  the  base  and  the  device  displays  typically  low  satura¬ 
tion  resistance.  At  somewhat  higher  values  of  V  g  the  collector  conductivity 
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is  not  fully  modulated  -nd  the  device  is  in  partial  or  quasi-saturation. 

Hahn  has  shown  that  a  line  drawn  through  the  points  on  common  emitter 
characteristics  at  which  operation  changes  from  quasi-saturation  to  active 
region  will  be  straight  and  that  its  slope  corresponds  to  the  equilibrium 
resistance  of  the  undepleted  collector.  This  line  is  shown  as  I in  Figure  3. 
When  greater  than  V^e  we  have  V^cV  I2  can  be  described  by 


I,  =  q  N  y  V  ,  A  /Wc 

2  dc  n  cb  c  (1) 

where  the  symbols  have  their  usual  meanings.  Below  this  line  the  transistor 
is  operating  in  the  active  region  while  above  it  it  is  operating  in  full  or 
partial  saturation. 

Also  shown  on  Figure  3  is  the  current  I  .  This  current  is  defined  as 
the  current  that  would  flow  if  all  the  carriers  in  the  epitaxial  region 
moved  with  saturated  velocity.  Thus 
I.  =  q  N,  v  A 

1  dc  s  c  (2) 

where  vg  is  the  saturated  velocity  of  the  electrons  in  the  collector. 

We  define  a  voltage  as  the  collector-to-base  voltage  at  which 
1^  =  I 2*  Thus,  from  (1)  and  (2)  we  find 


V  =  W  —  . 

Icy 

n  (3) 

The  ratio  of  v  to  y  is  of  the  order  of  10^  volts/cm  so  the  order  of 

o  XI 

magnitude  of  is  about  one  volt  per  micron  of  epitaxial  layer  thickness. 

Consider  operation  of  the  transistor  at  constant  V  ,  (  V  )  less  than 

cb  ce 

along  an  operating  path  such  as  that  labeled  A  in  Figure  3.  When  the 
collector  current  is  zero,  the  collector-base  junction  is  reverse  biased  and 
a  depletion  region  extends  into  the  epitaxial  layer.  As  the  collector 
current  is  increased,  the  reverse  bias  at  the  intrinsic  junction  must  decrease 
until,  wt  -n  the  operation  crosses  I2,  the  intrinsic  junction  becomes  forward 
biased.  Thus  as  the  collector  current  is  increased  alone  A,  the  deoletion 
region  must  collapse  and  give  way  to  injection  of  minority  carriers  into  the 
collector  at  currents  greater  then  I^. 

An  analytic  treatment  of  the  dependence  of  depletion  region  width  on 
collector  current  is  given  in  the  Appendix.  A  simple,  general  result  cannot 
be  obrained  because  of  the  dependence  on  electric  field  of  the  velocity  of 
electrons  moving  through  the  depletion  region.  We  can,  however,  obtain 
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analytical  results  which  apply  to  the  limiting  cases  of  constant  (saturated) 
velocity  and  constant  mobility. 

In  the  constant  velocity  case  we  find 


1  -  Vr2 


1  1/2 


1  -  Vh 


Here,  x2  and  xq  were  defined  in  Figure  1,  1^  and  I 2 

(2),  and  I  is  the  collector  current.  We  note  that 
c 

as  Kirk's  result. 


(4) 

were  defined  in  (1)  and 
(4)  is  in  the  same  form 


I2  may  be  more  or  less  than  I  depending  on  the  device  makeup  and  applied 
voltage.  If  12  is  less  than  1^  Vcb  will  be  less  than  V  (c.f.  Figure  3), 
and  (4)  shows  that  as  the  collector  current  is  increa-ed  the  depletion  region 
will  collapse,  x2  finally  reaching  zero  at  the  critical  current  I2>  This 
confirms  the  qualitative  discussion  above.  In  obtaining  (4)  we  have  neglected 
diffusion  current  near  x2;  therefore  we  do  not  expect  it  to  be  exact  therw, 
nor  do  we  expect  it  to  describe  the  final  collapse  of  the  depletion  region 
in  detail. 

If  I2  is  greater  than  I ^  Vcb  will  be  greater  then  and  (4)  indicates 
that  the  depletion  region  width  will  increase  with  incicasing  collector 
current,  x2  utimately  reaching  W^.  We  discuss  this  situation  further  in 
part  D  of  this  section. 

In  the  constant  mobility  case  we  obtain  in  the  Appendix  the  following 
approximate  relationship  between  the  edge  of  the  depletion  region  and  the 
collector  current: 


x„  2 

r— 

if  I 

- 

2 

1  -  c 

i  „  i 

t 

1  +  B  .c 

A 

o 

L  J 

*2 

f°r  Ic  I2"  B  is  a  cons; -nt  of  the  order  of  0.25  and,  like  (4),  this  result 
neglects  diffusion  current  near  x2 . 

Equation  (5)  shows  that  as  collector  current  is  increased,  the  depletion 
region  width  will  decrease,  reaching  zero  when  the  critical  current  I2  is 
reached.  Referring  again  to  Figure  3  we  have  now  shown  that  in  the  region 
of  the  characteristics  to  the  left  of  V  ,  increasing  the  collector  current 
from1  zero  results  in  the  collepse  of  the  depletion  region  whether  constant 
velocity  or  constant  mobility  is  assumed.  There  is  no  reason  to  suppose  that 
the  behavior  would  be  otherwise  when  the  velocity  is  between  these  two 
extremes  and  so  we  assume  that  collapse  of  the  depletion  region  is  to  be 
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expected  when  the  applied  voltage  is  less  than  V^,  which  we  will  heareafter 

refer  to  as  V  .  . 

crit 

D.  OPERATION  AT  HIGH  COLLECTOR  VOLTAGES 

In  Figure  4  we  turn  our  attention  to  operation  at  higher  values  of 
applied  voltage  than  we  have  so  far  considered.  We  repeat  the  features 
shown  on  Figure  3  and,  in  addition,  now,  define  the  voltage  V2.  This  corres¬ 
ponds  to  the  V  ^  at  which,  in  the  absence  of  collector  current,  the  epitaxial 


layer  is  fully  depleted.  The  value  of  this  voltage  is: 

V,  1  q  N ,  w  2 12  ee  . 

2  n  dc  c  o 


12 


(6) 


In  order  to  make  it  possible  to  consider  distinct  regions  of  operation,  we 


restrict  our  attention  to  devices  for  which  V, 


V  This  amounts  to 

crit 


imposing  the  following  relationship  between  collector  doping  and  collector 
geometry : 

N^c  >  1.02  x  lO^/W^  (microns) 


(7) 


for  silicon  devices.  This  is  not  a  severe  restriction. 

In  the  preceding  analysis  we  have  shown  that  when  the  applied  bias  is 
greater  than  Vcrij_  and  electrons  move  through  the  depletion  region  at  vg  the 
depletion  region  expands  rather  than  collapses.  Since  Vcrit  is  the  applied 
bias  which  makes  the  average  field  across  the  collector  large  enough  to  cause 
the  drift  velocity  to  approach  its  saturated  value  it  is  reasonable  to  assume 
that  for  biases  larger  than  V  this  situation  will  occur.  Thus  if  we 

consider  operation  aolong  a  path  such  as  B  in  Figure  4,  the  depletion  region 
will  increase  in  width  as  the  current  is  increased  until  it  reaches  the  end 
of  the  epitaxial  layer.  Thure  will  be,  then,  a  locus  of  currents  at  which 
the  depletion  region  has  expanded  to  Such  a  locus  is  qualitatively 

illustrated  in  Figure  4. 

For  currents  higher  than  these  values,  the  electric  field  distribution 

in  the  collector  will  change  as  discussed  in  part  B  of  this  section.  We  can 

calculate  the  current  ^c_jlcc  at  which  space-charge-limited  flow  will  attempt 

to  start,  from  equation  (12)  of  the  Appendix  using  the  boundary  condition 

that  the  electric  field  must  be  zero  at  the  collector-base  junction.  Again, 

there  is  analytical  difficulty  due  to  the  field  dependence  of  the  mobility. 

However,  since  we  have  taken  the  view  that  for  biases  above  V  it  is 

crit 

a  reasonable  assumption  that  the  carrier  velocity  is  v  ,  we  make  that 

s 
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assumption  and  find 


I  ,  =  qv 

C-ncc  s 


2ee  V  ,/qW  ' 

O  CD  C 


+  N 


dc 


Ac  ■  h  +  h 


crit 


(8) 


Calculated  values  of  I  ,  are  shown  on  Figure  4  for  the  case  V  /V  =0.25. 

c-ncc  crit  2 

In  the  region  of  the  characteristics  where  the  applied  voltage  is  greater 
then  V2,  the  collector  is  fully  depleted  in  the  absence  of  collector  current 
and,  as  the  current  is  increased  from  zero,  there  is  no  expansion  phase  but 
the  process  described  above  sets  in  at  once  and  (8)  again  gives  I 

c-hcc . 


3.  DISCUSSION 

The  results  of  the  preceding  section  may  be  summarized  as  follows. 

If  the  applied  potentials  are  such  that  operation  is  below  V  .  the 

crit 

depletion  region  will  collapse  with  increasing  collector  current  and  at 
a  critical  current  given  by  (1)  will  enter  full  or  partial  saturation.  If 
the  applied  voltages  put  the  operating  point  above  V  ,  the  depletion 
region  will  expand  with  increasing  collector  current  until  it  reaches  the 
end  of  the  epitaxial  layer,  following  which  further  increases  in  current  will 
result  in  the  onset  of  the  mode  associated  with  space-charge-limited  flow. 

A.  INTERPRETATION  OF  DEVICE  PERFORMANCE  DATA 

To  illustrate  the  usefulness  of  this  description  in  assessing  the 

behavior  of  8  and  f^ ,  we  show  in  Figure  5a,b,c  results  which  are  typieal 

of  measurements  of  those  quantities  made  on  a  number  of  identically  processed 

devices*.  In  Figure  5d  we  show  the  calculated  values  of  V  .  .  I  and  I 

crit’  1  2 

for  these  devices,  as  well  as  the  loci  of  the  points  at  which  the  measurements 
were  made.  The  $  measurement  of  Figure  5a  was  made  along  path  A  of  Figure  5d 
and  since  the  collector  current  has  been  taken  well  above  the  critical  current 
I2  we  would  expect  to  see  the  falloff  that  does,  in  fact,  occur.  The  8  measure¬ 
ment  of  Figure  5b  was  made  along  path  B  of  Figure  5d.  In  this  case,  the 
collector  current  is  always  well  below  and  we  do  not  expect  any  substantial 
falloff  in  8  .  This  expectation  is  confirmed  by  the  measurements.  Hower13 
has  developed  an  expression  from  which  one  may  calculate  the  8  falloff  above 
the  critical  current.  His  development  is  based  on  a  one-dimensional  geometry. 

Points  calculated  using  his  approach  approach  are  shown  in  Figure  5a  and  are  seen  to 
agree  well  with  the  measurements.  Thus  the  8  falloff  for  this  device  appears 
to  be  satisfactorily  explained  on  the  basis  of  one-dimensional  effects  due  to 
high  currents  in  the  collector. 

*  These  devices  and  the  measurements  of  f  were  supplied  by  M.  Dickman, 

Tektronix,  Inc.,  Beaverton,  Oregon, 
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Further  confirmation  that  a  one-dimensional  analysis  suffices  for  this 
device  is  found  in  the  fact  that  the  value  of  the  collector  resistance 
calculated  using  the  emitter  mask  area  agree  with  the  collector  resistance 
determined  by  Hahn's  method  described  in  II, C.  Thus  emitter  and  collector 
areas  are  substantially  the  same  and  there  is  no  suggestion  of  lateral 
carrier  flow. 

Turning  now  to  the  data  shown  in  Figure  5c,  we  observe  that  f^  peaks  at 
collector  currents  of  about  50  ma,  and  has  dropped  appreciably  at  currents 
which  are  well  below  the  critical  current,  I ^  for  the  measurements  on  path  B, 
and  Ic_.hcc  (which  will  be  greater  than  I  )  for  the  measurements  on  path  C. 

Thus  we  conclude  that  the  F^  falloff  for  this  device  cannot  be  attributed 
to  the  effects  of  interest  here,  and  other  explanations  must  be  sought. 

G 

As  a  second  illustration,  we  refer  to  the  paper  by  Whittier  and  Tremere  . 

They  interpreted  their  results  to  mean  that  their  measurements  were  made  on 

devices  which  were  operating  in  what  we  have  called  the  mode  associated  with 

space-charge-limited  flow.  That  is,  V  ,  was  above  V  .  ,  and  collector 

cb  cnt 

current  was  above  Ic_hcc*  They  concluded  that  the  fT  falloff  they  observed 
was  deminated  by  lateral  injection  of  the  type  proposed  by  van  der  Ziel  and 
Agouridis^. 

In  Table  I  we  compare,  for  ".he  devices  used  by  Whittier  and  Tremere, 

Vcrit  as  calculated  by  (3)  and  the  V  ^  actually  used  in  their  measurements. 

It  will  be  seen  that  in  almost  all  cases  the  voltage  or  voltages  used  are 
well  below  V  ^  .  This  suggests  an  alternative  interpretation  of  their  data; 
namely,  that  their  measurements  were  made  almost  entirely  in  the  region  of 
operation  where  falloff  of  f^,  is  due  to  the  collector  current  exceeding  I^. 

One  should  then  expect  the  falloff  to  be  of  the  type  characterized  by  Beale 

3 

and  Slatter's  one-dimensional  analysis,  rather  than  of  a  type  which  is  to 
be  expected  when  1^ _hcc  exceec*ec*  anc*  multi-dimensional  effects  must  surely 

be  present. 


B.  IMPLICATIONS  FOR  DEVICE  MODELING 

Emphasis  has  recenlty  been  given  to  the  point  that  in  selecting  device 

models  for  computer-aided  circuit  simulation,  one  should  choose  from  those 

11  14 

available  the  simplest  one  that  will  give  the  required  accuracy  ’  .  In  the 

hierarchy  of  models  which  are  now  used,  the  one-dimensional  Ebers-Moll  model 
is  the  least  complex.  Also  one-dimensional,  but  incorporating  considerably 
more  detail,  is  the  Gummel-Poon^  model.  In  it,  an  approximation  to  one 
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multi-dimensional  effect  is  possible  through  the  value  assigned  to  a  parameter 
which  is  included  to  characterize  emitter  crowding.  The  fullest  account  of 
multi-dimensional  effects  can  be  taken  by  using  one  of  the  modified  forms  of 
the  Ebers-Moll  model  which  employ  tabulated  values  of  current  gain  and  charge 
control  time  constants  rather  than  single  numbers17. 

Our  results  suggest  that  if,  in  a  circuit  simulation,  one  determines 
that  a  particular  device  is  operating  in  the  quasi-saturation  mode,  then 
the  one-dimensional  Gummel-Poon  model  will  probably  give  adequate  results 
up  to  currents  and  gains  typical  of  usual  circuits.  On  the  other  hand, 
if  the  cirauit  simulation  shows  that  the  device  is  operating  in  the  region 
where  space-charge-limited  flow  may  be  important,  then  there  are  certainly 
multi-dimensional  effects  entering  and  it  would  be  necessare  to  use  one  of 
the  modified  Ebers-Moll  models. 
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MtetWS 


APPENDIX 


BEHAVIOR  OF  THE  DEPLETION  REGION  EDGE 

Referring  to  Figure  1,  we  assume  that  the  collector  is  depleted  for 
®  ^  x  £  x2  anc*  that  for  ^  x  <  there  is  ohmic  conduction.  We  write 

Poisson  s  equation  in  the  depletion  region 

~5x  =  ee“  (qNdc  '  VVnAc)  (9) 

o 

where  vn  is  the  effective  electron  velocity,  and  integrating  it  obtain  an 
expression  for  the  voltage  drop  across  the  depletion  region 

r*2 

Vcb'  J  £  dX  *  (10 

0 

* 

The  expression  for  the  voltage  drop  across  the  undepleted  part  of  the 
collector  is 


Vr  =  ~  x9)/qNj  y  A 

r  c  c  2  dc  n  c 


(ID 


We  then  make  use  of  the  fact  that  the  applied  collector-to-base  voltage 


v  u  ■  V  ,  ,  +  V 
cb  cb  r 


(12) 


to  obtain  a  solution  for  x^. 

A.  SATURATED  VELOCITY  CASE 

Letting  v^  —  v^ ,  which  implies  that  the  diffusion  current  is  everywhere 
negligible,  (9)  can  be  easing  integrated  subject  to  the  boundary  condition 
that  the  electric  field  at  x ^  must  be  the  field  required  by  the  ohmic  conducLioi 
in  the  undepleted  region.  That  is 

Ex=x2=  ’  VqNdcynAc-  (13) 
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From  (12)  we  then  obtain 


1/2 

1  -  Ic/l9 

1  -  Ic/IxJ  («) 

which  was  given  in  the  text.  and  I2  are  defined  in  (1)  and  (2). 

B.  CONSTANT  MOBILITY  CASE 

Letting  vn  =  and  proceeding  us  above,  we  find  that  it  is  not 
possible  to  obtain  a  meaningful  solution  by  using  the  boundary  condition  (13). 
We  attribute  this  problem  to  the  neglect  of  the  diffusion  current  which  may, 
in  this  case,  be  appreciable  in  the  vicinity  of  It  is  possible,  however, 

to  obtain  an  approximate  result. 

Instead  of  (13),  we  use  a  general  form  of  boundary  conditon;  namely, 

£x=x'  =  £'  (14) 


where  x'  is  an  arbitrary  point  at  which  the  field  is  e\  The  result  of 
Integrating  (9)  under  this  conditon  and  using  (12)  can  be  put  in  the  form 


(x 


e  -  e 


e 

o 


en 


e/e  +  I  n 

o  c  1 


e'/e  +  I  /I, 
o  c  2 


(15) 


which  relates  normalized  position  to  normalized  electric  field.  Here 


(x-x')N  -  (x-x')/l/2  xq2/Wc  ) 


E 

o 


/W  . 

c 


(16a) 

(16b) 


The  boundary  conditon  we  would  like  to  match,  given  by  (13)r,  can  be 
written 


(£/£o)x=x2  -  - 


(17) 


and  we  are  interested  in  cases  for  which  0  ^  ljl2  <:  1  (cf.  Figure  3).  Note 

that  for  an  npn  device  *  the  field  in  the  collector  will  be  negative  and  will 
have  its  maximum  magnitude  at  the  collector-base  junction. 
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In  Figure  6  we  have  used  (15)  to  plot  normalized  electric  field  vs. 

normalized  position  with  Ic/I2  as  parameter.  The  boundary  condition  (14) 

chosed  for  this  plot  is  e '/eQ  =  -2  at  x' .  Since  from  (17)  the  lowest  possible 

value  of  e'/eo  at  x2  is  -1,  this  chose  assures  that  x'  is  in  the  depletion 

region;  i.e.  to  the  left  of  x2  in  Figure  1.  We  have  also  plotted  in  Figure  t 

a  line  representing  the  value  of  the  field  in  the  undepleted  portion  of  the 

collector  (i.e.  the  boundary  condition)  for  the  choice  I  /I  =  0.5. 

c  2 

For  this  value  of  1^,  one  would  expect  to  get  a  straightforward;  graphical 
solution  from  the  intersection  of  the  line  representing  the  boundary  condition 
when  Ic/^I2  =  0,5  and  the  curve  of  position  vs.  electric  field  for  the  same 
value  of  Ic/I2  less  than  that  corresponding  to  the  boundary  condition  will 
make  intersections,  whereas  those  for  values  of  I  /l2  equal  to  or  greater 
than  that  corresponding  to  the  boundary  condition  will  not.  In  order  to  get 
an  approximate  result,  we  make  a  linear  approximation  to  the  family  of  field- 
position  curves.  This  assures  us  of  an  intersection  and  allows  us  to  calculate 
a  value  of  x2-  The  approximation  is  illustrated  by  the  dotted  extension  of 
the  curve  for  Ic/I2  =  0.5  on  the  figure. 

The  form  of  the  linear  approximation  to  the  f ield-sotision  curves  is 


e/e 

O 


e'/e  + 

O 


(x-x1)  N 
1  +  B  I 

c 


(18) 


where  B^r0.25,  and  exact  value  depending  on  e1 /eQ.  This  can  be  integrated 
and  combined  with  (12)  to  obtain  the  result  for  x2 

(x2/xq)2  (1-IC/I2)  (l+BIc/I2) 

which  was  given  in  the  text. 
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IV  •  Quantum  Yield  of  P-I-N  Phototiodes  (Sheng  S.  Li  and  F.  A.  Lindholm) 


1,  Introduction 

Most  previous  studies  of  the  spectral  dependence  of  photocurrent  and 
quantum  efficiency  of  silicon  photodiodes^”^  were  concerned  with  the 
photocurrent  and  the  spectral  response  of  the  quantum  efficiency  for  solar 
cells  operating  in  the  photovoltaic  mode.  The  present  paper  will  deal 
with  p-i-n  silicon  diodes  employed  as  photodetectors. 

The  operating  principles  of  a  p-n  junction  solar  cell  and  a  p-i-n 
photodetector  are  different  in  that  the  former  uses  a  low  resistivity 
substrate  while  the  latter  uses  a  high  resistivity  substrate.  The  maximum 
power  requirement  of  the  solar  cell,  dictates  that  the  resistivity  of  silicon 
substrate  be  as  low  as  possible^  to  minimize  resistive  power  loss  as 
current  is  drawn  from  the  device.  As  a  consequence  of  the  low  resistivity, 
the  depletion  layer  is  very  thin,  usually  of  the  order  of  100°A.  For 
wavelengths  of  light  of  interest,  the  absorption  of  radiation  takes 
place  over  distances  large  compared  to  the  width  of  the  depletion  layer 
to  be  separated.  The  contribution  of  photocurrent  in  a  solar  cell  thus 
comes  mainly  from  both  the  surface  diffused  layer  and  the  base  region 
of  the  p-n  diode.  The  contribution  from  the  depletion  layer  region  is 
usually  neglected. 

This  contrasts  with  the  behavior  of  reverse  biased  p—i— n  photodiodes, 
which  are  the  subject  of  this  paper.  Since  a  p-i-n  photodiode  operates  in 
the  depletion  mode  so  that  high  speed  and  sensitivity  can  be  achieved 
simutaneously,  the  depletion  layer  (i.e.,  the  "i"  region)  for  such  device 
is  usually  very  wide,  and  the  contribution  of  photocurrent  is  mainly  from 
the  "i"  region.  Therefore,'  the  existing  theory  for  solar  cells  will  not 
properly  describe  the  spectral  response  of  p-i-n  photodetectors. 

The  analysis  of  the  spectral  response  for  a  depletion  mode  p-n  diode 
was  first  reported  by  Gartner  and  subsequently  by  Jordan  and  Milnes."* 
Gartner's  treatment  ignored  completely  the  diffused  surface  layer;  and 
his  result  thus  failed  to  predict  correctly  the  shortwave  spectral 
response.  Jordan  and  Milnes  focused  their  attention  solely  upon  the 
transient  photoeffect  in  diffused  p-n  junction  diodes.  Recently,  McIntyre7 
and  Sprigings  et  al.  calculated  and  measured  the  quantum  efficiency  of 
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the  silicon  p-i-n  photodiode  with  special  attention  given  to  the  behavior 
at  long  wavelengths.  The  magnitude  and  drop  of  the  quantum  efficiency  at 
long  wavelengths  seems  to  be  well  explained  by  these  papers. 

Thus  expressions  for  the  quantum  yield  previously  derived  for  a 
depletion  mode  photodiode  fail  to  predict  correctly  the  spectral  dependence 
observed  in  metal-semiconductor  and  p-i-n  photodiodes  for  short  wavelength 
excitation.  The  short  wavelength  regime  is  defined  as  pertaining  to 
frequencies  of  incident  photons  for  which  &W  >>  1,  where  a  is  the  absorption 
coefficient  of  the  semiconductor,  and  W  is  the  depletion  layer  width. 

For  the  metal— semiconductor  photodiode,  Li  et  al.  have  recently  derived 
a  general  expression  for  the  quantum  yield  that  takes  into  account  the 
effect  of  an  inversion  layer  near  the  metal-semiconductor  interface.  Using 
this  expression,  they  were  able  to  explain  quantitatively  the  sharp 
decrease  of  the  quantum  yield  observed  in  the  short  wavelength  regime 
(i.e.  UV  regime)  for  the  Au-nSi  Schottky  barrier  photodiodes.  For  the  p-i-n 
photodiode,  William10  has  recently  reported  experimental  results  taken  on 
several  silicon  p-i-n  photodiodes  with  different  diffused  layer  widths. 

His  measurements  of  the  spectral  dependence  of  the  quantum  yield  suggest 
the  quantum  yield  at  the  short  wavelength  side  of  th  ;  spectrum  depends 
mainly  on  the  diffused  surface  layer  thickness.  But  no  theory  presently 
exists  to  explain  this  dependence. 

In  this  paper,  we  shall  derive  a  general  expression  for  the  quantum 
yield  of  a  p-i-n  photodiode  that  takes  into  account  the  effect  of  the 
diffused  surface  layer  thickness.  It  is  demonstrated  that  for  short 
wavelength  photons  this  component  predominates  in  determining  the 
photocurrent  observed  at  the  diode  terminals. 

The  resulting  expression  is  applied  to  silicon  p-i-n  photodiodes. 
Theoretical  predictions  and  experimental  data  are  then  compared  to 
demonstrate  the  accuracy  of  the  predicitons  over  the  whole  spectral 
range  of  interest. 
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2.  Derivation  of  Quantum  Yield  for  a  p-i-n  Photodiode 

Fig.  1  (a)  shows  a  schematic  representation  of  a  p-i-n  photodiode 
under  reverse  bias  conditions.  The  p-type  region  has  thickness  x 

o 

and  the  intrinsic  region  has  thickness  W,  with  W  x  .  Let  us  assume 

o 

that  a  monochromatic  radiation  imp  inges  upon  the  surface  at  x  =  0. 

The  rate  of  generation  of  photocarriers  is  given  by  11 


g(x)  =  $o  (l-R)ae 


-ax 


(1) 


as  is  illustrated  by  Fig.  1(b).  In  Eq.  (1),  $  is  the  photon  flux 

o 

density,  R  is  the  reflection  coefficient  at  the  surface  of  p-region,  and 
a  is  the  absorption  coefficient. 

Under  steady  state  conditions,  the  total  photocurrent  density  J 

ph 

produced  by  the  incident  photons  is  the  superposition  of  electron  and 
hole  current  components: 


=  J  +  J 

P 


This  applies  at  any  plane  x  in  the  photodiode.  Thus 


(2) 


(3) 

(A) 

is  the  component  of  hole  current  density  arising  from  carrier  generation 

in  the  intrinsic  region.  To  find  the  functional  dependence  of  the 

quantum  yield  U,  therefore,  we  need  only  express  J  (x  ) ,  J  (W) ,  and  J 

n  o  P  i 

as  functions  of  $  ,a  ,  W  and  x  . 

o  o 

The  derivation  of  these  expressions  is  based  on  the  following 
assumptions:  (1)  photo  inject  ion  is  small  in  both  ri-type  and  p-type 
regions  (i.e.  AP  <<  n^  and  An  <<  p  in  the  n— type  and  p— type  regions, 
respectively);  (2)  reverse-bias  is  small  enough  that  no  avalanche 
multiplication  occurs  in  the  intrinsic  region;  (3)  excess  carrier 
lifetime  is  extremely  short  at  the  illuminated  surface  so  that 
fin(0)  =  0  holds;  (4)  the  excess  carrier  density  at  the  border  of  the 


J  =  J  (x  )  +  J  (x) 
ph  n  o  p  o 

=  J  (x  )  +  J  (W)  +  J. 
no  p  i 


where 


J.  =  J  (x  )  -  J  (W) 

1  p  o  p 
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intrinsic  region  is  small  enough  that  An(xQ)  ~  0  holds;  (5)  the  effect 
of  electric  field  in  the  diffused  layer  (i.e. ,  p— region)  is  neglected; 
(6)  recombination  in  the  depletion  region  is  neglected. 

Assumption  (3)  is  valid  for  most  silicon  p-n  junction  photodiodes, 
since  the  carrier  lifetime  near  the  surface  of  the  diffused  layer  is 
usually  several  orders  of  magnitude  smaller  than  in  the  bulk.  2,10 
This  is  due  to  the  fact  that  the  surface  impurity  concentration  of  the 
diffused  layer  is  normally  several  orders  of  magnitude  higher  than  that 
in  the  bulk.  As  a  result  of  this  high  impurity  concentration,  the 
surface  lifetime  is  expected  to  be  much  shorter  than  in  the  bulk.  Thus 
carriers  generated  near  the  surface  will  recombine  before  they  are  able 
to  diffuse  toward  the  junction. 

For  completeness,  a  rigorous  treatment  which  includes  the  effect 
of  electric  field  and  arbitrary  surface  recombination  velocity  in  the 
diffused  layer  is  given  in  the  Appendix.  However,  the  simplified  model 
described  above  gives  agreement  within  ten  percent  between  the  present 
theory  and  the  experimental  data  reported  for  Si  p-i-n  photodiodes. 

We  shall  consider  separately  the  three  regions  that  contribute  to 
the  photocurrent,  as  follows: 
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(i)  p-type  region  (0  <  x  <  x^) 

In  this  region  the  contribution  to  the  photocurrent  is  mainly  due 
to  the  electron  diffusion  current  directed  toward  the  intrinsic  region 
near  x  =  xq.  The  component  is  obtained  by  solving  the  continuity 
equation 


n  9  An  An  , 

°n  9  2  ‘  T~  =  -<I,0(1-R)ae 

x  n 


-ax 


The  general  solution  of  Eq.  (5)  can  be  written  as 

x  -x 

~ )  +  B  cosh  i 
_n 


X  -X 

o 


An  *  A  sinh(~Y — )  +  B  cosh(^Jl)  -  $oa(1~R)Tn  -ax 

n  L 


2  3  e 
<“  Ln2-1) 

The  constants  A  and  B  in  Eq.  (6)  are  determined  from  the  boundary 
conditions  that  An  =  0  at  x  =  0  and  An  =  0  at  x  =  x  .  Thus 


<J)0a(i“R)Tt, 

2  2  f1“cosh(-~)e 

(uL^-l)s  inh  (~) 

Li 

n 


x  -ax 
o. 


n 


(5) 


(6) 


(7) 


B 


$oa(1“R)Tn  ~a^ 

3  3  “  e  L 

(a  Ln  -l) 


(8) 


The  electron  diffusion  current  density  at  x  =  x  ,  which  consists  of 

o 

electrons  directed  toward  the  intrinsic  region,  is  obtained  by  solving 
Eqs.  (6)  -  (8)  under  the  assumption  a Ln  »  1.  This  procedure  yields 


J  (x  )  =  qD 
n  o 


3An 
n  3x 


x=x 


=  q$o (1-R) 


-ax 


aL  sinh(— ) 
n  L 

n 


x  -ax 

- —  (1-cosh  (j-^-)e  °) 

o,  n 


(9) 


which  shows  the  functional  dependence  of  J  (x  )  on  the  absorption 

no 

coefficient  a  and  thickness  of  p-type  region  x  . 
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(ii)  intrinsic  region  (x^  <  x  <  W) 

The  drift  current  density  produced  by  carriers  generated  in  this 
region  is  given  by:'*' 

V7+x 

Ji  =  q  /  °  g(x)dx 

X 

o 

=  q$o(l-R)(e"aW-e  “S  (10) 

where  we  have  assumed  that  W  »  x  ,  so  that  W  +  x  k  W 

°  o 

(iii)  n-type  region,  (x  >  W) 

In  this  region,  holes,  the  minority  carriers,  flow  mainly  by  diffusion 

toward  the  reverse  bias  junction.  Thus  the  pertinent  continuity  equation 
is 

d  ^  -  -*  a-a)ae'“x  (11) 

p  3  2  T  o 

r  x  p 

where  and  tp  denote  the  hole  diffusion  constant  and  lifetime,  respectively. 
From  the  solution  of  this  equation,  one  finds  that1 


aL 

Jp(H)  ,  -(,*#)<i-r>  t_£7  .-«■ 

P 


(12) 


From  Eqs.  (9),  (10)  and  (12)  ,  we  obtain  the  total  photo-current  density 

_-ctW 


Jph  ■  ’V1-10 


and  the  quantum  yeild 


,  x  -ax 

(l-cosh(72-)e  °) - 


n 


sinhfc^)  (14oLp) 

n 


(13) 


n  = 


J  , 

jsh 

q$ 

o 

=  (1-R) 

i  x  -0£X 

(  )  (l-cosh (^)e  °)  - 


-aw 


(1+aL  ) 

P 


(14) 
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which  reduces  to  Gartner's  expression  of  n  if  x  and  R  are  set  equal 
zero. 

We  shall  now  use  Eq.  (14)  to  compute  the  quantum  yi.ald  of  a  silicon 
P-i-n  photodiode,  and  compare  the  results  with  experimental  data. 

3.  Results  and  Discussions 

For  a  practical  silicon  p-i-n  photodiode,  the  p— type  region  is 
usually  very  thin  (of  the  order  of  one  micron)  and  the  intrinsic  region 
(depletion  region)  is  relatively  wide.  Thus,  as  Fig.  1(b)  implies, 
the  absorption  of  short  wavelength  photons  gives  rise  to  the  generation 
of  excess  carriers  confined  mainly  within  the  p— type  material. 

In  calculating  the  quantum  yield  for  silicon  p-i— n  photodiodes, 
we  use  the  data  for  the  absorption  coefficient  of  silicon  given  by 
Dash  and  Newman  and  the  carrier  diffusion  length  deduced  from  the  life¬ 
time  data  of  silicon  given  by  Ross  and  Madigan.13  Fig.  2  shows  the 
quantum  yield  versus  photon  wavelength  for  silicon  p-i-n  photodiodes 
with  the  depletion  layer  width  W  as  a  parameter  and  with  x  =  1.5  ym. 

The  solid  lines  are  the  results  calculated  from  Eq.  (14)  with  W  =  0.32, 

0.08,  0.04  and  0.01  cm  respectively,  and  the  dash  lines  are  the  published 
data  given  by  Melchior  et.al.  The  reflection  coefficient,  R,  for 
silicon  has  been  assumed  constant  and  equal  to  0.3  for  the  present 
calculation. 

The  sharp  decrease  of  the  quantum  yield  on  the  short  wavelength  side 

of  the  spectrum  can  be  explained  as  follows.  The  carriers  generated  by 

photons  of  short  wavelength  appear  mainly  in  the  p-type  region,  as  implied 

by  Eq.  (5)  and  Fig.  1(b).  Thus  the  photocurrent  resulting  from  short 

wavelength  excitation  must  come  mainly  from  carriers  generated  in  the 

p-type  region.  However,  only  near  the  plane  x  =  x  does  diffusion  direct 
.  o 

electrons  toward  the  intrinsic  region.  Thus  only  a  small  fraction  of 

the  carriers  generated  by  short  wavelength  photons  in  the  p-type  region 

will  contribute  to  the  photocurrent  JR(x)  .  The  resulting  photoresponse 

contrasts  with  the  predictions  of  Gartner  that  the  quantum  yield  is  unity  in 

the  short  wavelength  regime.  The  functional  dependence  of  quantum  yield 

is  displayed  by  Eq.  (9). 

To  test  further  the  adequacy  of  the  present  theory,  we  measured  the 
quantum  yield  of  an  HP  5082-4207  silicon  p-i-n  photodiode,  and  compared 
the  result  with  out  theoretical  computation;  the  results  are  shown  in 
Fig.  3.  Again  excellent  agreement  is  obtained. 
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To  see  the  relationship  borne  between  thickness  ::  of  the  p-type 
region  and  quantum  yield,  we  computed  the  quantum  yield  versus  photon 
wavelength  for  W  =  0.01  cm,  and  xq  =  0.4,  0.8  and  1.2  ym,  respectively, 
using  material  constants  appropriate  for  silicon  p-i-n  photodiodes. 

The  result,  displayed  in  Figure  4,  shows  that  quantum  yield  will  increase 
with  decreasing  xq  in  the  short  wavelength  regime.  This  dependence  is 
in  harmony  with  the  measured  spectral  response  of  silicon  p-i-n  photo¬ 
diodes  reported  by  William.10 

Two  points  deserve  emphasis.  First  our  study  reveals  that  over  the 
whole  optical  spectrum  the  quantum  yield  of  a  silicon  p-i-n  photodiode  can  be 
improved  by  simultaneously  increasing  depletion  layer  width  W  and  decreasing 
the  junction  depth  Xq.  This  is  consistent  with  experimental  observations  of 
the  quantum  yield  of  silicon  p-i-n  photodiodes.  Second,  we  have  shown 
quantitatively  how  the  drop  in  quantum  yield  occurring  on  the  short  wavelength 
side  of  the  spectrum  depends  upon  the  diffusion  layer  thickness. 
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Appendix 


Spectral  Response  in  p-type  Surface  Region;  General  treatment 

The  continuity  equation  for  the  p-region  in  the  steady  state  is 


i  dJ  A 

1. _ n  An 

e  dx  x . 

n 


a<I>o(l-R)e 


-ax 


(A.l) 


where  £n  =  n-np  with  denoting  the  electron  density  in  thermal  equilibrium. 
The  electron  current  density  is 


J 

n 


eD 


dn 
n  dx 


eu  n 
Hn 


(A.  2) 


where  is  the  electron  mobility  which  is  related  to  the  diffusion 
coefficient  by  yn  =  0^)0^  e  is  the  built-in  electric  field  in  the  p-type  region 
Substituting  Jn  given  by  (A. 2)  into  (A.l)  and  approximating  e  by  a  suitable 
constant,  we  obtain 


d2An 


dx 


/£  x  1  dAn 

L  dx 
c  n 


An 


n 


-ad>Q  (l-R)e-ax 
D 

n 


where 


.  (iSl)  _i 

e°  6  Ln  is  defined  as  the  critical  field. 
The  general  solution  of  Eq.  (A. 13)  has  the  form: 


(A.  3) 


-x/L, 


An  =  A.^e 

a$>  (1-R)t 
o  n 


-x/L 


+  A^  e 


2  +  C  e'OT 


where  C  =  -  . . 

(a  l  2+aL  I— 1-1) 
of  Eq.  (A.3),  n  nV 


(A. 4) 

,  as  determined  from  the  particular  solution 

(A.  5) 


A^  and  are  constants  to  be  determined  from  the  boundary  conditions 


at  x  =  0  and  x  =  x  . 

o 


(A.  6) 
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The  boundary  contitions  at  x  =  0,  the  surface,  is  given  by 


Dn(dAn/dx)  =  An(pne+s) 


(A.  7) 


and,  at  x  =  xq,  the  junction, 

n  -  n  ~  o 
P  - 


(A. 8) 


where  s  is  the  surface  recombination  velocity. 

Before  applying  boundary  contitions  (A. 7)  and  (A. 8)  in  determining 
constants  A^  and  A£  from  (A. 4)  it  is  noted  that  for  practical  diffused 
p-n  junction  silicon  diodes,  the  electric  field  in  the  p-region,  | e.  I,  is 
usually  quite  high  (i.e.ej^,  2000  V/cm)^,  and  at  room  temperature 


is  less  than  a  few  V/cm.  Thus 

and  =  0. 

L2 


:m)  ^ ,  and  I  e  Jj 
>>i,  and  tA.( 


1  1  j  ®  I 

6)  reduces  to  - —  =  - —  — 
L,  L  I  £ 

1  n  c 


Combining  Eqs.  (A. 4)  through  (A. 9),  we  obtain  constants  A^  and  A^: 


1  o  “Ctx 

A1  =  C^L~  +  D~^e  °  ~  (a  +  1  1  s 

1  n 


Cm, 


2  o  “ax 

a2  *  c[(rr  +  f-)e  0 

1  n 


(A. 10) 


"Xo/Ll 


L,  D 


1  n 


=  Cm, 


(A. 11) 


where  C  is  given  by  Eq.  (A. 5),  and  and  m2  are  constants  representing 
bracket  in  (A. 10)  and  (A. 11). 

To  obtain  the  short-circuit  photocurrent  in  the  p-region,  we  use  Eqs. 

(A. 2)  and  A. 8)  and  evaluate  J  at  x  =  x.  which  gives 

n  0 

J  (x)  =  q  D 
n  n  dx 


®1  “x/L  _ 
=  -qD  C  [■— =■  e  0  1 

L1 


-ax 


+  ae  u] 


(A. 12) 


The  quantum  yield  in  the  p-region  is  thus  given  by 
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J  (x  ) 
n'  o 


q$ 


o2L2(1-R) 


[“V+  ‘“Ll 
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A-  /Vh 


-ax 


+  e  ]  (A. 13) 


Eq.  (A. 13)  is  the  general  expression  for  the  quantum  yield  in  the  p-region 
taking  into  account  the  effect  of  the  surface  recombination  and  the  built-in 
electric  field  in  this  region. 
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Fig.  1(a)  Schematic  representation  of  a  p-i-n  photodiode  under 
reverse  biased  condition.  Light  is  impinging  upon 
the  p-type  region  of  the  diode. 

(b)  Spatial  dependence  of  the  photocarrier  generation  rater 

A  denotes  long  wavelength  photons  and  A  short  wavelength 
photons.  ^ 
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Wavelength  (/xm) 

Fig.  2  The  quantum  yield  versus  photon  wavelength  for  different  depletion  layer  widths.  The  solid 
lines  are  comouted  from  (14)  and  the  dashed  lines  are  from  reference  6  for  silicon  p-i- 
photodiodes.  We  have  assumed  the  reflection  coefficient  R  =  0.3. 


Wavelength  (/xm) 


Fig.  3  The  quantum  yield  versus  photon  wavelength  for  x  =  0.4  ym  and 
for  W  =  lOym.  The  solid  line  was  calculated  from  Eq,  (14) 
and  the  dashed  line  displays  measured  values  for  an 
HP  5082-4207  silicon  p-i-n  photodiode. 
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Fig.  4  The  quantum  yield  versus  photon  wavelength  for  different  values  of  x  and 
and  for  W  =  0.01  cm  for  silicon  P-i-n  photodiodes. 


V.  Cluster  Model  Calculations  of  Silicon  Defects  (H.  F.  Schaake  and  L.  L.  Hench) 


Introduction 

That  the  properties  of  a  material  are  determined  by  the 
positions  and  types  of  the  atoms  in  it  has  been  known  virtually 
since  atomic  theory  was  first  enunciated.  In  two  areas,  the 
problem  of  a  periodic  solid  and  the  problem  of  a  molecule, 
substantial  progress  has  been  made  in  recent  years  in  develop¬ 
ing  tractable  theories  from  which  all  the  properties  of  the 
material  may  be  calculated.  These  theories  have  had  a  profound 
influence  on  progress  in  these  areas;  for  example,  the  entire 
semiconductor  technology  has  been  built  upon  band  theory. 

A  large  number  of  problems  do  not  fall  into  either  of 
these  two  areas;  these  problems  are  collectively  known  as  the 
defect  solid  state.  Some  of  these  problems  have  been  success¬ 
fully  treated  using  perturbation  techniques  and  the  band  theory 
of  periodic  solids.  However,  when  the  basic  solid  deviates 
greatly  from  a  periodic  structure,  such  as  in  an  amorphous 
solid,  or  when  the  defect  cannot  be  treated  as  a  perturbation, 
such  as  a  vacancy,  calculational  methods  have  not  been  avail¬ 
able. 

A  logical  starting  point  for  the  solution  of  defect  solid 
state  problems  is  from  either  of  the  two  existing  frameworks; 
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i.e.,  band  theory  or  molecular  theory.  In  a  band  theoretic 
approach,  a  suitably  small  number  of  atoms  encompassing  the 
defect  and  its  surroundings  would  be  taken  as  the  unit  cell  of 
a  solid;  periodic  boundary  conditions  would  repeat  this  cell 
throughout  all  space.  In  the  molecular  approach,  a  suitably 
small  number  of  atoms  encompassing  the  defect  is  treated  as  a 
large  molecule . 

In  this  work,  we  report  on  the  development  of  a  molecular 
approach  to  the  determination  of  the  properties  of  defect 
solids.  In  the  first  section  the  basic  method  is  developed 
(a  modification  to  the  semi-empirical  CNDO/2  method) anj 
justified  using  crystalline  defect-free  clusters.  In  the 
second  section,  the  method  is  applied  to  defects  in  silicon 
which  result  from  radiation  damage;  i.e.,  impurity-vacancy  com¬ 
plexes.  In  the  third  section,  the  method  is  applied  to  a  study 
of  amorphous  SiC>2  clusters. 

Molecular  Cluster  Method 

The  Schrodinger  equation  for  a  system  of  n  electrons  in 
the  potential  field  of  n*  nuclei  is 


H?  =  E<F 


(1) 


whe  re 


the  Hamiltonian  H  is 

n  vj  n' 
H  *=  Z  -  y1  -  I 
i=l  A= 1 


given  by  (in  atomic  units*) 


n  n 

Z  Z 

j  > i  i  =  1 


1 

r  .  . 
ij 


r 


ij 


(2) 


*Note:  Atomic  units  are  generally  used  in  atomic  calcula¬ 

tions.  This  system  is  defined  by  fi  =  h/2ir  =  m  =  e  =  1. 
Distances  are  measured  in  Bohr  radii  (i.e. ,  the  radius  of  the 
hydrogen  atom):  1  Bohr  radii  =  0.529  K.  Energies  are  mea¬ 
sured  in  Hartrees:  1  Hartree  =  27.21  eV  =  2  Rydbergs. 
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The  first  term  of  the  Hamiltonian  is  the  kinetic  energy 
of  the  electrons,  the  s  'ond  the  potential  energy  of  the  elec¬ 
trons  in  the  field  of  the  nuclei,  and  the  third  term  the  poten¬ 
tial  energy  of  the  interaction  of  the  electrons  with  each 
other.  The  wave funct ions  Y  have  the  following  restrictions 
placed  on  them. 

Normality  - 

/v*h3?1  d3?2  ...  d3;n  =  i  (3a> 

and  antisymmetr icity  - 


¥(...  r  j  ,  rJ+1  ...)  --  -*'( 


rj+l’  rj 


(3b) 


^  the  following,  no  attempt  will  be  made  to  derive  the 
equations  used  in  the  calculations.  These  derivations  are 
abundantly  available  in  the  literature.  Particularly  clear 
discussions  are  given  in  references  4  and  5. 


A .  Closed-shell  Systems 

In  general,  the  form  of  ¥  is  extremely  complex.  However, 
an  adequate  approximation  for  many  problems  is  to  consider  the 
simplest  analytic  form  satisfving  equations  (3a)  and  (3b)  above, 
i.e.,  a  single  determinant  of  one-electron  wave f un c t i ons : 


Y 


1 

/ir 


♦i<*i)  <f> i^2) 

^2(^l)  'P2(^2)  *2(*3) 
^3  (  r  1 )  <P3(r2)  <P3(r3) 


(4) 
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The  one-electron  functions  which  best  fit  equations  (1) 
and  (2 ;  (i.e.,  which  give  the  lowest  energy  E), provided  that 

there  are  an  equal  number  of  electrons  of  opposite  spins,  are 
the  solutions  to  the  Hartree— Fock  equations: 


F*i  - 


i  =  1  ,  ...  n 


where  the  Fock  operator  F  is  defined  as 


(5) 
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F  =  [H 


(5a) 


core 


+  ^  <2J,  -  K,)] 
j  J  J 


with  the  core  Hamiltonian 

P  2  e  A  rpA 

being  the  energy  of  an  electron  in  the  bare  nucleus;  the 
Coulomb  operator 


(5b) 


VV  ■  4  vV 


(5c) 


is  the  repulsion  energy  between  an  electron  at  r  ^  and  an  elec¬ 
tron  in  the  state  i/»  ,  while  the  exchange  operator  is  defined  as 

Kj(r1)^i(r1)  =  [/^(?2)  ^(r^d3^]^  (r^  (5d) 

and  arises  from  the  Pauli  exclusion  principle  requiring  the 
wavefunction  Y  to  be  antisymmetric  under  interchange  of  coordi¬ 
nates  (3b)  . 

The  solution  to  these  equations  is  also  very  complex,  and 
they  have  been  treated  only  for  extremely  simple  molecules. 

An  alternative  formulation  of  the  problem  is  to  expand  the  one- 
electron  wavef unctions  in  a  convenient  set  of  basis  functions, 
the  most  convenient  being  atomic  orbitals,  <{> 


*i(V  ■  2  %iV?i  -  V 


k=  1 

Since  the  AO's  are  fully  specified,  only  the  coefficients 
cui  need  be  determined.  For  the  lowest  energy  E  in  equation 

(1),  the  cyi  are  the  solutions  to  the  Ha r t r e e- Fo ck-Ro o t haan ^ 7 ^ 
equations  : 


2(F  -e.S  ) c  .  =  0 
v  MV  i  MV  vi  u 


where  the  Fock  matrix  is  given  by 


i  =  1  ,  .  .  .  n 


(6) 


76 


F  =  H 
y  v  y  v 


core 


+  Z  PAa  I Aa)  "  y(u A | va) 


A  a 


with  the  core  Hamiltonian  matrix 

n  1  7 

..core  .  ,  .  .  1  2  A  -+ 

\v  ■  ;*p<rl><-2  h  -  I  7THv(rl)d  rl 

A  1 A 

the  density  matrix 


(6a) 


(6b) 


PAa  =  2  ^  c Ai  cai 


the  AO  two-electron  integrals 


(6c) 


(Mv|Xo)  '  "  V?1>V?2>  7^  ■*,X(7i)'*,a^2)d3?id3?2  (6d) 

and  the  overlap  integral  matrix 


S 


yv 


( 6e  ) 


These  equations  are  algebraic  equations  of  the  third  order 
m  ;yi  and  are  completely  equivalent  to  the  Hartree-Fock  equa¬ 
tions  provided  the  basis  set  is  complete.  They  may  be 
solved  by  the  following  iterative  procedure:  assume  a  set  of 
cpi;  find  PAa  and  thus  Fy y ’  solve  equations  (6)  for  a  new  set 
of  cpi;  etc’»  until  the  new  set  of  Cyi  differs  little  from  the 
old  set.  Such  a  solution  is  called  a  self-consistent  field 
solution.  It  is  most  usual  to  iterate  until  the  total  elec¬ 
tronic  energy  E  changes  by  a  small  amount.  The  total  electronic 
energy  is  given  by 


E  =  2  E  tiCOre  +  E  E  (2 J . .  -  K  ) 
i  11  i  j  ^  iJ 


=  E 


y  A 


uXHmv  +  7  E  -  i  (uj|vo)] 


2 


(7) 


with  the  Coulomb  exchange  and  are  Hamiltonian  integrals 


J  .  . 
ij 


E 

y  Ava 


CyiCAjCviCaj(liVlACT> 


77 


£ 

yX  vio 


cmcxjcvic0j(l,xlvo) 


H 


co  re 
1  i 


MV 


c  c  HCOre 
Mi  vi  yv 


(9) 

(10) 


The  total  energy  of  the  system  relative  to  the  ion  cores 
and  free  valence  electrons  is  given  by 


ETotal 


E  +  £  £ 
A>  B 


ZaZ0R 
A  B  AB 


(ID 


i.e.,  the  electronic  energy  plus  the  core  repulsion  potential 
energy . 

The  resulting  one-electron  wave  functions  \ p  are  known  as 
molecular  orbitals.  (Note  that  in  these  equations  Greek  sub¬ 
scripts  refer  to  AO's,  while  the  Arabic  subscripts  always 
refer  to  MO ' s ! ) 

The  one-electron  eigenvalues  correspond  roughly  to  the 
energy  required  to  remove  an  electron  in  the  state  i  from  the 
molecule  (Koopman's  theorem). ^  They  are  given  by 


Hicr +  5  (2jd  -  *±J> 


i  =  1  ,  .  .  .  n 


(12) 


If  the  basis  set  contains  m>n  orbitals,  then,  in  addition  to 
the  n  occupied  orbitals,  there  will  be  (m-n)  unoccupied  orbi¬ 
tals  known  as  virtual  orbitals.  These  correspond  to  states  in 
which  an  additional  "test"  electron  is  added  to  the  molecule, 
and  hence  not  exactly  to  excited  states.  The  energies  of  the 
virtual  orbitals  may  be  corrected  to  yield  approximate  excita¬ 
tion  energies,  as  follows. 

The  ground-state  electronic  configuration  is  given  as 

2 


<P2  *3 


n 


i.e.,  one  in  which  each  of  the  lowest  lying  orbitals  are  occu¬ 
pied  by  one  electron  with  spin  a,  and  one  with  spin  B.  If  we 
look  at  an  excited  configuration 


2  2  2 
r2  ^3 


*n  ♦„+! 
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there  are  four  ways  of  coupling  the  spins:  act,  66,  a6  and  6a. 
These  configurations  give  rise  to  magnetic  quantum  numbers 
Ms  =  1,  -1,  0  and  0,  respectively.  The  first  two  and  a  combi¬ 
nation  of  the  second  two  belong  to  a  triplet  state,  S  =  1, 

Mg  ■  1,  0,  -1,  while  the  remaining  combination  of  the  second 

two  belong  to  a  singlet  state  S  =  0,  M  =0.  The  energy  of  the 

s 

triplet  state  above  the  ground  state  may  be  calculated  directly 
from  the  excited  configuration  act. 


AE„ 


e£  ‘ 


-  J 


Jin 


Z  >  n 


(13) 


The  coulomb  integral  arises  since  an  electron  has  been  removed 
from  molecular  orbital  n  and  therefore  should  appear  only 

once  in  equation  (12).  The  exchange  integral  K  in  equation 

Xj  XI 

(12)  is  correct,  since  the  electron  remaining  in  orbital  n  and 
the  electron  in  orbital  Z  have  the  same  spin,  and  therefore 
have  an  exchange  interaction. 

The  energy  of  the  configurations  a6  or  a6  above  the  ground 
state  will  be  given  by 


AE 

c 


e 

n 


Jin 


+  K 


Jin 


Z  >  n 


(14) 


The  coulomb  integral  arises  as  before  since  an  electron  has 

been  removed  from  state  n.  The  exchange  term  K  should  not 

Jin 

appear  in  equation  (12)  for  this  configuration,  since  the  elec¬ 
trons  in  states  Z  and  n  have  opposite  spin  and  hence  no  exchange 
interaction.  This  term  is  removed  by  adding  K  in  equation 

Xj  XI 

(14).  The  sum  of  the  configurational  energies  for  a6  and  6a 
will  he  equal  to  the  sum  of  the  singlet  and  triplet  energies; 
hence  for  the  singlet  state 


AE  =  2 AE  -  AE_ 
s  c  T 


=  eji  - 


n 


-  J 


Jin 


+  2  K 


Jin 


Z  >  n 


(15) 
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Note  that  since  both  J  and  K  are  positive,  the  singlet  state 
has  a  higher  energy  than  the  triplet  state,  in  accordance  with 
Hund ' s  rule .  ^  ^ 

1 •  Open  Shell  Systems 

In  the  foregoing  it  was  assumed  that  in  the  ground  state 
each  orbital  is  occupied  by  two  electrons  of  opposite  spin. 

In  a  large  number  of  problems,  such  as  impurity  states  in  semi¬ 
conductors,  dangling  bonds  in  amorphous  semiconductors,  and 
magnetically  ordered  materials,  this  is  not  the  case.  In  these 
cases,  since  the  number  of  electrons  with  a  spin  is  not  equal 
to  the  number  of  electrons  with  0  spin  (we  shall  consider 

na  >  n0^»  the  wavefuncti°ns  for  a  states  will  differ  from  the 
wave  functions  for  0  states. 

For  the  ground  state,  the  Hart ree-Fo ck  Roothan  equations 
become 


r  (f“  -  e?  S  )c“  =  0 

yu  i  yv  vi 

I  (F0  -  ej  S  )A  =  0 

^  yv  i  pv'  vi 


(16a) 

(16b) 


where  the  Fock  matrices  are  given  by 


yv 


Xa 


Xa 


Xo 


(17) 


and  a  similar  matrix  for  0.  These  equations  are  solved  in  a 
manner  similar  to  the  closed  shell  system:  An  initial  guess 
is  made  for  F^  and  F  ;  equations  (16)  are  solved  for  new 
coefficients  c^  and  and  so  on  until  self-consistency  is 

obtained . 

The  eigenvalues  are  given  by 


e“  =  +  l 


( J  a  \  j.  T  ot  0 

Uij  Kij}  +  Z  Jij 


(18) 


while  the  total  electronic  energy  is  given  by 
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na  ng 
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i  j  1J  i  j 


jaS  +  riB  Ta  jSot 
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a  act  ngng  oo 

-l  £  K?“  -  £  £  Kpp) 
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(19) 


The  number  of  electrons  in  the  system  is  equal  to  n  +n 
The  difference  between  these  two  occupation  numbers  may  be 
specified  in  terms  of  the  multiplicity: 


Multiplicity  =  2S+1 
nQ-ng  =  Multiplicity-1 

Excited  states  for  the  open  shell  calculations  will  be 
given  by 


(20a) 

(20b) 


AE£n  e£  "  en  +  J£n  *  K^n  H  >  n  (21) 

if  the  initial  orbital  n  and  final  orbital  £  have  the  same  spin, 
and 


AE£n  =  G£  "  en  +  J£n  £  >  "  (22) 

if  they  have  different  spins. 

2 •  Further  Comments  on  Excited  States 

The  way  in  which  we  treat  excited  states  is  approximate. 
This  problem  may  be  illustrated  most  strongly  by  considering 
a  closed  shell  system  to  be  calculated  by  the  open  shell  method. 
In  such  a  calculation,  the  mult iplicity=l ,  and  therefore  n  =n  . 
Under  these  circumstances,  the  two  Fock  matrices  will  be  iden¬ 
tical  at  selr-consistency,  e“  =  eE,  and  the  two  sets  of  eigen¬ 
vectors  cvi  will  be  identical.  Yet  we  will  not  get  the  singlet 
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and  triplet  excitations  as  derived  in  the  closed  shell  section 
The  problem  lies  in  treating  the  molecular  orbitals  as  the 
states  of  the  systems  rather  than  configurations.  Thus,  in 
the  problem  under  consideration  there  will  be  two  identical 
excited  configurations,  a&  and  Bet.  These  configurations  inter 
act  (i.e.,  <ip  6 1 H  |  ip  ^  0),  resulting  in  a  splitting  of  the 

levels  into  a  singlet  state  and  a  component  of  the  triplet 
state.  The  mathematical  process  is  known  as  a  configuration 
interaction  and  a  simple  form  of  it  was  used  to  derive  the 
singlet  and  triplet  levels  in  the  closed  shell  section. 

The  use  of  configuration  interaction  changes  the  Hartree- 
Fock  approximation.  The  single  de terminantal  approximation  to 
the  many  electron  state  function  ignores  correlation  which 
results  from  the  fact  that  the  charge  on  the  electrons  cause 
them  to  repel  each  other,  and  hence  the  probability  of  finding 
two  electrons  close  together  is  less  than  finding  them  far 
apart.  Letting  the  electrons  correlate  their  motion  reduces 
the  energy.  Correlations  may  be  introduced  into  the  many  body 
wavefunctions  by  adding  other  single-determinant  wavef unctions 
in  an  amount  to  minimize  the  total  electronic  energy. 

A  convenient  way  of  generating  these  additional  determin¬ 
ants  is  to  consider  excited  configurations  having  the  proper 
multiplicity.  Thus,  each  additional  determinant  has  one  row 
replaced.  Both  the  ground  state  and  the  excited  state  levels 
are  improved.  In  fact,  since  the  excited  wavefunctions  corres¬ 
pond  to  the  cluster  with  an  extra  electron  and  hence  are  too 
diffuse,  the  excited  energy  levels  will  be  improved  more  by  a 
configuration  interaction  than  by  the  ground  state. 

3 .  Semiempirical  Approximations 

The  above  framework  is  a  theory  in  which  calculations  at 
the  ab  initio  level  can  and  have  been  carried  out.  The  dis¬ 
advantage  of  such  a  procedure  is  the  large  number  of  integrals 
which  must  be  evaluated,  many  numerically,  and  the  size  of  the 
basis  set  which  must  be  used,  as  all  electrons  are  considered 
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in  the  calculation.  Such  a  procedure  is  time  consuming,  and 
even  for  small  clusters  makes  extensive  demands  on  computer 
storage.  Since  during  the  course  of  study  of  a  defect  problem, 
many  defect  clusters  must  be  considered,  the  cost  of  ab  initio 
calculations  makes  them  unattractive. 

A  potential  way  around  this  quandry  is  to  adopt  semi- 
empirical  techniques.  That  is,  rather  than  attempt  to  evalu¬ 
ate  every  integral  and  use  a  large  basis  set,  many  of  the 
integrals  which  are  difficult  to  evaluate  and  which  are  known 
to  be  small  are  discarded,  others  are  evaluated  from  empirical 
data,  and  the  basis  set  is  generally  limited  to  the  valence 
atomic  orbitals.  Such  an  approach  definitely  has  its  short¬ 
comings;  to  name  two:  1)  a  particular  set  of  assumptions  may 
work  well  on  one  cluster,  but  fail  miserably  on  another,  and 
2)  a  set  of  assumptions  may  be  adjusted  to  give  good  proper¬ 
ties  one  type  and  strongly  mispredict  another.  On  the  other 
hand,  the  very  nature  of  the  way  in  which  defect  solid  state 
problems  must  be  studied,  that  is,  calculations  must  be  made 
on  a  large  number  of  similar  clusters,  each  one  differing  only 
slightly  from  another,  suggests  that  a  semiempirical  technique 
which  basically  follows  the  framework  of  an  ab  initio  calcula¬ 
tion  may  successfully  predict  differences  between  clusters  once 
the  parameterization  has  been  adjusted  to  put  the  property  of 
interest  in  the  right  ball  park.  This  is  the  approach  taken 
in  this  work,  and  the  remainder  of  this  section,  as  well  as 
other  sections  indicates  that  this  is  a  not  unreasonable 
assumption . 

A  wide  variety  of  semiempirical  techniques  have  been 
developed  in  the  past  thirty  years.  However,  one  that  fits 
our  requirements  for  a  theory  which  follows  the  framework  of 
an  ab  initio  calculation  and  which  has  been  successfully  applied 
to  a  large  range  of  molecular  problems  is  the  CNDO/2  method  of 
Pople,  Sentry  and  Segal, ^  ^  and  is  the  theory  which  we  have 
adopted.  For  a  full  account  of  the  development  of  the  approxi¬ 
mations  and  the  parameterization,  the  reader  is  referred  to 
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the  references  cited  above;  here,  we  shall  merely  enumerate 
the  approximations  made  in  the  CNDO/2  method. 

(1)  Complete  Neglect  of  Differential  Overlap  (CNDO) . 
All  two-electron  atomic  orbital  integrals  are  set  equal  to 
zero  except  for  those  of  the  type  (p y | XX) .  These  integrals 
are  further  reduced  to  one  per  pair  of  atoms  (or  per  atom): 
Clip  I  XX)  =  YAB4>y  on  A,  <j>^  on  B. 

(2)  Replacing  the  overlap  matrix  with  the  identity  in 
the  Hartree-Fock-Roo thaan  equations,  and  neglecting  overlap 
in  normalizing  the  molecular  orbitals. 

(3)  The  diagonal  elements  of  the  core  Hamiltonians  are 
approximated  as 


H 


core 

yy 


-  2  <VV  +  l  VaB 


on  A . 


1^  is  the  ionization  potential  for  orbital  y  on  atom  A. 

A^  is  the  electron  affinity  for  orbital  y  on  atom  A. 

is  the  core  charge  (atom-valence  electrons)  of  atom  B. 
(A)  Off  diagonal  elements  of  the  core  Hamiltonian  are 
set  equal  to 


H 


core 

yv 


6AB  Syv 


where  sab  ■  y  «;+•;> 


<j>  on  A,  <b  on  B. 
y  v 


unless  either  A  or  B  (or  both)  are  second  row  elements,  whence 


6ab  ■  f  <»;♦>!> 


^A  an<*  ^B  are  emPiricai  parameters  adjusted  to  give  good 
results . 

(5)  Y AB  =  <nSAnSA I ri2 I nSBnSB>  *  where  nS^  is  the  valence 
S  orbital  of  atom  A,  and  nSg  the  valence  S  orbital  of  atom  B. 
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For  hydrogen,  the  valence  basis  set  consists  of  one  S 
orbital;  for  first-row  elements  one  S  and  three  p  orbitals, 
and  for  second-row  elements  one  S,  three  p  and  five  d  orbitals. 
All  orbitals  are  Slater-type  orbitals  (STO)  for  calculation 
simplicity . 

Using  these  approximations,  the  Fock  matrix  may  be  writ¬ 
ten  as 


MM 


-  7  .+*,.)  +  .)Y 


AA  2  MM  AA 


+  bJa(pbbyab‘zbyab) 


FPV  ■  7  kab<s:+®b>sUv  -  7  pm^ab 


M  ^  V 


(23a) 

(23b) 


^Ag  =  1  or  3/4  as  explained  in  (4);  is  on  atom  A,  <j>  on 

atom  B;  P  =  IP... 

A A  A  on  A  A \ 

The  self-consistency  is  started  by  an  initial  Huckel-type 
approximation  by  placing  YAA=YAB“0  for  all  atoms.  This  initial 
Fock  matrix  is  diagonalized,  a  density  matrix  constructed,  the 
total  electronic  energy  found,  a  new  Fock  matrix  generated 
(from  now  on  YAB»YAA  have  their  calculated  values),  etc.,  until 
the  energy  changes  by  less  than  a  predetermined  amount.  This 
procedure  may  be  executed  by  the  program  listed  in  the  Appendix.* 
The  extension  to  open  shells  is  straightforward  and  will  not  be 
considered  here.  The  program  in  the  Appendix  will  also  perform 
open  shell  calculations. 


4 •  Cluster  Boundary  Conditions 

Early  in  this  investigation  it  was  found  that  a  cluster 
of  17  carbon  atoms  in  a  crystalline  arrangement  did  not  have  a 
well-defined  energy  gap  --  levels  were  spread  throughout  a 
rather  wide  range.  This  problem  resulted  from  the  fact  that 
the  coordination  was  improper;  instead  of  a  coordination  of 
four,  the  surface  atoms  had  a  coordination  of  one.  The  result 
is  an  extremely  large  concentration  of  "dangling  bonds"  —  a 

♦Available  upon  request  from  the  Dept,  of  Materials  Science 
and  Engineering,  University  of  Florida,  Gainesville,  Florida  32601, 
c/o  L.  L.  Bench. 
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feature  not  typical  of  the  solid.  This  problem  has  apparently 
been  recognized  by  several  other  workers ^10"13^  who  have 
derived  several  ways  around  them:  saturating  the  surface  with 
hydrogen,  imposing  periodic  boundary  conditions  or  partial 
periodic  boundary  conditions,  or  immersing  the  cluster  in  a 
potential  well.  These  boundary  conditions  have  various  weak¬ 
nesses.  The  surface  hydrogen  atoms,  while  resulting  in  a 
physically  plausible  molecule,  introduce  a  new  chemical  species 
into  the  system,  and  under  certain  circumstances,  the  surface 
requires  two  hydrogen  atoms  to  be  separated  by  about  1  Angstrom 
a  muc.i  smaller  distance  than  found  in  any  molecule.  Addition¬ 
ally,  the  17  atom  carbom  cluster  requires  36  surface  hydrogens 
for  saturation,  resulting  in  the  C-H  and  H-H  interactions 
assuming  an  extremely  important  role  in  determining  the  elec¬ 
tronic  structure.  Periodic  boundary  conditions  are  perhaps  the 
best  conditions  since  the  cluster  is  then  surrounded  by  an 
array  of  atoms  similar  to  that  which  is  found  in  the  solid. 
However,  the  geometrical  constraints  placed  on  the  cluster  in 
order  to  have  translational  periodicity  are  rather  severe,  and 
generally  require  a  large  number  of  atoms. 

An  alternative  approach  may  be  formulated  by  considering 
the  problem  of  the  surface  in  terms  of  the  Fock  matrix.  Con¬ 
sider  first  the  cluster  embedded  in  the  solid  surrounded  by 
all  of  the  other  atoms  in  the  solid.  The  Fock  matrix  will  con¬ 
tain  off  diagonal  elements  between  all  orbitals.  If  a  trans¬ 
formation  can  be  found  which  minimizes  off-diagonal  elements 
between  orbitals  on  atoms  designated  as  surface  atoms  of  the 
cluster,  and  orbitals  on  atoms  outside  the  cluster,  then  the 
Fock  matrix  can  be  reduced  to  the  following  form: 


Under  these  circumstances,  the  cluster  Fock  matrix  is  sepa¬ 
rated  from  the  remainder  of  the  total  Fock  matrix  and  may  be 
diagonalized  by  itself  to  yield  the  behavior  of  the  wavefunc- 
tions  within  the  cluster.  Such  a  transformation  is  in  general 
impossible  in  ionic  and  covalent  materials;  however,  a  trans¬ 
formation  to  hybridized  surface  orbitals  comes  surprisingly 
close . 


First  of  all,  it  should  be  realized  that,  in  typical 
solids,  the  magnitude  of  off-diagonal  density  matrix  elements, 

Pnv’  are  one  to  two  orders  of  magnitude  smaller  for  orbitals 
on  second  neighbor  atoms  than  for  orbitals  on  nearest  neighbor 
atoms.  Furthermore,  the  magnitude  of  these  elements  decreases 
very  rapidly  as  the  atoms  to  which  the  orbitals  are  attached 
increase  further.  If  the  basis  is  transformed  to  a  hybridized 
basis,  the  overlap  integrals  are  substantially  larger  for  near¬ 
est  neighbor  atoms  than  for  more  distant  neighbors.  In  dia¬ 
mond,  for  example,  with  an  sp3  hybrid  basis,  the  overlap  inte¬ 
gral  with  the  nearest  neighbor  atom  to  which  the  hybrid  is 
oriented  is  0.94,  to  other  nearest  neighbor  atoms,  0.01,  and  to 
second  neighbor  atoms,  0.15.  This  means  that  in  a  hybridized 
basis  there  is  one  large  off-diagonal  element,  and  the  rest 
are  fairly  small. 

This  transformation  to  a  hybrid  basis  then  approximately 
satisfies  the  requirements  set  forth  previously,  and  we  can 
anticipate  fairly  good  results  if  we  terminate  the  Fock  matrix, 
keeping  orbitals  with  large  off-diagonal  elements;  i.e.,  the 
hybridized  orbitals  on  the  surface  atoms  which  are  oriented 
toward  the  interior  atoms  of  the  cluster. 

The  use  of  surface  hybrids  has  a  firm  intuitive  base. 
Chemical  saturation  of  the  system  is  maintained  by  placing  an 
atom  with  one  (or  more)  hybrids  on  the  surface.  The  hybrid 
forms  a  bond  with  an  interior  atom,  but  since  there  are  no 
hybrids  oriented  to  the  exterior  of  the  cluster,  there  are  no 
dangling  bonds."  This  approach  requires  that  we  add  an  elec¬ 
tron  for  every  hybrid,  and  that  in  order  to  maintain  charge 
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neutrality  (i.e.,  in  a  purely  covalent  structure,  the  charge 

on  an  ion  is  0)  the  core  charge  must  also  be  set  equal  to  the 

number  of  hybrids.  This  approach  will  work  with  Column  IV 

elements  of  the  periodic  table,  C  and  Si,  which  normally  have 

3 

tetrahedral  coordination  and  sp  hybridization.  Only  these 
elements  have  been  used  for  surface  hybrids  in  this  work. 

The  program  in  the  Appendix  is  capable  of  performing 
calculations  as  described  above.  Only  surface  atoms  are  hybrid¬ 
ized,  leading  to  a  mixed  basis  set.  The  diagonal  portion  of 
the  core  Hamiltonian  must  be  entered  externally  — —  for  an  sp 

hybrid  it  has  the  value  (-7HC°re  +  yHCOre). 

4  ss  4  pp 

5 .  Parameterization  for  Energy  Level  Diagrams 

While  the  CNDO/2  method  leads  to  good  values  for  bond 

lengths  and  charge  distributions  in  molecules,  it  is  known  to 

predict  a  value  for  the  optical  absorption  energy  that  is  too 

large.  For  example,  the  molecule  neo-pentane  (C,.H12,  central 

C  atom  plus  4  methyl  groups  in  tetrahedral  coordination)  is 

calculated  to  have  a  lowest  excitation  energy  of  about  15  eV 

in  comparison  with  the  experimental  value  of  7.5  eV.  This 

problem  has  been  examined  by  a  number  of  investigators,  most 

notably  Jaffe  and  DelBene  ,  (14-16 )  and  Salahub  and  Sandorfy .  ^ 7'20^ 

They  have  introduced  modifications,  the  chief  ones  being  to 

reduce  the  electron-electron  interactions  by  reducing  y  and 

to  alter  the  values  of  3°.  In  the  calculation  of  Energy  Level 

Diagrams  we  have  used  an  extension  of  the  method  which  Sandorfy 

and  Salahub  used  for  saturated  hydrocarbons,  with  the  value  of 

6°  adjusted  to  give  reasonable  energy  gaps.  The  electron- 

electron  interaction  integrals  y  are  evaluated  by  the  method 

(21-22)  Aa 

of  Pariser  and  Parr'  rather  than  by  the  CNDO/2  method. 

Furthermore,  configuration  interactions  are  not  performed; 
this  is  a  rather  time-consuming  procedure.  The  excitation 
energy  is  taken  as  the  configurational  energy,  equation  (14). 

This  leads  to  uncertainties  in  the  calculated  energy  gaps  of 
about  x  j  or  about  0.05-0.2  eV,  depending  on  the  size  of 
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the  cluster  (as  the  cluster  increases  in  size,  K . .  decreases). 
Considering  the  approximate  nature  of  the  calculation,  this  is 
considered  to  be  an  acceptable  uncertainty. 

6 .  Population  Analysis 

Having  cast  the  problem  in  a  Linear  Combination  of  Atomic 

Orbital  (LCAO)  framework,  it  is  extremely  easy  to  get  extensive 

information  on  the  chemical  bonding  in  the  cluster.  This 

information  is  most  easily  obtained  from  a  population  analysis, 

(2  3) 

introduced  by  Mulliken.  Consider  the  LCAO  approach  to  the 

hydrogen  atom.  The  ground  state  MO  will  be 


{C1*1  +  c2*2} 


The  electron  density  as  a  function  of  position  for  a  single 
electron  will  be  given  by 

IU2  *  cl*l  +  (C1C2  +  C2C1)<^,JL^2  +  C2*2 

This  distribution  may  be  given  the  following  interpretation: 

2  2  2  2 
c^<J>2  *s  t*ie  charge  distribution  arising  from  atom  1,  c2^2 

from  atom  2,  and  (c^c2+c2cl^i^2  t*ie  charge  distribution  aris¬ 
ing  from  the  covalent  bond.  Note  that  if  the  two  atoms  are 
separated  by  an  infinite  distance,  <J)  1 4>  2  =  0 ;  i.e.,  there  is  no 
bond,  as  expected.  Integrating  each  of  these  terms,  we  come 
to  the  following  interpretation,  since  /^Qdt =1 ,  and  remembering 
that  two  electrons  occupy  the  orbital: 


2^sll 

is 

the 

charge  on 

atom 

1, 

2C2S22 

is 

the 

charge  on 

atom 

2  ,  and 

4cl  c2 

S12 

is 

the  charge 

in  the  bond 

We  can  therefore  define  the  orbital  population  matrix  as 
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(24) 


B  =  P  S 
yv  yv  yv 

and  the  atom  bond  population  as 


AB 


E 

y  on  A 


v  on  B 


P  S 
y  v  y  v 


(25) 


®AA  rePresents  the  charge  on  atom  A  and  2BAfi  the  charge  in  the 
bond  between  atoms  A  and  B.  If  B^  is  negative,  then  the  bond 
is  anti-bonding. 

For  open  shell  calculations 


yv 


yv 


+  P' 


yv 


We  can  also  define  a  Spin  Population  matrix 


(26) 


B 


s 

yv 


I 

y  on  A 


E 

v  on  B 


(27) 


which  tells  where  the  electrons  contributing  to  the  excess 
spin  are  distributed. 

In  the  CNDO  method  the  molecular  orbitals  are  not  normal¬ 
ized.  Therefore,  the  populations  do  not  correspond  exactly  to 
charges.  Nevertheless,  they  do  qualitatively  indicate  the 
distribution  of  charge. 


7 .  Results  on  Crystalline  Clusters 

We  are  now  in  a  position  to  discuss  the  results  of  calcu¬ 
lations  on  crystalline  clusters.  The  purpose  of  this  section 
is  to  demonstrate  that  the  modified  CNDO  method  as  outlined 
above  does  in  fact  yield  reasonable  results  for  solids.  Later 
sections  will  treat  specific  defect  problems. 

a •  Diamond  and  carbon  structures 

A  large  amount  of  work  has  been  done  on  carbon-based  mole¬ 
cules;  this  work  forms  a  convenient  starting  point. 

Consider  first  the  molecule  ethane  (Figure  1).  Using  the 
Sandorfy  and  Salahub  hydrocarbon  parameterization  (Table  1)  ,  (^“^0) 
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the  singlet  excitation  energy  is  found  to  be  9.01  eV,  versus 
an  experimental  value  of  9.4  eV.  If  we  remove  one  hydrogen 
atom  and  the  associated  sp3  hybrid  on  the  carbon  atom,  the 
singlet  excitation  increases  slightly  to  9.19  eV.  Going  even 
further  and  plucking  off  three  hydrogen  atoms  surrounding  a 
single  carbon  atom  and  removing  the  three  associated  sp3  hy¬ 
brids  yields  a  singlet  excitation  energy  of  10.48  eV  (the 
singlet  excitation  energy  of  methane  CH^  by  comparison  is 
9.7  eV). 


Table  1 

Sandorfy  and  Salahub  Parameterization 
for  Saturated  Hydrocarbons 


^H-H  “  ”9,5  eV 

0C-H  "  -11*0 
8C-C  *  ’15*° 

Y  C  *  11*1 


This  trend  of  increasing  excitation  energy  with  decreasing 
cluster  size  is  well  known  as  a  chemical  shift.  The  bond 
populations  for  these  three  clusters  is  illustrated  in  Table  2. 
(In  this  table  and  in  all  remaining  tables,  atoms  having  less 
than  the  full  complement  of  orbitals  are  designated  by  a 
prime.)  It  will  be  noticed  that  removing  atoms  has  little 
effect  on  the  bonds  formed  with  other  atoms.  The  charge  dis¬ 
tributions  are  also  little  affected. 

The  effect  of  surface  hybridized  orbitals  can  be  compared 
to  surface  saturation  by  hydrogen  by  considering  the  molecule 
neo-pentane  and  the  cluster  C^-C^*  These  clusters  consist  of 
a  central  carbon  atom  and  four  carbon  atoms  in  tetrahedral 
coordination.  The  second  neighbor  coordination  shell  consists 
of  12  hydrogen  atoms  for  neo-pentane  or  12  C'  atoms  (in  <220> 
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type  locations),  each  containing  a  single  hybrid  (Figure  2). 

The  energy  gap  is  found  to  be  8.35  eV  for  neo-pentane  and  8.99 
eV  for  C ^ C^2  (Figure  3) .  Nearest  neighbor  bond  populations 
are  shown  in  Table  4.  It  will  be  noticed  that  there  is  very 
little  difference  between  the  two  clusters.  This  is  an  indica¬ 
tion  that  for  diamond,  surface  saturation  with  hydrogen  is  not 
a  particularly  bad  approximation.  There  are  additional  small 
differences ;  reflecting  primarily  the  shorter  hydrogen  carbon 
distance,  the  second  neighbor  C— H  and  H— H  bond  populations  are 
of  the  order  of  30-100  times  greater  than  the  corresponding 
C-C'  and  C'-C'  bond  populations.  The  parameterization  used 
in  this  report  is  summarized  in  Table  3. 


Table  2 

Ethane  Bond  Populations  and  Charge  Distribution 


Bond  Population 

C2H6 

CC'H5 

cc'h3 

C-C 

0.7204 

— 

— 

C-C' 

-- 

0.6877 

0.6496 

C-H 

0.6733 

0.6751 

0.6805 

C'-H 

-- 

0.6674 

— 

Charge 

C 

4.01 

■c* 

o 

o 

4.05 

C' 

— 

3.09 

1.00 

H 

c 

0.997 

1.00 

0.98 

H’ 

c 

-- 

0.961 

— 

Table  3 

CNDO  Parameterization  for 

Element  6  (eV) 

Electronic  Structure 

Y±i  (eV) 

C 

-15.0 

11.1 

B 

-17.0 

— 

0 

-45.0 

— 

Si 

-7.0 

8.0 

P 

-10.0 

— 
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Table  4 


Bond  Populations  and 
Neo-Pentane  and 

Charge 
C5  C12 

Distribution 

Cluster 

Bond  Population 

C5H12 

C5  C12 

C-C 

0.694 

0.686 

C-H 

0.672 

— 

C-C’ 

— 

0.641 

Charge 

C1 

3.910 

3.902 

C2 

4.023 

4.0U7 

H 

1.000 

-- 

C' 

-- 

1.006 

Note:  Ci  is  the  central  atom,  C2 

one  of  four  carbon  atoms  in 
tetrahedral  coordination 
about  Ci* 


Adding  additional  shells  of  carbon  atoms  causes  the  energy 
gap  to  reduce  as  is  illustrated  in  Figure  4.  The  CC^  cluster 
is  a  methane-like  arrangement  of  atoms;  C^C^  consists  of  a 
central  carbon,  its  full  complement  of  four  nearest  neighbor 
and  12  second  neighbor  carbon  atoms,  12  third  neighbor  C' 
atoms  occupying  <3U>  positions  and  having  two  hybrids  oriented 
towards  <220>  type  carbon  atoms,  and  12  C’  atoms  in  fifth 
neighbor  <  33  L>  positions,  each  having  a  single  hybrid.  The 
energy  gap  is  seen  to  decrease  from  13.23  to  8.99  to  7.33  eV 
as  the  cluster  size  increases.  The  experimental  indirect  band 
gap  in  diamond  is  5.33  eV.  Additional  comparison  with  crystal¬ 
line  diamond  can  be  made  with  the  valence  band  width.  Band 
structure  calculations  indicate  that  this  is  from  20-25  eV . 

The  cluster  calculations  show  an  increase  from  15.24  to  28.78 
to  38.19  eV  as  the  cluster  size  increases.  These  last  two 
large  values  aris.e  from  a  few  states  lying  abnormally  low. 

All  but  the  lowest  four  states  in  the  large  cluster,  for 
example,  lie  in  a  range  of  24  eV.  Some  indications  as  to  why 
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this  occurs  can  be  seen  by  considering  the  electron  distribu¬ 
tion  in  these  lowest  orbitals.  The  lowest  orbital  is  an  A 
orbital  (totally  symmetric),  and  in  the  large  cluster,  an 
electron  in  this  orbital  is  distributed  as  follows:  Cq  m  0.14, 
C1  -  0.08,  C2  -  0.03,  C*  -  0.013,  and  C*  -  0.002.  Thus,  this 
orbital  is  localized  toward  the  interior  of  the  cluster,  a 
feature  introduced  by  the  finite  size  of  the  cluster,  since 
this  state  in  the  infinite  solid  would  be  uniformly  distributed 
across  all  atoms.  Since  the  total  charge  distribution  is  what 
would  be  expected  for  this  agglomeration  of  atoms  (Table  5), 
this  means  that  there  are  orbitals  which  are  somewhat  localized 
on  the  surface  —  again  an  artifact  of  the  finite  size  of  the 
cluster . 

The  nearest  neighbor  bond  populations  and  the  atomic 
charge  distribution  of  these  clusters  are  shown  in  Table  5. 
Again,  it  will  be  noted  that  these  are  fairly  consistent  for 
clusters,  although  some  small  variations  are  observed. 


Table  5 

Bond  Populations  and  Charge  Distributions 
in  Diamond  Clusters 


Bond  Population 

cc; 

4 

C5C12 

C17C41 

CCTC1 

— 

0.686 

0.674 

crc2 

-- 

-- 

0.684 

C-C' 

Charge 

* 

0.641 

0.648 

0.635 

(<311>C ' ) 
(<331>C '  ) 

co 

4.010 

3.902 

4.005 

C1 

-- 

4.007 

3.934 

C2 

-- 

-- 

3.985 

C' 

0.998 

1.006 

2.023 

1.014 

(<311>C ' ) 

(<331>C ' ) 

*Not  calculated. 
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b .  Silicon 

The  electronic  structure  of  large  and  small  clusters  of 
silicon  are  illustrated  in  Figure  5.  These  clusters,  Si,.Si{2 
and  Si^.,Si^^  are  analogues  of  the  diamond  clusters  C^Cj^  and 
C17C41‘  The  Parameterization  used  in  these  calculations  is 
shown  in  Table  3;  the  value  of  0  was  adjusted  to  yield  a  reason 
able  energy  gap.  The  energy  gap  Is  seen  to  decrease  from  3.08 
to  2.41  eV  as  the  cluster  size  is  increased,  compared  to  an 
experimental  indirect  band  gap  of  1.17  eV  at  0°K.  The  width 
of  the  valence  states  is  11.42  eV  in  the  snail  cluster  and 
15.12  eV  in  the  large  cluster.  This  compared  favorably  with 
a  value  of1  17  eV  from  band  structure  calculations.  The  nearest 
neighbor  bond  populations  and  charge  distributions,  listed  in 
Table  6,  are  again  consistent  with  a  covalent  structure. 

Table  6 

Bond  Populations  and  Charge  Distribution 
of  Silicon  Clusters 


Bond  Population 

Si5Sii2 

Si17Si24 

Sio-Sii 

0.752 

0.732 

Sil'Si2 

-- 

0.749 

Si-Si' 

0.710 

0.718  (<311>Si') 

0.705  (<331>Si  '  ) 

Charge 

Sio 

3.939 

4.102 

S11 

3.983 

4.010 

Si2 

-- 

3.958 

Si' 

1.011 

2.027  (<311>Si') 

1.004  (<331>Si ' ) 
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c .  B-Silicon  Carbide 


An  acid  test  of  the  semiempirical  approach,  and  especially 
the  use  of  surface  hybrids  is  in  6-silicon  carbide  (sphalerite 
structure),  a  solid  which  contains  two  types  of  atoms.  Two 
small  clusters  were  used,  one  containing  a  central  carbon  atom 
surrounded  by  four  silicon  atoms  and  twelve  C*  surface  hybrids; 
the  other  containing  a  central  silicon  atom  surrounded  by  four 
carbon  atoms  and  twelve  Si'  surface  hybrids.  The  electronic 
structure  of  these  two  clusters  is  illustrated  in  Figure  6; 
the  energy  gap  of  the  cluster  with  carbon  at  its  center  is 
3.63  eV,  that  of  the  cluster  with  Si  at  its  center,  3.30  eV, 
while  the  experimental  indirect  band  gap  is  2.68  eV  at  0°K. 

Even  more  remarkable  is  the  agreement  in  bond  populations 
and  atom  charges  (Table  7),  indicating  that  the  surface  hybrids 
actually  appear  to  the  interior  of  the  cluster  as  the  atoms 
they  represent. 


Table  7 

Bond  Populations  and  Charge  Distribution 
of  6-SiC  Clusters 


C  at 

Origin 

Si 

at  Origin 

Bond 

Population 

Bond 

Population 

Si-C 

0.663 

Si-C 

0.66  7 

Si-C* 

0.618 

C-Si  ’ 

0.612 

Charge 

Charge 

C 

4.508 

C 

4.416 

Si 

3.421 

Si 

3.210 

C' 

1.151 

Si* 

0.927 

d .  SIO^  Clusters 

The  electronic  structures  of  the  silica  clusters  SiO.Sij! 

4  5 

and  Si^O^Si^  are  shown  in  Figure  7.  The  small  cluster  con¬ 
sists  of  a  central  silicon  atom  surrounded  by  four  oxygen  atcms 
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in  tetrahedral  coordination  and  four  silicon  surface  hybrids 
such  that  the  Si-O-Si'  bond  angle  is  180°.  The  large  cluster 
contains  a  central  silicon  atom,  four  oxygen  atoms  in  tetra¬ 
hedral  arrangement,  four  silicon  atoms,  twelve  oxygen  atoms— 
three  each  attached  to  the  four  silicon  atoms  in  tetrahedral 
coordination--  and  twelve  silicon  surface  hybrids:  all  Si-O-Si 
bond  angles  are  180°.  These  configurations  correspond  to  the 
high  crlstoballlte  configuration  of  silica.  The  energy  gap  of 
the  small  cluster  is  10.02  eV,  while  that  of  the  large  cluster 
is  10.25  eV .  This  compares  favorably  with  the  experimental 
value  of  10  eV .  The  nearest-neighbor  bond  populations  and 
atomic  charge  distributions  are  shown  in  Table  8.  Again,  the 
results  are  good.  The  charge  distribution  indicates  that  the 
structure  is  about  63%  covalent — not  far  from  the  value  of  50% 
predicted  by  Pauling's  formula. 


Table  8 

Bond  Populations  and  Charge  Distributions 
for  Silica  (High  Crlstoballlte)  Clusters 


Bond 

SiO, Si! 

- 4 — 4 

—5—12—12 

slo-°i 

0.596 

0.591 

S1l-°1 

— 

0.590 

S1l-°2 

— 

0.596 

Si'-O 

0.545 

0.545 

Charge 

S10 

2.489 

2.610 

S11 

-- 

2.515 

°1 

6.678 

6.707 

°2 

— 

6.674 

Si' 

0.700 

0.702 
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e .  Summary 

The  preponderence  of  evidence  presented  in  the  foregoing 
sections  is  that  the  CNDO  method  in  conjunction  with  surface 
hybrids  does  reasonably  describe  the  electronic  structure  of 
a  variety  of  solids.  Furthermore,  deviations  are  predictable 
with  the  given  parameterization;  for  example,  the  energy  gaps 
are  consistently  too  large.  Charge  distributions  undergo  a 
small  fluctuation  from  atom  to  atom;  however,  results  are  con¬ 
sistent  to  within  about  3%.  It  thus  appears  that  the  method 
is  ready  to  be  applied  to  defect  clusters. 


B .  Impurity  and  Vacancy  Complexes 

In  this  section  we  shall  consider  the  electronic  structure 
of  clusters  containing  impurity  atoms  and  vacancies.  Of  par¬ 
ticular  interest  in  understanding  the  effects  of  radiation 
damage  in  crystalline  semiconductors  is  impurity  clustering 
and  impurity-vacancy  complexing.  It  is  of  interest  to  deter¬ 
mine  the  structure  of  these  defects  in  order  to  better  under¬ 
stand  the  resulting  properties  of  radiation-damaged  devices. 

We  shall  therefore  first  consider  single  boron  and  phosphorus 
impurities  in  silicon.  While  these  are  shallow  acceptors  and 
donors,  respectively,  and  can  therefore  be  treated  by  perturba¬ 
tion  methods,  a  preliminary  study  is  necessary  by  CNDO  tech¬ 
niques. 

The  electronic  structure  of  the  small  silicon  cluster  with 
a  boron  atom  replacing  the  central  silicon  atom  is  shown  in 
Figure  7.  Boron  is  seen  to  give  rise  to  a  shallow  acceptor 
state  0.21  eV  above  the  valence  state;  the  energy  gap  of  the 
cluster  is  4.46  eV  (the  energy  gap  of  the  pure  silicon  cluster 
is  3.08  eV,  and  the  experimental  band  gap  1.17  eV).  If  the 
energies  scale  linearly,  this  would  correspond  to  an  acceptor 
state  0.055  eV  above  the  valence  band  in  the  crystal--the 
experimental  value  being  0.05  eV. 
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The  electron  distribution  Is  shown  In  Figure  8.  The 
charge  deficiency  (or  excess)  is  defined  as  the  difference 
between  the  charge  on  the  defect  cluster  and  the  charge  on 
the  pure  silicon  cluster.  Because  of  the  CNDO  approximations, 
there  is  an  error  of  the  order  of  ±0.04  electrons  per  atom 
(or  bond)  relative  to  each  other  in  these  calculations.  This 
discrepancy  arises  from  termination  of  the  self-consistency 
iterative  procedure.  Since  the  convergence  is  slower  in  the 
open  shell  calculation  than  in  the  closed  shell  computation, 
termination  after  a  specified  increment  in  energy  results  in 
the  open  shell  calculation  being  further  from  self-consistency 
than  the  closed  shell  calculation.  A  second  source  of  error 
is  the  normalization  of  the  molecular  orbitals.  Since  the 
orbitals  are  normalized  assuming  zero  overlap,  the  charges 
derived  have  only  a  semiquantitative  meaning.  For  this  reason, 
all  charge  deficiencies  (or  excesses)  should  be  considered  as 
qualitative  results. 

The  excess  spin  distribution  is  sensitive  to  the  self- 
consistency  of  the  cluster,  and  therefore  differences  of  ±0.04 
electrons  should  be  considered  insignificant.  One  puzzling 
result  is  the  rather  large  deficiency  on  one  surface  hybrid 
(Figure  8c)  of  -0.296  and  a  surplus  of  0.266  on  a  neighboring 
silicon.  This  will  probably  not  result  in  any  great  error, 
since  this  results  in  only  about  1/4  or  0.025  eV  change  in 

the  interaction  energy.  This  unusual  situation  was  noted  in 
several  of  the  open  shell  calculations.  Note  that  this  is  a 
shift  in  spin  density  only;  the  overall  charge  distribution 
(Figure  8a)  remains  uniform. 

The  results  of  open  shell  calculations  for  the  large 
cluster  with  the  boron  atom  in  both  a  central  and  noncentral 
(111)  position  are  shown  in  Figures  9-12.  The  position  of  the 
acceptor  level  is  seen  to  be  relatively  insensitive  to  atom 
position,  with  the  boron  in,  the  central  position,  the  acceptor 
state  is  0.40  eV  above  the  valence  states,  with  an  energy  gap 
of  3.20  eV  (Figure  9),  while  with  the  boron  off-center,  the 
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acceptor  state  is  0.54  eV  above  the  valence  states,  with  an 
energy  gap  of  2.45  eV .  The  total  energies  of  these  two  clus¬ 
ters  differ  by  1.5  eV. 

Concerning  the  electron  distribution  in  the  clusters  the 
following  generalizations  can  be  made:  (1)  The  atoms  remain 
fairly  neutral ;  that  is,  the  boron  "soaks  up"  few  electrons 
from  the  silicon  (Figures  8a,  10a,  12a) .  Since  the  acceptor 
state  is  not  strongly  peaked  on  the  boron  atom  (Figures  8d, 

10c,  12c),  this  means  that  lower-lying  occupied  orbitals  have 
smaller  occupancies  on  the  boron  than  on  the  silicon  atoms. 

(2)  The  acceptor^ state  (Figures  8d,  10c,  12c)  is  fairly  spread 
out.  The  peaking  which  does  occur  appears  to  be  more  a  result 
of  the  finite  size  of  the  cluster  than  of  the  impurity  itself 
(compare  Figures  10c  and  12c,  in  particular).  (3)  The  Si-B 
bond  is  considerably  weaker  than  the  Si-Si  bond  (Figure  8b). 

The  charge  deficiency  in  this  bond  is  not  completely  unexpected, 
since  the  boron  atom  is  electron  deficient.  (4)  The  excess  spin 
distribution  roughly  corresponds  to  the  charge  distribution  in 
the  acceptor  state  (Figures  8c, d,  10b, c,  12b, c)  This  is  not 
totally  unexpected,  for  if  we  can  speak  of  a  hole  in  silicon, 
then  it  must  have  spin.  The  spin  must  be  opposite  in  sign  to 
and  have  the  distribution  of  the  electron  which  is  missing. 

(5)  There  is  an  indication  that  the  cluster  would  like  to  under¬ 
go  a  Jahn-Teller  distortion.  The  clusters  with  the  atom  at  its 
center  have  Td  symmetry,  yet  the  electron  distributions  deviate 
from  this.  No  attempt  was  made  to  introduce  these  distortions 
into  the  calculation.  (6)  Considering  the  magnitude  of  the 
energy  gaps  in  comparison  with  those  of  the  pure  silicon  clus¬ 
ter,  it  appears  that  there  should  be  a  tendency  for  the  band 
gap  to  widen  in  the  vicinity  of  a  boron  impurity  in  the  crystal. 

The  results  on  clusters  containing  a  single  phosphorus 
impurity  are  shown  in  Figures  13-16.  During  the  self-consistency 
iteration  on  the  large  cluster  with  the  phosphorus  at  the  center 
of  the  cluster,  the  energy  diverged  (i.e.,  initially  decreased, 
then  increased  by  progressively  larger  amounts).  The  reason 
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for  this  divergence  is  not  understood;  consequently,  the  results 
on  this  cluster  are  not  presented. 

The  phosphorus  atom  is  seen  to  give  rise  to  a  donor 
state  lying  0.39  eV  below  the  conduction  states,  with  an  energy 
gap  of  3.51  eV  in  the  small  cluster  (Figure  13),  and  a  donor 
state  lying  0.50  eV  below  the  conduction  states  with  an  energy 
gap  of  2.20  eV  in  the  large  cluster  (Figure  15).  It  should  be 
remembered  that  the  energy  gap  is  defined  as  an  excitation  from 
the  top  of  the  valence  states  into  the  conduction  states,  and 
consequently  the  J  and  K  integrals  in  equation  (21)  are  com¬ 
puted  between  the  impurity  state  and  conduction  states  for  the 
depth  of  the  donor  level,  and  the  top  of  the  valence  states  and 
the  conduction  state  to  determine  the  energy  gap.  This  some¬ 
what  roundabout  method  is  partially  a  result  of  considering  a 
finite  cluster,  and  partially  a  result  of  the  use  of  a  Hartree— 
Fock  procedure,  as  discussed  in  the  first  section. 

The  discussion  for  the  boron  impurity  is  also  applicable 
to  the  phosphorus  impurity.  The  only  surprising  result  is 
that  the  bonding  between  phosphorus  and  silicon  is  weaker  than 
between  silicon  atoms  (Figure  14b).  On  an  intuitive  basis,  the 
opposite  would  have  been  expected,  since  the  phosphorus  atom 
has  an  excess  of  electrons.  In  fact,  the  phosphorus  -silicon 
bonds  are  charge  deficient,  as  are  the  neighboring  silicon- 
silicon  bonds.  This  is  probably  the  explanation  for  the  limited 
solid-solubility  of  P  in  Si. 

It  is  interesting  to  compare  the  results  of  the  hydrogenic 
perturbation  treatment  of  an  impurity  state  with  the  CNDO 
treatment.  In  the  hydrogenic  model,  a  point  charge  is  placed 
on  a  lattice  site.  The  surrounding  medium  is  polarized  by  the 
charge ,  reducing  the  effective  charge  of  the  impurity  center 
from  unity  to  1/K  where  K  is  the  dielectric  constant.  The 
effective  charge  then  acts  as  the  nucleus  of  a  hydrogen  atom 
with  an  electron  having  effective  mass  m*  in  orbit  about  the 
center.  The  hydrogenic  model  makes  no  attempt  to  explain  what 
is  happening  close  to  the  impurity  center;  in  fact,  the  pertur¬ 
bation  technique  breaks  down  close  in  to  the  center.  It  will 
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be  a  good  description  provided  that  (1)  the  close-in  effects 
are  short  range,  and  (2)  that  the  impurity  state  is  spread 
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tion  on  an  impurity  state  far  from  the  defect,  we  must  have 
(1)  a  large  cluster,  preferably  with  periodic  boundary  condi¬ 
tions,  and  (2)  a  high  degree  of  self-consistency,  since  elec¬ 
tron  differences  are  small.  On  the  other  hand,  we  can  obtain 
accurate  information  close  in  to  the  defect  without  either  of 
the  above.  Confining  the  carrier  to  a  small  volume  makes  it 
become  strongly  influenced  by  the  structure  close  to  the 
center,  and  gives  rise  to  something  closer  to  an  impurity 
bound  excitonic  state.  The  impurity  states  should  therefore 
lie  somewhat  deeper  than  in  reality,  as  is  in  fact  found. 

In  the  hydrogenic  model,  the  charge  of  the  impurity  center 
is  reduced  by  the  appearance  of  a  virtual  charge  caused  by 
polarization  of  the  medium.  In  the  CNDO  model,  the  charge  of 
the  impurity  center  is  reduced  by  a  real  charge  shift  caused 
by  alteration  of  the  populated  molecular  orbitals. 

The  CNDO  model  should  therefore  be  best  when  considering 
localized  defects.  A  particular  example  of  this  is  in  the 
chemical  bonding,  between  the  impurity  atom  and  the  lattice. 

The  detailed  structure  of  the  impurity  state  far  from  the  cen¬ 
ter  will  be  largely  masked  by  the  finite  size  of  the  cluster, 
and  the  limited  degree  of  self-consistency  in  the  calculation. 
However,  semiquantitative  results  can  still  be  obtained,  as 
has  been  demonstrated. 

1 .  Vacancy  in  Silicon 

The  problem  of  a  vacancy  is  one  which  is  beyond  the  scope 
of  perturbational  treatments.  The  vacancy  in  diamond  has  been 
studied  by  Swenson  using  the  Extended  Huckel  Theory  (a  non- 
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S6lf (13) SlStent  fleld  method)  and  surface  saturation  with  hydro 
gen.  In  the  present  work,  we  will  use  a  similar  procedure 

with  the  CNDO  method. 

A  vacancy  was  created  in  the  large  silicon  cluster  by 
removing  the  central  atom.  The  remaining  atoms  maintained  the 
Td  symmetry  of  the  initial  cluster.  The  atoms  surrounding  the 
vacancy  are  shown  in  Figure  17a.  When  the  central  atom  is 
removed,  these  four  atoms  will  each  have  one  unsaturated  sp3 
hybrid  oriented  toward  the  center.  This  configuration  is  not 
stable;  we  expect  the  cluster  to  distort  so  that  covalent  bonds 
between  the  atoms  can  form.  Experimentally,  this  distortion  is 
known  to  have  D2<1  symmetry .  (25  >  A  coordinate  system  describing 
all  possible  distortions  from  Td  symmetry  to  D2d  symmetry  may 
be  constructed  as  in  Figure  17b.  The  positions  of  the  other 
three  central  atoms  will  be  related  to  the  one  shown  through 
the  symmetry  operations  of  the  group  D2d.  Only  the  four  cen¬ 
tral  atoms  were  moved. 

To  determine  the  equilibrium  (lowest  energy)  position  of 
these  atoms  a  potential  surface  is  constructed  showing  the 
energy  of  the  cluster  as  a  function  of  R  and  0.  Such  a  project 
is  rather  formidable  even  in  the  semiempirical  CNDO.  (Typically 
10  minutes  360/65  computer  time  or  15-20  minutes  370-155  time 
with  a  368K  partition  are  required  for  each  point.)  Therefore, 
limited  sections  through  the  potential  surface  were  made,  based 
on  the  assumption  that  the  vacancy  is  qualitatively  similar  to 
that  in  diamond.  These  section",  are  shown  in  Figure  18.  An 
additional  section  was  attempted  with  0  =  145°;  however,  the 
initial  Huckel  approximation  mispredicted  the  ordering  of  the 

levels,  and  self-consistency  still  had  not  been  achieved  after 
nine  iterations. 

In  accordance  with  general  practice,  the  CNDO/2  parameter¬ 
ization  of  Pople ,  Santry  and  Segal(1_4)  was  used  to  determine 
the  potential  surface,  since  this  parameterization  is  known  to 
give  good  equilibrium  bond  lengths.  Some  difficulty  was 

encountered  with  the  ooinf-Q  fnr  n  -  n  c  .  , 

v  iuis  tor  k  -  u.o,  the  electronic  energy 
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initially  decreased,  then  increased.  The  energy  is  therefore 
not  at  self-consistency;  these  points  may  consequently  lie 
somewhat  higher  or  lower.  The  minimum  energy  in  Figure  18  is 
found  to  occur  at  approximately  R  =  0.35,  9  *  110°.  Because  of 
the  limited  nature  of  the  potential  surface  constructed,  this 
value  should  be  considered  as  approximate. 

It  is  interesting  to  note  that  the  electronic  structure 
in  the  vicinity  of  the  energy  gap  was  relatively  insensitive 
to  the  tetrahedral  distortion.  For  example,  using  the  Pople, 
Santry  and  Segal  parameterization,  the  energy  from  the  highest 
populated  state  to  the  lowest  unpopulated  state  is  0.23  eV, 
and  the  energy  from  the  highest  populated  state  to  the  conduc¬ 
tion  states  is  2.8  eV  for  the  cluster  with  R  =  0.08,  0  =  70°. 

For  the  cluster  with  R  =  0.35,  9  =  110°,  these  values  are 
0.42  eV  and  3.06  eV,  respectively. 

The  energy  structure  of  the  cluster  using  the  parameteri¬ 
zation  of  Table  3  is  illustrated  in  Figure  19.  It  will  be 
noticed  that  a  vacancy  acts  as  an  acceptor  state.  Comparing 
this  structure  with  that  from  the  pure  Si  cluster  (Figure  5), 
it  will  be  seen  that  the  vacancy  splits  off  an  occupied  level 
from  the  valence  states  and  introduces  tvo  acceptor  states 
slightly  above  this  level.  The  details  of  how  this  structure 
fits  into  the  overall  band  structure  cannot  be  determined;  the 
calculation  indicates  the  accepter  states  lie  0.27  eV  above 
the  populated  states.  Additionally,  the  overall  structure 
appears  to  indicate  that  the  occupied  state  split  off  from  the 
valence  states  lies  fairly  deeply  in  the  band  gap.  If  simple 
energy  scaling  is  valid,  the  structure  would  appear  as  in 
Figure  20.  Further  calculations  with  impurities  added  to  the 
vacancy  cluster,  which  will  be  discussed  later,  indicate  that 
this  is  a  not  unreasonable  'model. 

The  electron  distribution  for  the  vacancy  cluster  is  shown 
in  Figure  21.  There  is  seen  to  be  some  depletion  of  charge 
from  the  atoms  surrounding  the  vacancy  (Figure  21a).  The  charge 
piles  up  on  surface  hybrids  (not  shown).  This  may  be  a  result 
of  the  surface  boundary  conditions. 
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The  bond  populations  (Figure  21b)  show  the  formation  of 
a  bond  between  pairs  of  silicon  atoms  surrounding  the  vacancy. 
With  the  formation  of  these  bonds,  there  is  a  slight  weakening 
of  bonds  with  neighboring  silicon  atoms.  The  equilibrium 
positions  of  the  atoms  are  those  in  which  the  various  bonds 
are  maximized. 

The  charge  distribution  in  various  orbitals  in  the  vicinity 
of  the  energy  gap  are  illustrated  in  Figures  21c-f.  Figure  21c 
shows  that  the  highest  populated  orbital  is  fairly  peaked  cn 
the  four  atoms  surrounding  the  vacancy  (60%  on  these  atoms); 
we  are  thus  justified  in  speaking  of  this  orbital  as  arising 
from  the  vacancy.  Figure  21d  shows  one  of  the  pair  of  doubly 
degenerate  acceptor  states  lying  in  the  energy  gap.  This  state 
is  also  rather  peaked  on  a  pair  of  central  silicon  atoms  (50% 
on  a  pair  —  the  other  molecular  orbital  of  this  set  is  peaked 
on  the  other  pair  of  silicon  atoms),  and  may  also  be  definitely 
attributed  to  the  vacancy.  On  the  other  hand,  states  higher 
or  lower  in  energy  than  the  above  states  are  spread  out  through¬ 
out  the  cluster  (Figures  21e,f)  and  may  therefore  be  spoken  of 
as  conduction  and  valence  states. 

The  rather  large  degree  of  localization  of  the  states 
within  the  energy  gap  is  a  fairly  good  indication  that  they 
lie  rather  deep  within  the  energy  gap.  This  observation  lends 
support  to  the  energy  diagram  sketched  in  Figure  20. 

2 .  Vacancy-Impurity  Complexes 

laving  considered  the  pure  vacancy,  the  next  step  is  to 
add  impurities  to  the  silicon  cluster  to  see  how  impurities 
and  vacancies  interact. 

The  first  problem  we  shall  consider  is  a  single  bcron 
impurity  added  to  a  vacancy  cluster.  Two  such  problems  were 
considered;  in  the  first  the  boron  atom  was  placed  on  one  of 
the  four  positions  surrounding  the  vacancy  (nearest  neighbor 
model);  in  the  second,  the  boron  atom  was  placed  on  one  of  the 
neighboring  silicon  atoms  (next  nearest  neighbor  model). 
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additional  distortions  were  considered,  although  considering 
the  change  of  charge  in  the  vicinity  of  the  vacancy,  a  small 
additional  distortion  does  seem  likely. 

The  energy  structure  of  the  nearest  neighbor  boron  model 
is  shown  in  Figure  22.  It  will  be  noted  that  the  boron  atom 
behaves  as  we  intuitively  might  have  expected.  Since  it  is  a 
shallow  acceptor  in  silicon,  and  since  the  vacancy  has  split 
off  a  pair  of  electrons  from  the  valence  band  fairly  deeply 
into  the  energy  gap,  then  it  would  appear  that  on  adding  a 
boron  atom,  the  boron  atom  will  be  ionized,  accepting  an  elec¬ 
tron  from  the  vacancy  complex.  This  is  in  fact  observed. 

The  boron  has  an  additional  effect,  and  that  is  to  cause 
the  doubly  degenerate  pair  of  acceptor  levels  arising  from  the 
vacancy  to  split  apart.  This  "chemical  splitting,"  as  we 
shall  see,  is  dependent  upon  the  separation  of  the  vacancy  and 
the  boron  atom.  Again  assuming  that  we  can  scale  the  energies, 
the  resulting  energy  level  diagram  is  shown  in  Figure  23. 

The  charge  distribution  for  this  cluster  is  shown  in 
Figure  24.  The  cluster  is  seen  to  be  close  to  being  neutral 
(Figure  24a),  as  in  the  impurity  clusters.  The  bond  popula¬ 
tions,  however,  have  been  altered  (Figure  24b).  In  particular, 
the  bond  which  formed  between  atoms  2  and  4  during  distortion 
of  the  vacancy  is  seriously  weakened  (bond  population  from 
0.30  to  0.12)  while  that  between  atoms  1  and  3  is  weakened 
slightly.  The  weakening  of  these  bonds  is  most  likely  due  to 
the  charge  deficiency  on  the  boron  atom,  which  allows  one  less 
electron  to  take  part  in  the  bonding  process. 

The  distribution  of  the  excess  spin  (Figure  24c)  is  rather 
interesting  in  that  it  lends  support  to  our  previous  argument 
that  the  ionized  boron  acceptor  is  spread  out,  while  the 
occupied  vacancy  levels  are  localized.  In  fact,  the  molecular 
orbital  of  the  highest  populated  state  (the  unionized  vacancy 
state)  is  found  to  have  a  distribution  close  to  that  of  the 
exces?  spin. 
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Moving  the  boron  to  the  next  nearest  neighbor  position 
results  in  the  energy  structure  shown  in  Figure  25.  Again 
the  boron  acceptor  is  ionized  by  an  electron  from  one  of  the 
vacancy  states.  The  chemical  splitting  of  the  vacancy  acceptor 
states  is  seen  to  be  substantially  less.  A  scaled  energy  level 
diagram  is  shown  in  Figure  26.  We  can  anticipate  that  as  the 
boron  atom  is  moved  further  from  the  vacancy,  the  chemical 
splitting  will  become  less  and  less,  even  though  the  boron 
acceptor  state  remains  ionized. 

The  charge  distribution  with  the  boron  atom  in  the  next 
neighbor  position  is  illustrated  in  Figure  27a.  Charge  neu¬ 
trality  is  still  nearly  maintained.  The  weakening  of  the 
vacancy  bonds  is  still  observed  (Figure  27b),  although  not  to 
as  marked  degree  as  when  the  boron  was  a  nearest  neighbor  atom. 

The  excess  spin  distribution  (Figure  27c)  is  similar  to 
that  which  was  found  in  the  nearest  neighbor  boron  cluster, 
and  is  indicative  of  a  localized  orbital  in  the  energy  gap. 

It  is  interesting  to  note  that  the  atom  with  most  excess  spin 
remains  the  same;  that  is,  the  silicon  atom  which  is  paired 

with  the  atom  closest  to  the  impurity  atom  through  the  vacancy 
bond.. 

The  properties  of  a  phosphorus  -vancacy  complex  have  also 
been  determined.  An  attempt  to  calculate  the  structure  of  such  a 
complex  when  the  phosphorus  atom  was  one  of  the  four  atoms  sur¬ 
rounding  the  vacancy  (nearest  nieghbor  model)  was  unsuccessful; 
the  energy  was  found  to  diverge  during  the  self-consistency 
procedure.  This  divergence  resulted  from  a  large  amount  of 
charge  moving  on  to  the  atoms  surrounding  the  vacan  y. 

A  calculation  of  a  cluster  in  which  the  phosphorus  atom 
was  a  next  nearest  neighbor  was  successful,  however.  The  energy 
structure  of  this  complex  is  shown  in  Figure  28.  This  combina¬ 
tion  too  behaves  as  one  would  intuitively  expect.  The  phos¬ 
phorus  atom,  being  a  shallow  donor,  becomes  ionized,  yielding 
its  electron  to  one  of  the  vacancy  acceptor  states.  A  scaled 
energy  level  diagram  is  shown  in  Figure  29. 
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The  charge  distribution  (Figure  30a)  again  shows  that  the 
cluster  approximately  maintains  neutrality.  The  bonding  is 
similar  to  that  which  was  observed  with  boron;  that  is,  a 
weakening  of  the  vacancy  bonds.  Interestingly,  the  populations 
of  the  vacancy  bonds  are  smaller  in  the  phosphorus  cluster 
than  they  were  in  the  boron  cluster;  this  even  though  the  phos¬ 
phorus  atom  has  an  excess  of  electrons. 

The  excess  spin  distribution  (Figure  30c)  is  most  inter¬ 
esting,  in  that  it  shows  signs  of  spirt  ordering.  The  reason 
for  this  behavior  can  be  seen  by  considering  the  three  populated 
orbitals  with  levels  in  the  energy  gap.  Two  of  these  have  a 
spin;  one  of  these  is  peaked  sharply  on  atom  3  and  more  weakly 
on  atom  1;  the  other  is  peaked  sharply  on  atom  4  and  more 
weakly  on  atom  1.  The  other  orbital  has  B  spin;  it  is  peaked 
sharply  on  atoms  1  and  2.  Consequently,  the  spins  on  atom  1 
cancel,  leaving  a  net  8  spin  on  atom  2  and  a  spins  on  atoms 
3  and  4,  as  the  excess  spin  populations  show.  That  this  occurs 
can  be  attributed  to  a  high  degree  of  localization  of  the 
states.  For  several  electrons  confined  to  a  few  atoms,  the 
Coulomb  integral  is  rather  large.  The  magnitude  of  this 

integral  may  be  reduced  if  the  states  of  the  system,  rather 
than  being  spread  out  over  all  atoms,  are  confined  to  just  a 
few.  This  phenomenon  is  known  variously  as  a  Mott  transition  ^2^ 
or  Wigner  crystallization.^27^ 


Multi-Impurity  Clusters 


Let  us  now  turn  our  attention  to  clusters  in  which  there 
are  two  or  more  impurities  present.  These  clusters  will  be 
representative  of  regions  in  the  crystalline  semiconductor 
where  impurity  clustering  occurs  due  to  heavy  doping,  diffusion, 
or  radiation  damage . 

The  energy  structure  of  a  cluster  containing  two  boron 
atoms  which  are  nearest  neighbors  is  shown  in  Figure  31.  A 
cluster  in  which  two  boron  atoms  are  separated  by  a  second 
neighbor  distance  is  presented  in  Figure  32.  In  both  of  these 
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clusters,  a  shallow  acceptor  state  capable  of  holding  two  elec¬ 
trons  is  formed.  The  charge  distribution  in  this  acceptor 
state  is  depicted  in  Figure  33c  for  the  nearest  neighbor  borons, 
and  in  Figure  34c  for  the  next  nearest  neighbor  boron  cluster. 

In  both  these  cases  there  appears  to  be  a  tendency  for  the 
electron  in  the  acceptor  state  to  avoid  the  boron  atom  itself, 
as  was  the  case  for  a  single  boron  atom  (Figures  10  and  12). 

Two  important  questions  are:  (1)  which  of  the  two  arrange¬ 
ments  of  boron  atoms  is  the  most  stable,  and  (2)  is  a  pairing 
of  the  boron  atoms  stable  with  respect  to  completely  isolated 
atoms  in  the  silicon  crystal?  Replacing  a  silicon  atom  with 
a  boron  in  the  central  position  results  in  an  increase  in  the 
total  cluster  energy  of  0.9407  Hartrees  ,  while  replacing  it  in 
the  off-center  position  results  in  an  increase  of  0.8865  Har¬ 
trees.  Thus,  if  there  were  no  interactions  between  boron 
atoms,  the  energy  of  the  nearest  neighbor  complex  should  be 
1.8272  Hartrees  above  the  pure  silicon  cluster,  while  that  of 
the  second  neighbor  arrangement  should  be  1.7729  Hartrees  above. 
In  actuality,  the  nearest  neighbor  arrangement  has  an  energy 
1.6525  Hartrees  above  the  pure  silicon  energy,  while  the 
second  neighbor  arrangement  of  boron  atoms  has  an  energy  1.6706 
Hartrees  above.  This  means  that  the  interaction  energy  between 
two  boron  atoms  in  the  nearest  neighbor  configuration  is 
-0.1747  Hartrees,  while  that  of  the  next  nearest  neighbor 
arrangement  is  —0.1023  Hartrees.  Thus,  the  nearest  neighbor 
arrangement  of  boron  atoms  is  about  4.75  eV  more  stable  than 
having  them  separated  an  infinite  distance,  and  about  2.0  eV 
more  stable  than  the  next  nearest  neighbor  arrangement.  It 
appears  that  pair  clustering  of  boron  impurities  in  silicon 
is  favored  over  a  random  distribution,  at  least  in  the  present 
model.  It  should  be  noted  that  here  we  are  treating  only  a 
chemical  effect,  since  no  elastic  distortions  have  been  added. 
Distortions  may  alter  this  chemical  effect  and,  in  addition, 
will  introduce  an  elastic  interaction  energy  which  will  probably 
be  repulsive  at  all  distances  (the  boron  atom  is  smaller  than 


109 


the  silicon,  hence  it  should  introduce  a  dilatational  center 
of  strain) . 

The  energy  structure  of  a  silicon  cluster  containing  two 
phosphorus  atoms  separated  by  a  second  neighbor  distance  is 
shown  in  Figure  35.  There  is  a  rather  deep  donor  state  (doubly 
occupied).  The  distribution  of  charge  in  both  the  highest 
occupied  orbital  and  the  lowest  unoccupied  orbital  is  somewhat 
localized  (Figure  36c, d),  but  a  good  portion  of  this  localiza¬ 
tion  is  probably  due  to  cluster  boundary  conditions  (compare 
with  Figure  14).  This  seems  to  argue  for  the  two  phosphorus 
atoms  compensating  each  other;  that  is,  the  occupied  donor 
levels  lie  not  far  from  the  top  of  the  valence  band,  while  the 
unoccupied  levels  lie  not  far  from  the  bottom  of  the  conduction 
band.  This  is  in  contrast  to  the  case  of  the  two  boron  impuri¬ 
ties  where  there  was  little  interaction  between  the  two  impurity 
atoms . 

The  two  phosphorus  atom  complex  is  found  to  be  more 
stable  than  the  isolated  atoms  by  approximately  2.2  eV.  An 
attempt  was  made  to  perform  a  calculation  on  a  complex  with  the 
phosphorus  atoms  separated  by  a  nearest  neighbor  distance,  but 
the  energy  was  found  to  diverge.  Thus,  it  is  not  known  if  the 
nearest  neighbor  arrangement  would  be  more  stable  than  the  next 
nearest  arrangement. 

The  electronic  structures  of  silicon  clusters  having  both 
single  boron  and  single  phosphorus  impurity  atoms  separated  by 
a  nearest  neighbor  distance  and  a  next  nearest  neighbor  dis¬ 
tance  are  shown  in  Figures  39  and  40.  In  both  clusters  compen¬ 
sation  is  evident,  the  phosphorus  impurity  state  having  been 
ionized  into  the  boron  acceptor  level.  Comparing  these  figures 
with  the  energy  structures  of  the  two  single  boron  impurity 
clusters  in  the  center  and  off-center  position  (Figures  9  and 
H)»  it  will  be  noticed  that  the  differences  are  mainly  due  to 
the  cluster  surface  conditions. 

The  relative  stability  of  these  two  clusters  can  be  deter¬ 
mined  as  follows.  On  replacing  an  off-center  silicon  atom  with 
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a  phosphorus  atom  the  total  energy  of  the  cluster  decreases 
by  2.5125  Hartrees.  Using  the  values  previously  given  for  the 
boron  atoms,  with  no  interaction  energy  between  the  phosphorus 
and  boron,  the  nearest  neighbor  cluster  total  energy  will  be 
1.5718  Hartrees  below  that  of  a  pure  silicon  cluster,  while 
that  of  the  next  nearest  neighbor  cluster  will  be  1.6261 
Hartrees  lower.  In  fact,  these  clusters  are  1.7204  and  1.7904 
Hartrees  lower.  This  means  that  the  second  neighbor  cluster 
is  about  4.5  eV  more  stable  than  isolated  impurities,  and  about 
0.5  eV  more  stable  than  the  nearest  neighbor  cluster.  Hence, 
the  chemical  interaction  between  boron  and  phosphorus  impuri¬ 
ties  is  attractive  at  large  distances  and  repulsive  at  close 
distances . 

The  question  as  to  why  the  boron-phosphorus  interaction 
is  repulsive  at  close  distances,  while  that  of  the  boron— boron 
is  attractive,  may  be  seen  in  the  following.  In  both  types  of 
clusters,  the  covalent  portion  of  the  bonding  decreases  as  the 
impurities  are  separated  (the  sum  of  the  bond  orders  decreases 
by  about  0.1,  see  Figures  33b,  34b,  39b  and  40b).  The  ionic 
portion  of  the  bonding  of  both  types  of  clusters  improves  as 
the  impurities  are  separated.  A  simple  calculation  of  the 
nearest  neighbor  contribution  to  the  ionic  bonding  (in  units 
of  products  of  the  charges)  shows  that  in  the  two  boron  atom 
cluster  the  ionic  bonding  increases  from  0.16  to  0.26  while 
that  of  the  boron-phosphorus  cluster  increases  from  0.02  to 
0.18.  Thus,  in  the  boron-phosphorus  cluster,  the  increase  in 
ionic  bonding  offsets  the  decrease  in  covalent  bonding  as  the 
atoms  are  separated.  The  rather  low  value  of  ionic  energy  for 
the  nearest  neighbor  BP  cluster,  and  the  reason  for  the 
sharper  increase  of  the  ionic  energy  in  the  BP  clusters  as 
the  impurities  are  separated  is  a  chemical  interaction  of 
these  atoms  at  close  distances. 

It  is  also  interesting  that  boron-phosphorus  complexing 
is  predicted  to  be  more  stable  when  both  species  are  present, 
at  least  considering  the  present  models.  Two  boron-phosphorus 
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complexes  would  be  about  2  eV  more  stable  than  one  phosphorus- 
phosphorus  and  one  boron-boron  complex.  It  should  be  empha¬ 
sized,  however,  that  the  only  consideration  made  at  arriving 
at  these  conclusions  is  a  chemical  interaction.  Consideration 
of  distortions  and  elastic  interactions  may  alter  these  results 
A  final  impurity  cluster  which  was  considered  was  one 
containing  two  phosphorus  impurities  and  one  boron.  The 
energy  structure  of  this  cluster  is  shown  in  Figure  41,  and  a 
Scaled  Energy  Level  Diagram  in  Figure  42.  The  charge  distribu¬ 
tions  for  this  cluster  are  shown  in  Figure  43.  It  will  be 
noted  that  there  is  a  highly  localized  state  about  one  of  the 
phosphorus  atoms  (Figure  43c, d).  This  is  consistent  with 
assigning  the  unionized  donor  state  to  a  level  deep  in  the 
energy  gap.  It  also  supports  the  previous  argument  that  there 
is  a  substantial  phosphorus  -phosphorus  interaction  which 
strongly  alters  the  donor  levels. 
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Figure  l  Ethane  Molecule* 
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Figure  2 


Neo-pentane  Molecule  and  CcC' 
Cluster.  Shaded  atoms  are°H1?n 
neo-pentane. 
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Figure  8C  Excess  Spin  Distribution  Figure  8D  Charge  Distribution, 

on  Atoms,  Small  Cluster;  Impurity  Acceptor  State 

B  Impurity.  Small  Cluster;  B  Inpurity 


16S1J4  “EXCITAriON  TROM-LEVEC  -5ir 

in  central 


c 
o 

lO  C/3  C  4J 
4  CQ  O  rl 

u  u  ta 

2  4-10  0 
4  O  XI  CV 


)0 

-4 


n^'OCNCfVJ^N 
f\j  m 


r 


0  0  0  0* 

o  *  o  o  o 

o  O  O  O  0 

O  0  0  O  0 

O  O  O  0  0 

o  o  o  o  o 

o  o  o  o  o 

0  0  0  0  0 

O  O  0  o  o 

0  0  0  0  0 


I 

0  O  O  0  o 
0  o  o  o  o 
ft  o  0  o  o 
O  O  O  O  0 
0  O  O  O  0 
0  0  0  0  0 
0*000 
0  o  o  o  o 
O  O  0  o  o 
*  O  0  w-  * 


i 

^  ^  « 

_ _ i 

i 

(M  n  4  m 

>0 

l 

| 

WWW 

W  w 

1 

1 

i 

!  ! 

<  n  co  aa  < 

< 

(\jnoN!M 

< 

O'  N  N  X  X 

j 

1  I 

N  IT)  O  rO  vfl 

in 

no--* 

in  n 

fU  1 

• 

•  •  • 

1  •  • 

•  t 

rv  rvj  n  m  4 

LO 

I  <  I 

i 

ii  i 

i 

i 

jl 

1  ! 

33  2)  X  <  D 

o  ui  ►-  r-  c 

< 

i 

i 

0>NNX  X 

| 

in 

1 

1 

O'  x  rvj  O'  x 

i 

to 

nonipn 

rv 

•  •  • 

•  • 

• 

rvj  rv  rn  n  4 

in  | 

i  a  i 

1 

i  a 

i*  1 

»  !  !  | 

4  O'  rv  O' 

O'  x> 

1 

:  i 

o 

n  x  x 

rv  o'  <■  o 

VC  4  O' 

c 

rv 

• 

•  •  • 

•  •  • 

•  • 

— •  rv  rvj  rn  rn  <r  <• 

in 

1 

i  1  i 

I 

i  i  i 

■i  ! 

i  i 

i 

T  <  X  T>  < 

«441 

■j  (\jCf  j-  c 

<C  X 

m 

©  e 

ir-NCCJHP  | 

!  l  1 

rvj  vc 

^  7  x  n  x  n 

•  • 

•  •  • 

•  •  • 

•  • 

r\j  no  n 

*5  4  4  in  1 

i  i 

i  1  l 

1  1  1 

- 1 

T  < 

ft  <  < 

r  <  v 

i 

< 

-">  <*»  (M  1/1  — 

d  o'  ^  in 

OOO'SN 

iH  j 

o  in  in  in 

1 

in  4  in  o  1- 

m  v£ 

- 

^  X  4  O'M 

-  <  r- 

•  • 

•  •  •  i 

•  •  • 

• 

-2 

-3 

-3 

o  4  4 

i  i 

in 

i 

X  < 

D  <  <  [ 

I  cc  < 

D  < 

rv  c 

xl 

■»  r*fc 

n  m 

CCICNCJ 

d  v£  in 

n  in 

<t  in  so  ©  in  *- 

<  in  o 

O  O' 

4  X 

/  cmr  4s 

^  in 

•  • 

•  •  • 

»  •  • 

•  • 

r-  cc  O'  o 


rvj 


□0<2}3)lSa?COD 
inrorr'^r-.acin© 
4  ^  n  n  n  m  f\i 


4rnrvrv*-c\jin.n 

«X  o  h-  i-i  lO  ©  vTI  AJ 


O'  O  O  H  rH  (\J  (>J  ,”1 

|  <—•  W  W  W 

I  I  I  I  I  I  I 


T  <<<<<<<< 

in  sj  4’rnrr)rr)'\jojr\j 

o>  in  ji  o>  m  4  -r  in  o 

r- 

•  •••••••• 

x  f> 

I  |  M  ^  t-4  M  *-t  M 

i  I  I  I  I  I  I  I 


oT'LT^r^nor^  *  O' 
in444nnmr\jrvj** 


rvjooT'Ttrccr^o'  c 
X  rj  O'  —  r*.  ^  in  ©  in  m 

X  O'  O  ©  <r  w  f\j  <V  m 

I  I  I  i— •  ^  •-*  *-•  w  ^ 

I  I  I  I  I  I  il  I 


<r'rfr<<r<r<r<j<r 
fvx  -J  "  x  >J  cr  r-  0s 
m  <  <  nrnfVM^ 


n  it  (v  o  c'-cmcT'©  — 
®  ^  cr  ^  x  huc  o  sf 

T  O'  a  c  t>  —  rv  n. 

I  |  ,| 

. . 


<<TTO<IT<I<X 

^0'r)cr^r>0'«cnx 

n  4  4  ^  nri  MMfv^ 


I  rv  r- 

i  lc 


O'f^lx-ornirtn^l- 

X  J‘0'fVX^vC'r^lc< 

•  •  . . .  • 

x  O'  o  ©»-.►-  rvj  f\j  rn 
i  I  ;!*■*•-*►*»■**-•►-«*■* 
I  1  SI  I  I  II  I 


lr<<<<TrD<r(r< 
rv©ac4C'rrvo>cn<x 
in  in  4  4‘Jr^rnrvrvrv^ 

i  r  r  r 

oino'or^xa/o.ncttn 

mi74C7rVO'4f^rVf0  4 

!  •  •  •!•  •  •  •  •  •  i  •  • 

it  r  O'  O'  c  o  ^  ^  i\j  ^  n 

>1  I  I  I  t-t  r-i  *-»  r* 

I  I  il  I  I  I  I 


■NNtt  NvDlOO 

1  I  •— 


r 


o  o 
o  o 
o  o 
o  o 
o  0 
0  0 
o  o 
o  o 
0  o 
0  0 


o  o 
0  0 
o  o 
o  o 
0  0 
o  o 
o  o 
0  o 
o  o 
0  0 


0  O  O  0  0  0  0 

o  o  o  o  o  o  o 
0  0  0  0  0  0  0 
0  0  ft  0  0  O  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
0  0  0  0  0  0  0 
o  o  O  0  0  0  0 
0  0  o  0  0  0  0 
ooooooo 


L. 

+J 

u. 

Lt 

u. 

-J 

u. 


1 

; 

j 

, 

r-^ 

| 

1 

j 

< 

< 

! 

! 

< 

1 

1 

i  r>i 

I  "" 

1 

j 

i  r 

I 

' 

rn 

n* 

i  O 

in 

• 

j 

,  • 

< 

I  < 

, 

‘O 

w 

f— 

| 

1 

1 

j 

1 

1 

I 

j 

03 

T 

• 

CO 

1 

i 

vO 

n 

; 

I 

►— 

I 

i 

| 

i 

1  O. 

| 

io 

| 

i  0 

in 

j 

• 

• 

• 

i 

< 

1  o 

i 

X 

*-• 

i  r 

r* 

j 

1 

i  1 

i 

i 

i 

T 

1 

i  < 

< 

< 

in 

fXJ 

© 

•0 

*"* 

1  «-• 

! 

o 

rv  nj 

> 

r-  rv 

in 

i 

• 

• 

■  • 

• 

! 

< 

* 

fs. 

r 

r— • 

*-* 

W 

1 

1 

.1 

1 

< 

r  oj 

1 

IT 

<  r 

in 

-VC 

<  n 

w 

r* 

'C 

mn 

O' 

n  < 

CC 

r-  rv 

in 

O'  © 

• 

• 

• 

l  • 

•  • 

< 

<C  n- 

JX 

©  -* 

r-» 

p* 

rv  rv 

1 

1 

1 

’’ 

l  i 

T 

< 

<  T 

ir 

X  T 

< 

>  X 

p 

<  rv 

vO 

c 

c  0 

~  . 

r>  so 

X 

X  n  n-» 

it 

O'  c 

• 

• 

•  • 

• 

•  • 

< 

ON  CO 

CC 

o  •— 

rv  rv 

1 

.  1 

1  1 

1 

i  I 

< 

X  < 

< 

<  < 

< 

7  X 

in 

n  rv 

< 

1 

rv  rv  ^ 

4 

©  N* 

O' 

.  X 

o  rv  in 

©  © 

• 

•  1 

•  • 

• 

•  • 

*  ! 

«  N  T 

Ti 

W  M 

rv  rv 

i 

1 

1  1 

1 

>o 

t  (V 

C 

r>  R 

w 

W 

■»  w 

- 

w  w 

0 

* 

0 

0  0 

0 

0 

0  0 

0 

0  < 

0  0 

o  o 

0 

0  s 

0 

0  0 

0 

0  a 

0  0 

0  0 

0 

0  4 

0  0 

0  0 

0 

0  0 

• 

i 

0  0 

0 

0  0 

*  * 

0  0 

0 

0  0 

*  f1 

0  0 

0 

0 

0 

! 

0  0 

!  O 


© 

o 


ID 

I 


•— «  j 


PC  ~ 


IfJ  M 
O'  ►* 
I  •  M 

w  w 

IM 

•  H 


127 


o 

in 


o 

O  :■ 

•  I 

o 

<V  I 


in 

rv 

I 


TERM  rTEVEC~=  =Fi  59  EWEHGf^GAT 


si,  0.006  Si  0.022 


128 


Figure  ioa  Total  Charge  Deficiency  Figure  iob  Excess  Spin  Distribution 

on  Atoms;  Large  Cluster,  on  Atoms;  Large  Cluster, 

B  Impurity  In  Center.  B  Impurity  in  Center. 
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Figure  ioc  Charge  Distribution  of 
Impurity  Acceptor  State; 
Large  Cluster,  B  Impurity 
in  Center. 
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Figure  12c  Charge  Distribution, 

Impurity  Acceptor  State 
B  Impurity  Off-Center. 
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Figure  14A  Total  Charge  Excess  Figure  mb  Total  Charge  Deficiency 

on  Atoms,  Small  Cluster,  in  Bonds,  Small  Cluster, 

F  Impurity.  P  impurity. 
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Figure  16c  Charge  Distribution, 
Impurity  Donor  State, 
Large  Cluster,  P  Impurity 


Figure  17A  Atoms  Surrounding  a  Vacancy. 


Figure  17B  Coordinate  System  for  Tetragonal 
Distortion  of  Atoms  Surrounding 
Vacancy. 
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Figure  20  Scaled  Energy  Level  Diagram  of  Vacancy. 
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Figure  21E  Distribution  of  Charge  Figure  21F  Distribution  of  Charge 

in  Highest  Valence  State:  in  Lowest  Conduction  State: 

Vacancy  Cluster.  Vacancy  Cluster. 
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Figure  23  Scaled  Energy  Level  Diagram:  Vacancy  Plus 
Boron  Impurity:  Nearest  Neighbors. 
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Figure  24c  Excess  Spin  on  Atoms; 

Vacancy  +  B  Impurity; 
Nearest  Neighbors. 
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Impurity:  Next  Nearest  Neighbors. 
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Figure  27c  Excess  Spin  on  Atoms: 

Vacancy  +  B  Impurity: 
Next  Nearest  Neighbors. 
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Figure  29  Scaled  Energy  Level  Diagram:  Vacancy  Plus 

Phosphorus  Impurity:  Next  Nearest  Neighbors. 
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Figure  34A  Charge  Distribution  Figure  34b  Bond  Populations:  Silicon 

on  Atoms:  Silicon  Cluster  Cluster  with  2  Boron 

with  2  Boron  Impurities:  Impurities:  Next  Nearest 

Next  Nearest  Neighbors.  Neighbors. 


SI  0 >027 


163 


of  silicon  cluster 
with  2  phosphorus 
impurities  next 
nearest  neighbors. 


w 

t4 

U* 


a 

ON 

O 

ol 

• 

• 

•  • 

1 

m 

1 

in 

1 

•>  •> 

•» 

an  n- 

ro 

ro 

—*  \r 

o  • 

in 

•  • 

• 

• 

<M  <\i 

1  1 

co 

1 

in 

1 

* 

«►  •> 

•»  » 

C  If 

<r  o 

no  ro 

cc  o  cr  <\»  nC 


^  CO  f\j 

1  1  1 

ro  ^ 

1  '  ' 

*3“  in 

1  i 

*  •  • 

I 

•  • 

*  •> 

iin  ^ 

o  *c 

—  o 

la;  ro  a- 

^3-  O 

in 

i  •  •  • 

•  • 

•  • 

•  no  ro 

1  1  1 

m  sj- 

1  1 

•3*  m 

1  1 

-  •* 

>c  n 

ro  n* 

CC  o 

O'  ^  CO 

•O  ^ 

cd  a c 

•  •  • 

•  • 

•  • 

•— *  fV  <v» 

1  1  1 

ro 

1  1 

>r  in 

l  1 

»  •> 

•  •> 

O  vf  vD 

r-  a :  ro 

O  ro  r— 

^  r-  O 

n-  in , 

•  •  • 

•  •  • 

•  • 

no  r\i  ro 

1  1  1 

ro  ro  nJ- 

1  1  1 

>3  IT, 

1  1 

_ 

^  in  >r 

o 

sj-  03 

WWW 

*  *  it 

*  *  * 

* 

9  if-  it 

*  *  * 

*  *  H 

it  it  it  ; 

>  4>  H 

*  *  #  i 

*  ♦  *  it\i 

cc 

“II 


I  I  I 


o 

o 


**•■»•* 
it  *  *  * 
it  *  *  i* 


*  * 
#  * 
#  * 


«  lit 

*  I* 

#  l* 

*|* 


*  »•*> 
-rf  |4* 

it 


*  * 

*  * 
*  * 
it  * 

*  it 
it  # 

*  it 

*  it 

*  * 


• 

• 

i  • 

► 

tM 

jro 

s0 

1  >o 

h- 

■  ^ 

o 

• 

• 

• 

• 

O' 

ro 

ro 

\  1 

1 

i 

i 

|  1 

i 

1 . 

i 

» 

•» 

0?  O  >1 

I^CCH 

O 

m  03  O 

o.'  ^ 

n- 

ro 

•  •  • 

•  •  • 

• 

• 

9 

9 

0 

o  ^  ^ 

ro 

rO 

1  7 

i  7.7 

s-4 

1 

i 

.  .. 

•» 

*c  ^ 

r-  r-  ro 

cc 

ro  cn  o 

r*-  ro  ,m 

m 

•  •  • 

c  a 

o  —  ^ 

ro 

ro 

.  .  r 

i  i  i 

Y* 

»  • 

O  n- 

m  o 

r-  C'  C 

co  r- 

cc 

cr 

^  O'  ro 

C*  ro,  x 

cr  r- 

m  i 

•  • 

•  •  • 

«  • 

•  ! 

to  cr 

I  | 

c*  o  a 

rt  M 

ro  ro 

ro  f 

''ii 

.J 

1  1  1 

1  1 

1 

ro  Nf  r<v 

fo  co  sO 

WWW 

— 

— 

| 

*  * 

it  it  it 

it  it  it 

*  * 

*  u 

it  *  1  *  *  *  i 

'l  4:*  ■«•  i 

*  * 

* 

■tt  a  ut  it  -if  ns  ^  it  it  a- 

it  it  i 

»**!»**■' 

* » 

*  * 

it  it  \ 

#•  1 ir  <-  1 

*  it 

«■  1  + 

i'-  *  i  i:-  it  it  <»•  it  *  ■ 

:t  it 

*  :i 

*  *  !* 

*  **  i-u 


Us 


*  Uf 

*  i* 

*  i* 


*  # 


i!  I* 
Hi  i« 

*  i* 


it  :■* 
*  * 
it  lit 


i  r\j 


rO  vT 


* 

* 

!«• 

i« 

I* 

i* 

# 


ro!^* 

Ic* 


<r 

I 


•sf  isj  sj-  co  jrr  ^ 

!  O'  m  r-  j  O'  o 

•  i  •  •  •  i  •  •  • 

•c i*o  r-  r-  |x  co  O' 

•“ l*“ •  r-* 

III  I  I  I  I  I 


o 

o 


^■4 

I 


* 

<r  fr  I1* 

4r  ■{}■  i  -vf 

it  +  *it 
it  I  * 
■it  it  j  it 

it  it  '  it 
<r  i*  {* 
it  it  « 
it  |  * 


164 


*  * 

*  it 

#  * 
#  * 
«  * 

*  it 

#  * 
*  * 
*  * 
*  * 


fV  ro 

I  I 


*  * 
*  * 
*  » 
*  # 
it  it 
•»•  * 
*  # 
*  # 
»  » 
#  # 


iiiiiiiimiiiiinimiiiiiniimiiiiiiiinimmiinminmiuiiiiiiiiiiii 

FCSm.I  level  =  -7.92  ENER GY  GAP  =  0796 


SI3.917 


165 


Figure  36A  Charge  Distribution  Figure  36B  Bond  Population: 

on  Atoms:  Silicon  Cluster  Silicon  Cluster  with 

with  2  P  Impurities:  2  P  Impurities: 

Next  Nearest  Neighbors.  Next  Nearest  Neighbors. 
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Figure  36c  Charge  Distribution  Figure  3 6D  Charge  Distribution  in 

trt  Highest  Populated  Orbital:  Lowest  Unoccupied  Orbital; 

Silicon  Cluster  with  2  P  Silicon  Cluster  with  2  P 

Impurities:  Next  Nearest  Impurities.:  Next  Nearest 

Neighbors.  Neighbors. 
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Atoms:  Silicon  Cluster  with  B  and  P  Impurities: 
B  and  P  Impurities:  Nearest  Neighbors. 
Nearest  Neighbors. 
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Figure  40A  Charge  Distribution  on  Figure  40B  Bond  Populations; 

Atoms:  Silicon  Cluster  Silicon  Cluster  with 

with  B  and  P  rmpurities:  B  and  P  Impurities: 

Next  Nearest  Neighbors.  Next  Nearest  Neighbors. 
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igure  42  Scaled  Energy  Level  Diagram:  Silicon  Cluster  wit 
2  Phosphorus  and  1  Boron  Impurities. 
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Figure  43c  Excess  Spin  Distributions:  Figure  43D  Charge  Distribution  in 

Silicon  Cluster  with  Highest  Occupied  Orbital: 

2  P  and  1  B  Impurities.  Silicon  Cluster  with 

2  P  and  1  B  Impurities. 


VI.  Investigation  of  Recombination  and  Trapping  Processes  in  3As„Se„-2Sb„Se„ 

Amorphous  Films  by  Photoconductivity  Method  (By  S.  S.  Li  andP?  K.  (Thaudhari) 

1.  INTRODUCTION 

The  electrical  and  optical  properties  of  3As2Se3-25b2Se3amorphous 

films  have  recently  been  studied  by  Chaudhari  et.al.1-2^  Possible  band 

structures  were  proposed  by  them  and  detailed  information  on  the  basic 

parameters  such  as  mobility  gap,  electron  mobility,  absorption  coefficient, 

optical  gap  and  index  of  refraction  for  these  films  have  been  reported 

in  their  papers.^-  ^ 

3-4) 

Weiser  et.al.  have  recently  reported  the  experimental  results 
on  the  transport  and  photoconductivity  in  the  2As2Te3-As2Se3  amorphous 
system,  and  analyzed  their  photoconductivity  data  in  terms  of  "recombination 
edges"  in  the  CFO  model  for  amorphous  chalcogenides.  “*) 

In  this  paper  we  report  the  results  of  our  photoconductivity 
experiment  on  3As2Se3~2Sb2Se3  amorphous  films  in  the  temprature  range 
from  278°K  to  308°K.  Experimental  results  on  these  films  suggest  to  us 
that  their  photoconductivity  behavior  is  analogous  to  high  resistivity 
semiconductors  in  which  the  recombination  and  trapping  processes  can  be 
considered  as  taking  place  through  an  effective  trap  center.  A  simple 
model  is  established  to  account  for  the  recombination  and  trapping 
processes  of  the  photo-injected  carriers  in  these  films. 

2.  THEORY 

Results  of  the  photoconductivity  experiments  on  the  specimens  of 
3As2Se3~2Sb2Se3amorphous  films  indicate  that  high  injection  conditions 
prevail  over  the  light  intensity  ranges  used  in  the  present  experiment. 

This  means  that  the  light  intensity  is  strong  enough  so  that  the  condition 
An>>no  and  Ap»pQ  is  satisfied.  The  observed  dependence,  of  the.  photo¬ 
current  on  light  intensity,  Iq,  could  be  well  interpreted  by  using  a 
recombination  and  trapping  model  which  takes  into  account  the  trapping 
of  carriers  through  an  effective  trap  center  and  the  dependence  of  the 
carrier  lifetime  on  the  injection. 

A  generalized  express  ion  for  the  photoconductivity  in  these  films  may 
be  derived  under  the  assumptions  that: 

(1)  The  effect  of  trapping  through  an  effective  single  trap  center  can  be 

expressed  by 
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Ap  =  r/Sn 


(1) 


where 


r  -  J-  <  x 

n 


for  electron  trapping, 
for  hole  trapping 


and 


(2)  The  dependence  of  excess  carrier  lifetime  on  injected  carrier  density 
(or  light  intensity)  may  be  expressed  in  the  form: 


T  =  «An^ 
n 


(2) 


where  <  is  a  constant,  independent  of  An,  and  An  is  the  photoinjected 
excess  electron  density;  8  is  a  parameter  which  may  vary  between  —1  and  0, 
depending  upon  the  injection  range.  The  photocurrent  per  unit  electric 
field  is  given  by 


d 

ip.c.  =  q  1  W(^nAlv4ViP)dy 

o  v 

An 

=  qp  /  °W  (An  +  %  - 2^2 - 

n  b  An 

[2  f  DRdAn] 

o 


(3) 


where  D  =Dn(P+n^)  /  (p+bn)  is  the  generalized  diffusion  coefficient (3-a) 

and. 

An 

R  —  is  the  rate  of  recombination  (3-b) 

n 

Here  AnQ  is  the  excess  electron  density  at  the  illumination  surface  and  d  is 

the  film  thickness;  W  is  the  width  of  the  specimen,  b=  y  Au  .  Note  that  fot 

n  p 

large  injection,  p  =  Ap  and  n  =  An  in  (3-a). 

With  the  help  of  Eqs,  (1)  and  (2),  Eqs.  (3-a)  and  (3-b)  reduce  ^o: 
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and 


21)  r 

D  ^  — 2 — 

v  (r+b) 


_  /sn 
tcAn6 


-1 

K 


(4-a) 


(4-b) 


Substituting  (4-a),  (4-b)  and  (1)  into  Eq.  (3),  we  obtain: 

L/2 


Xpc  ~  ( 


2q 

2 


wn  YD„»:r(b+r)  a-ztf 

j_  -j  ) 


WAn 


(2+C)/2 


(5) 


Eq.  (5)  is  valid  for  different  injection  ranges, so  long  as  the  assumptions 


made  in  Eqs.  (1)  and  (2)  are  valid.  The  values  of  ^Inay  vary  from-  one  injection 
range  to  another  and  may  be  determined  from  the  slope  of  i  versus  I  (light 
intensity)  plot.  For  example,  0  -  0,  if  the  carrier  lifetime  is  constant  and 


independent  of  light  intensity;  this  is  indeed  the  case  in  the  low  and  high 
injection  regions.  Furthermore,  if  the  effect  of  trapping  is  negligible,  then 
F  in  Eq.  (5)  is  equal  to  unity. 

The  injected  electron  density,  AnQ,  at  the  illuminated  surface  is  directly 
related  to  the  incident  photon  flux  density,  Qo,  and  may  be  determined  from 
the  boundary  condition.  The  result  for  the  present  case  is  given  by ^ 


<  T)  1  *  “ 
Q  =  8  An  +  2  1  "TTT  Y 1  An 
o  o  o  *  K  (2-J»)  )  c 


2-0 


(  6) 

where  8q  is  the  surface  recombination  velocity  at  the  illuminated  surface. 

Qo  is  the  photon  flux  density  which  is  related  to  the  light  intensity  Iq (watt/cm2 
by  *  tty”  •  D  iS  the  effective  diffusion  constant,  and  is  given  by  Eq.  (4-a)/ 

In  general,  wefound  that  the  dependence  of  the  photocurrent  on  the 
light  intensity  in  semi-insulating  materials  (e.g.,  Cr-doped  GaAs)  could 
be  analyzed  by  means  of  Eqs.  (5)  and  (6),  using  a  monochromatic  light  source. ^ 

3.  EXPERIMENTAL  DETAILS 

Samples  on  Corning  0211  glass  were  employed  for  photoconductivity 
measurements.  An  Aminco  Grating  monochromator,  model  #4-8400, 

was  used  with  a  650.  watt  Tungsten  movie  projector  lamp  focused  on  the 
entrance  of  the  monochromator.  The  sample  was  mounted  inside  the  mono¬ 
chromator  housing  to ovoid  any  background  illumination.  Light  was 
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focused  on  the  back  of  the  semitransparent  substrate  in  order  to  eliminate 
any  contact  effect  (e.g.  photovoltaic  effect).  Ohmic  contacts  weremade 
on  these  films  by  evaporating  a  small  gold  dot  on  top  of  the  film. 

I-V  measurements  on  these  films  indicate  that  good  ohmic  contact  is 
obtained  for  the  test  specimens.  The  d.c.  photoconductivity  was  measured 
by  using  an  HP  413A  micro-micro  ammeter  and  a  Harrison  6112A  digital  power 
supply.  One  volt  bias  was  used  for  most  of  the  data.  The  photo-response 
was  measured  by  using  a  constant  slit  width.  The  films  shoved  considerable 
photoresponse  convering  the  whole  visible  spectral  range  (Fig.  1).  The 
relative  spectral  output  of  the  light  source,  monochromator  mirror,  and 
lens  combination  was  measured  by  an  Optics  Technology  model  #615  power 
meter. 

Other  measurements  such  as  conductivity,  absortption  coefficients  in 

these  films  were  reported  in  references  1  and  2,  and  will  not  be  repeated 
here. 

4.  RESULTS  AND  DISCUSSIONS 

The  specimens  used  in  the  present  photoconductivity  experiment  were 
identical  to  those  used  in  the  electrical  conductivity  measurements 
reported  by  Chaudhari  e.al.2)  The  film  has  a  dark  resistivity  of  about 
10  ft-cm  and  the  electron  mobility  estimated  from  optical  measurements , 1_2) 
is  less  than  5  cm  /V-sec.,  at  300°K.  The  present  photoconductivity 
experiment  is  performed  under  the  condition  that  the  photoinjected  excess 
electron  density  is  much  greater  than  the  dark  electron  density  (i.e. 

An>no)*  i-nvest igate  the  recombination  and  trapping  processes  of  the 
photoinjected  carriers  in  these  films,  we  first  measure  the  spectral 
response  of  the  photoconductivity.  The  result  is  shown  in  Fig.  1  for 
a  2.23y  thick  3As2Se^-2Sb2Se^  film.  Here  the  photocurrent  is  normalized 
to  the  per  unit,  incident  power  as  me cured  with  the  optical  power  meter. 

There  is  no  evidence  for  a  wavelength  at  which  the  response  vanishes 
completely*.  In  such  a  case  the  band  gap  is  defined  as  the  energy  at 
which  the  response  has  fallen  to  half  of  its  peak  value.  This  yields 
Ep^  ~  1*  70  (+0..02)  eV,  as  the  photoconductivity  gap,  which  is  referred  to 
as  an  intrinsic  gap  for  the  following  reasons;  (1)  the  absorption  due  to 
localized  states  would  not  give  rise  to  photoconductivity;  and  (2)  the 
recombination  time  for  transitions  from  the  tail  of  the  valence  band  to 
the  nonlocal ized  states  in  the  conduction  band  might  be  considerably 
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less  than  that  for  the  intrinsic  process.  It  is  clear  ,'rom  Fig.  1  that 

thin  films  of  3As2Se^-2Sb2Se2  are  sensitive  photoconductors  over  the 

range  of  visible  light.  The  fall-off  in  the  photoconductivity  response 

at  higher  energy  is  due  to  the  recombination  of  excess  carriers  at  the  surface. 

Optical  absorption  data  reported  in  reference  1  reveals  that  intrinsic 

absorption  coefficient  is  in  excess  of  lO^cm  '''for  photon  energy 

greater  than  1.67eV  (or  A<0.73  pm).^ 

In  investigating  the  recombination  mechanism,  we  measure  the  photo¬ 
current  versus  light  intensity,  I  ,  by  choosing  monochromatic  light 

o  O 

with  A  =  5850A  for  several  3As2Se^-23b2Se2  amorphous  films.  The  result 
for  a  1.35p  thick  film  is  displayed  in  Fig.  2  for  different  temperatures. 

Two  distinct  injectionregions  were  observed.  At  the  intermediate  light 
levels,  is  di.reocly  proportional  to  Iq.  This  result  may  be 

interpreted  in  terms  of  the  theoretical  expression  given  in  the  previous 
section. 

(a)  Intermediate  light  level: 

In  this  injection  range,  we  assume  that  the  effect  of  electron  trapping 
in  the  trap  centers  is  important  and  thus  the  relation  between  An  andAp  is 
described  by  Eq.  (1),  and  the  electron  lifetime  depends  on  the  light  intensity 
in  the  form  that 

1  -  (j.e.  with  (3  =  -1/3  in  Eq .  (2)) 

n 

(7) 

The  dependence  of  the  carrier  lifetime  on  the  injection  is  due  to  the 
fact  that  in  the  intermediate  injection  level  charge  fluctuation  in  the 
trap  centers  plays  an  important  role  in  controlling  the  carrier  recombination 
lifetimes.  The  carrier  lifetime  is  normally  decreased  with  increasing  light 
intensity.  If  the  injection  is  very  high  (strong  light  intensity),  the 
trap  centers  are  saturated  by  the  excess  carriers,  ind  in  this  case  the  carrier 
lifetime  reaches  a  constant  value  at  high  injection  level. 

If  the  inci ’ent  radiation  is  uniformly  absorbed  in  the  films,  and  thus 
the  bulk  term  in  Eq.(6)  is  the  dominant  term  that  contributes  to  the 
photocurrent,  then  the  result  for  photocurrent  in  the  intermediate  injection 
range  can  be  deduced  from  Eqs.  (5-6). 
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2(f) 

7k 


(8) 


which  shows  that  i  varies  in  proportion  with  Q  '  in  this  injection  range, 
pc  o 

Since  the  photon  flux  density  Q  is  linearly  related  to  the  light  intensity 

°  ,  0.71 

I  ,  we  conclude  that  i  in  this  range  is  directly  proportional  to  I  ’  , 
o  pc  ° 

which  is  in  excellent  agreement  with  the  observed  result  shown  in  Fig.  2. 

(b)  The  high  light  levels: 

If  the  light  intensity  is  further  increased,  we  see  that  the  slope  of 

i  versus  I  approaches  unity  (see  Fig.  2).  This  means  that  at  high  light 
pc  o 

levels  the  lifetime  of  excess  carriers  becomes  independent  of  light  intensity 
and  the  effect  of  trapping  is  negligible.  The  expression  of  the  photocurrent 
per  unit  electric  field  for  this  case  is  deduced  from  Eqs.  (5-6)  by  setting 
l  and  B=  0,  which  yields: 


i 

pc 


■  q"P  (1+b)  “V,,’ 1  /2‘V '  v(v  v1  /2> 


(9) 


where  D  =  ^n  ig  the  high  injection  ambipolar  difussion  constant. 

Eq.  (9)(Pre&lcts  that  at  high  light  levels,  the  photocurrent  is 
directly  proportional  to  (or  1^),  which  agrees  welll  with  the  observed 
data  shown  in  Fig.  2. 

Summing  up  the  above  analysis,  we  come  to  the  conclusion  that  the 
effect  of  trapping  is  obvious  in  the  intermediate  light  levels, and  the  carrier 
lifetime  is  a  function  of  light  intensity.  However,  at  high 
light  levels,  the  trap  levels  are  saturated  by  the  excess  electrons,  and 
therefore  the  carrier  lifetime  reaches  a  constant  value,  which  is  consistent 
with  our  obsevation.  The  present  results  on  the  recombination  and  trapping 
processes  of  the  photoinjected  carriers  in  these  high  resistivity  amorphous 
films  are  similar  to  that  observed  in  the  high  resistivity  semiconductors 


such  as  Cr-doped  GaAs.^ 

To  determine  the  trap  levels  in  these  films,  we  investigate  the 


temperature  dependence  of  the  photocurrent  at  a  constant  light  intensity. 

Fig.  3  shows  i  versus  1/T;  the  slope  of  this  semi-log  plot  yields  an 
pc 
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effective  trap  level  at  0.3eV  below  the  band  edge.  The  trap-depth  may 

correspond  to  a  maximum  in  the  statistical  distribution  of  discrete  levels. 

These  levels  may  be  due  to  structural  defects,  short  range  disorder 

fluctuations. ^  Results  of  Kolomiets  et  al.^  for  Tl2^e2-As2'Te2  and 

As„Se  glasses  also  show  evidence  of  trap  levels.  In  the  presence  of 

1  J  11) 

such  a  distribution  Rose  has  shown  that  the  photocurrent  varies  as  a 
fraction  power  (i^is  proportional  to  1^,  with  0.5<n<l)  of  the  light 
intensity.  This  conclusion  supports  our  theoretical  model  and  the 
experimental  findings  reported  here. 
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VII. 


DIELECTRIC  PROPERTIES  OF 
LITHIUM-DISILICATE  G L AS S - C E RAM I C S 


By 

W.  D.  Tuohig  and  L.  L.  Hench 


Introduction 

The  importance  of  microstructure  as  a  ' ariable  for  deter¬ 
mining  dielectric  properties  of  solids  has  been  recognized 
for  many  years.  Maxwell^  derived  an  expression  for  the 
dielectric  properties  of  a  layered  body  consisting  of  material 
having  two  different  dielectric  constants.  Wagner^  extended 
this  approach  to  materials  containing  dispersed  spheres  of 
material  A  having  a  dielectric  constant  in  a  continuum 
matrix  B  with  a  dielectric  constant  Kj.  Sillars(3)  extended 
the  model  to  include  dispersed  spheroids  and,  more  recently, 
Fncke  derived  a  relationship  for  ellipsoids.  All  of  these 
expressions  contain  the  conductivity,  dielectric  constant  and 
volume  fraction  of  both  phases  in  describing  the  dielectric 
properties  of  a  two-phase  composite  material. 

Most  commercial  ceramic  dielectrics  contain  a  vitreous 
phase  which  plays  a  key  role  in  determining  the  performance  of 
the  material,  both  from  an  electrical  and  a  mechanical  stand¬ 
point.  The  presence  of  a  glassy  matrix  in  porcelain  insulators 
renders  them  both  strong  and  impervious  to  moisture.  Because 
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of  high  dielectric  losses  usually  associated  with  glasses,  the 
extent  of  vitreous  phase  must  be  minimized  for  high  tempera¬ 
ture  dielectrics.  Therefore,  both  the  composition  and  the 
amount  of  glassy  phase  are  of  the  utmost  importance  in  deter¬ 
mining  the  suitability  of  a  ceramic  material  for  electronic 
applications.  A  general  review^^  of  mi  c  r  o  s  t  r  uc  t  ural  effects 
in  ceramic  dielectrics  also  points  out  that  porosity,  grain 
boundary  structures  and  other  mi c r o s t r uc t ur al  features  are 
important  parameters  for  ceramic  dielectrics.  This  study  also 
points  out  that  little  systematic  information  is  available  for 
characterized  mi c r o s t r uc t ur e s . 

Experimental  determination  of  the  dielectric  properties 
of  heterogeneous  materials  was  reported  by  Buchner^^  for 

TiO  in  various  matrices.  The  effects  of  porosity  on  dielec- 

L  (7) 

trie  properties  has  been  studied  by  Economos. 

The  dielectric  properties  of  glasses  has  been  the  subject 
of  study  by  a  number  of  workers.  Both  Owen  ’  and  Isandv 
explained  their  results  in  terms  of  a  type  of  heterogeneous 
dielectric  model.  Cha r le s ' ^ ^ ^ ^  study  of  dielectric  properties 
of  lithia-silica  glasses  indicated  the  importance  of  phase 
separated  mic ro s t r uc t u re s  upon  the  interpretation  of  dielec¬ 
tric  properties.  More  recently,  Kinser  and  Hench  have 

shown  the  importance  of  thermal  history  and  resulting  crystal¬ 
line  microstructure  upon  the  dielectric  properties  of  glasses. 
Their  study  further  indicated  that  the  metastable  crystalline 
phase  precipitated  as  the  result  of  isothermal  heat  treatment. 
It  was  established  that  this  metastable  phase  was  a  precursor 
to  equilibrium  phases  during  controlled  devitrification  in  tiie 
lithia-silica  system. 

Recent  studies  have  shown  that  a  controlled  volume  frac¬ 
tion  of  L^S^O,.  crystals  can  be  grown  in  a  glass  of  the  same 

composition  with  appropriate  nucleation  and  crystallization 

(12) 

treatments . 

It  is  the  objective  of  the  present  study  to  establish  the 
dielectric  behavior  of  the  33  mole%  R^O-b?  mole%  SiC^  glass 
as  a  function  of  the  volume  fraction  of  LigS^O,.  present. 
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Experimental  Procedure 


Specimens  for  this  study  were  prepared  by  melting  at 
1350°C  in  a  covered  Pt  crucible  5pm  Min-U-Sil  (Pennsylvania 
Glass  Sand  Corporation)  and  Li?CO^  in  amounts  to  yield  a 
33  mole%  Li20-Si02  glass  composition.  Initial  melting  was 
for  a  period  of  18  hours;  the  material  was  then  cast,  crushed 
and  re melted  for  4  hours  at  1350°C  prior  to  casting  disc¬ 
shaped  specimens  1.8  cm  in  diameter  by  approximately  0.5  cm 
thick.  Specimens  were  immediately  transferred  to  an  annealing 
oven  at  300°C  to  relieve  quenching  stresses. 

Vapor  deposited  gold  electrodes  in  double  guard  ring  con 
figuration  were  employed.  A  sample  chamber  allowed  measure¬ 
ments  to  be  made  under  a  lp  vacuum  at  controlled  temperatures. 

The ' apparatus ,  which  includes  Wayne  Kerr  B-221  and  B-601 

(11) 

bridges,  has  been  described  by  Kinser  and  Hench. 

Thermal  treatments  consisted  of  nucleation  at  475°C  for 
24  hours  with  subsequent  crystallization  at  575  C  for  times 
of  25  to  85  minutes.  The  respective  volume  fractions  of 

Li  Si  0  were  10,  25,  50  and  70%  ±  5%  as  inferred  from  the 

225  (12) 

data  of  Freiman  and  Hench. 

Scanning  microscopy  was  undertaken  on  fracture  surfaces 
etched  with  dilute  HF .  The  specimens  were  thoroughly  rinsed 
and  dried  and  plated  with  gold-palladium  by  vacuum  deposition. 

Results 

Figure  1  shows  the  frequency  dependence  of  the  me  as t  red 
ac  conductivity.  a  is  almost  independent  of  frequency  over 
the  audio  range  for  the  samples  with  10%  and  25%  crystals 
present.  However,  the  dependence  becomes  more  pronounced  with 
increasing  crystallinity. 

Figure  2  shows  the  dependence  of  conductivity,  AC  and  DC 
on  volume  fraction  of  Li2Si20^  crystals.  A  major  change  in 
dependence  occurs  when  =  .25  due  to  alkali  depletion  of  the 

glass  . 

Figure  3  is  a  plot  of  conductivity  at  1 0,000  Hz  vs.  1/T. 
Data  for  the  as-cast,  10%,  25%  and  50%  samples  exhibit 
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identical  slopes,  implying  that  there  is  no  change  in  activa¬ 
tion  energy.  The  70%  specimen  exhibits  a  somewhat  lower 
activation  energy.  a  decreases  continuously  with  V  as  ex¬ 
pected. 

Figure  4  shows  dielectric  constant  K  1  vs.  frequency  for 
the  few  volume  fractions  of  crystalline  phase  at  100°C.  K' 
decreases  as  Vy  increases  and,  for  Vv  =  0.70,  there  is  almost 
no  change  in  K'  with  frequency  over  the  range. 

Figure  5  shows  the  dielectric  constant  vs.  V  crystalline 
phase  at  500  and  20,000  Hz  measured  at  100°C. 

Figure  6  is  a  plot  of  cDC  vs.  1/T.  Slopes  of  as-cast 
10%,  25%  and  50%  specimens  are  identical,  implying  no  change 
in  activation  energy.  A  continuous  decrease  in  o  occurs  as 
the  concentration  of  free  carriers  (alkali  ions  in  vitreous 
network)  decreases.  The  70%  specimen  exhibits  higher  activa¬ 
tion  energy  for  DC  conduction. 

Figure  7  shows  the  loss  tangent  vs.  frequency  at  a 
measuring  temperature  of  100°C.  As-cast  specimens  show  usual 
low  frequency  dispersion.  The  loss  peak  observed  after  24 
hours  at  475°C  implies  the  presence  of  Li20-Si02  as  reported 
by  Kinser  and  Hench.  With  the  onset  of  equilibrium  crystal¬ 
lization,  loss  peaks  are  no  longer  present,  and  the  loss  factor 
decreases  markedly  with  increasing  volume  fraction  of  crystal¬ 
line  phase. 

Disc  uss ion 

Scanning  electron  micrographs  of  etched  fracture  surfaces 

are  shown  in  Figure  8.  The  formation  of  rod-like  crystals 

reported  by  Freiman  and  Dench  is  evident.  The  V  =  0.50  indi- 

v 

cates  that  crystal  growth  occurs  to  impingement.  It  appears 
that  residual  matrices  apparently  crystallize  in  a  dendritic 
fashion.  This  secondary  crystallization  apparently  produces 
.smaller  crystallites  within  the  interstices  of  the  impinged 
spherulites.  This  microstructure  is  no  longer  evident  in  the 
Vv  =  0.70  specimen,  which  has  an  appearance  of  an  intergranular 
fracture. 
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The  dielectric  properties  appear  to  reflect  the  continuous 
or  matrix  phase.  The  AC  and  DC  conductivity  of  specimens  in 
which  the  major  phase  is  vitreous  are  similar  in  terms  of 
magnitude  and  activation  energy.  As  L^S^O,.  becomes  the 
dominant  constituent,  the  activation  energy  is  clearly  altered. 
The  dielectric  constant  decreases  as  the  volume  fracture  of 
crystalline  phase  increases,  as  does  the  loss  factor.  The 
appearance  of  the  relaxation  peaks  at  2,000  Hz  after  nuclea- 
tion  and  its  subsequent  disappearance  is  in  agreement  with  the 
findings  of  Kinser,  and  with  unpublished  work  of  the  present 
authors . 
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VIII.  Discussion 


The  research  of  this  tenth  semiannual  period  has  provided  technical 
findings  that  bear  on  problems  of  interest  co  the  Department  of  Defense. 

The  performance  of  transistors  plays  a  central  role  in  determining 
the  operation  of  many  electronic  systems  used  in  military  applications. 

For  such  applications  as  detection  and  amplification,  for  example,  the 
gain  and  the  cutoff  frequency  are  especially  important  parameters. 

In  Section  III  of  this  report  we  demonstrate  that  the  bias  voltage  at 
the  collector  terminals  can  control  both  gain  and  cutoff  frequency  and 
determine  the  device  behavior  at  high  currents.  The  study  reported 
here  furnishes  an  understanding  that  may  enable  improved  design  of  devices 
and  aid  in  the  computer-aided  design  of  transistor  circuits. 

Section  IV  extends  the  present  understanding  of  the  quantum  yield 
of  p-i-n  photodiodes.  Such  devices  find  use  as  pbotodetectors  in 
application  of  interest  to  the  D.O.D. 

In  ...  .-fion  V  a  method  is  described  for  the  computer  calculation 
of  the  properties  of  silicon  containing  certain  types  of  defects. 

Among  these  defects  is  the  impurity-vacancy  complex  which  can  commonly 
result  from  radiation  diimage.  Thus  the  study  pertains  to  the  design  of 
radiation  hardened  systems. 

The  last  two  sections  describe  properties  of  amorphous  semi¬ 
conductors  and  glass-ceramics  potentially  critical  in  the  application 
of  these  materials  to  devices  and  military  systems.  Because  of  the 
hardness  they  exhibit  to  irradiation,  amorphous  semiconductors  are 
of  special  interest. 

This  is  the  final  report  in  the  series  describing  research  supported 
by  ARPA  Contract  F  19628-68-C-0058 .  A  major  part  of  the  research  program 
will  continue  under  the  support  of  ARPA  Contract  F  19628-72-C-0368. 
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